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ABSTRACT 

Gater, Roger Allen. Ph.D., Purdue University, January 
1969. A Fundamental Investigation of the Phenomena That 
Characterize Liquid-Film Cooling. Major Professor: 

Dr. Mel R. L'Ecuyer. 

An experimental investigation was conducted to deter- 
mine the net rate of mass transfer (evaporation plus en- 
trainment) from a thin liquid film to a proximate gas stream. 
Three different flow situations were investigated: (a) con- 

stant velocity, cold gas flow; (b) constant velocity, heated 
gas flow; and (c) strongly accelerated, heated gas flow. 

The experiments were conducted in a horizontal test section 
by passing an air stream over a liquid film formed on the 
lower surface. The liquid film was 2-in. wide and was ter- 
minated at a length of 10 in. by a knife-edge capture slot. 
The flow parameters for the experiments with a constant 
velocity gas stream were: gas stream temperatures of 40F, 

400F, and 600F; gas stream static pressures of 75 and 150 
psia; gas stream velocities of 40 to 400 fps; and liquid 
coolants of methanol, butanol, water and RP-1 (a hydro- 
carbon fuel). The flow parameters for the experiments with 
an accelerating gas stream were: stagnation temperature of 

400 F ; stagnation pressures of 50, 75, 100, 125 and 150 psia; 
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variation in gas stream velocity of 200 to 660 fps over the 
liquid film; and liquid coolants of methanol and water. 

The measured rate of mass transfer with a heated gas 
stream was found to be as much as 10 times greater than 
the rate of mass transfer predicted from the simple (or 
classical) theory for mass transfer. In most cases, the 
mass transfer due to the entrainment of bulk liquid into 
the gas stream was the most significant contribution to the 
overall mass transfer. In those cases where evaporative 
mass transfer was important, the "effective" roughness of 
the surface of the liquid film was found to increase the 
rate of evaporation by as much as 80 percent over that pre- 
dicted by the simple theory. 

All of the data for mass transfer (approximately 600 
datum points) was correlated in terms of relatively simple 
dimensional parameters, with one parameter characterizing 
the phenomenon of entrainment and a second parameter char- 
acterizing the phenomenon of film surface roughness. Of 
particular significance is the fact that a fifteenfold 
variation in the liquid viscosity was not found to have any 
noticeable effect on the rate of entrainment. 

An analysis was made to extend the results of the 
present investigation to the more general case of an arbi- 
trary film coolea length. The utility of the analysis was 
investigated by analyzing the data reported in three primary 
references. The agreement between the reported data and the 
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results predicted from the present theory was good, with the 
agreement being to within the accuracy of the experimental 
data in most cases. 

Photographs were taken of the surface of the liquid 
film so that the gross characteristics of the film surface 
could be studied as a function of the pertinent flow vari- 
ables for the liquid and gas phases. 

Measurements were made of the maximum temperature that 
a liquid film obtains for the subject flow situation. A 
theory was developed for the prediction of the maximum liq- 
uid temperature that was in good agreement with the experi- 
mental data. 

An analysis was made of the heat transfer characteris- 
tics for the gas-vapor cooled region immediately downstream 
of the liquid film. The analysis allowed for compressible 
gas flow in the presence of a favorable pressure gradient. 
Consideration of the limited experimental data that are 
reported in the literature for the subject problem indicated 
that the analysis was promising. 
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1. INTRODUCTION 
1.1 Genera] Discussion 

Liquid-film cooling refers to the introduction of a 
thin continuous liquid film onto a given surface for the pur 
pose of protecting that wetted surface from thermal damage 
by a proximate hot gas stream. This technique provides a 
means of protecting, for example, the internal surfaces of a 
rocket motor from the injurious effects of the hot pro- 
pellant gases. The motivation for developing liquid-film 
cooling technology for rocket motors is due to the use of 
higher energy chemical propellants, the trend toward higher 
combustion chamber pressures, and the development of nuclear 
rocket engines. These newer developments are characterized 
by a large increase in the heat flux from the working fluid 
to the walls of the combustion chamber and exhaust nozzle, 
such that the regenerative cooling of those walls becomes 
marginal (l, 2, 3)*. 

In the design of a liquid-film cooled system for a 
rocket motor, the two principal problems to be considered 
are (a) that of determining the rate of coolant consumption 
that will be required to wet, and thus thermally protect, a 


♦Numbers in parentheses refer to REFERENCES. 
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specified internal surface area, and (b) that of determining 
the effect that the cooling process has on the performance 
of the system. In the subject investigation, attention was 
concentrated on the first problem only; specifically, with 
the investigation of the phenomena that characterize the 
convective exchange of heat and mass between a hot gas stream 
and a thin liquid film. Although the principal motivation 
for this research was the interest in liquid-film cooling, 
the investigation is of a much broader interest because it 
relates to the basic problem of the interaction of a high 
velocity gas stream and a thin liquid film, a situation that 
is frequently encountered in engineering. 

1.2 A Statement of the Problem 

Figure 1 illustrates the process termed liquid-film 
cooling. Presented in the figure is a hot gas stream flow- 
ing along a surface A-B. A liquid coolant is introduced at 
station x = x^ and Is swept downstream over the surface by 
the gas, thereby forming a continuous protective liquid film. 
As a result of the evaporation of liquid In the film, to- 
gether with the possible entrainment of liquid droplets Into 
the gas stream, the liquid film will terminate a«. some point, 
say, x * Xg. The wall region wetted by the liquid will here- 
in be termed the liquid-film cooled region . 

Since the gas and vapor mixture which Is convected down- 
stream is relatively cool, the wall region immediately down- 
stream of the liquid film will be Insulated to some degree 
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from the hot gas stream. It is convenient to define the 
gas-vapor cooled region as that wall region over which the 
temperature of an adiabatic surface would increase from that 
for the liquid film at x = to the recovery temperature 
for the mainstream at, say, x - The extent of the gas- 

vapor cooled region so defined is thus a measure of the 
effectiveness of the gas-vapor cooling process. 

The basic analytical and design problem can be stated 
as that of determining, for given liquid and gas flow param- 
eters, (1) the rate of coolant injection required to estab 
lish a desired wetted area, and (2) the degree of the in- 
sulating effect of the gas-vapor layer downstream of the 
liquid film. The remaining portion of this discussion is 
directed toward a more detailed consideration of each of 
these two problems. 

1.2.1 The Liquid-Film Cooled Region 

The basic problem for the liquid-film cooled region is 
that of being able to predict the wetted surface area as a 
function of the rate of liquid injection and as a function 
of the pertinent flow parameters for both the liquid and gas 
phases. The solution to this problem Is difficult to obtain 
because of the complex phenomena that characterize the flow 
ot a high velocity gas over a liquid film. 

The subject problem is outlined more specifically as 
follows: For the case where a high velocity, hot gas stream 

Is flowing over a thin (of the order of a few thousandths of 
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an inch thick) liquid film, the mass transfer from that film 
is made up of two fundamental contributions: (a) the mass 

transfer that results from the transfer of energy to the 
liquid film and the resultant evaporation of the liquid; and 
(b) the mass transfer that results from a physical breakdown 
of the interfacial structure such that liquid droplets are 
entrained into the gas stream. 

The evaporative contribution to the total mass transfer 
rate (item (a) above) can be divided into yet two more dis- 
tinct contributions: (a-l) that contribution which would 

result if the liquid film surface presented a hydrodynami- 
cally smooth, zero velocity boundary to the gas stream; and 
(a-2) the remaining portion of the evaporative mass transfer 
which results because the surface of the liquid film is 
moving at a finite velocity and does r.ot, in general, pre- 
sent a "smooth" surface to the gas stream. (The surface of 
the liquid film generally exhibits a wavy or pebbled appear- 
ance which should serve to intensify the turbulent transfer 
of heat and mass near the gas-liquid interface.) 

Contribution (a-l) is the same as that prevailing in 

the case of simple mass transfer*. That part of the total 

evaporative contribution Is relatively amenable to analysis. 

Reference (4), for example, presents a recent review of the 

extensive analytical and experimental work performed on that 

problem of simple mass transfer. Although there still 

*The process termed herein "simple-mass-transfer cooling" 

Is equivalent to the process that is frequently termed 
"transpiration cooling" In the published literature. 
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remain some problem areas in the thorough treatment of simple- 
mass- transfer cooling, the effect of mass addition to a bound- 
ary layer upon the rate of heat transfer to the surface has 
been adequately described for the purposes of the present 
investigation. 

No systematic attempt has previously been made to ex- 
perimentally determine the individual contribution of each 
of the three aforementioned mass transfer mechanisms (a-i), 
(a-2), and (b) to the overall mass transfer rate. Although 
the fundamental problem is that of being able to predict the 
total mass transfer rate for a particular flow situation, it 
is evident that in order to develop a satisfactory corre- 
lation for that total mass transfer rate, a basic understand- 
ing has to be developed of the nature of each of the three 
contributions discussed above. 

1.2.2 The Gas-Vapor Cooled Region 

Figure 2 serves to introduce the analytical problem for 
the gas-vapor cooled region. Illustrated in Fig. 2 are the 
pertinent wall temperature profiles immediately downstream of 
the terminus of the liquid film. The upper. solid curve pre- 
sented represents the wall temperature profile for the case 
where the wall is adiabatic*. The lower solid curve repre- 
sents a characteristic p.ofi.e for the case where the wall is 
co.oled externally. As indicated in the figure, it Is con- 
venient to base the heat transfer coefficient for gas-vapor 

♦The term "adiabatic" means that the rate of heat transfer 
through the wall to the surroundings Is negligible. 
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cooling, denoted by h Qy , on the temperature difference 

(T -T ). When that is done, the rate of heat transfer to 
aw w 

an externally cooled wall is given by 

■ h gv < T aw-V f 1 - 1 ’ 

It follows, therefore, that the basic problem is to deter- 
mine, for given flow parameters, (l) the adiabatic wall tem- 
perature profile, and (2) the local heat transfer coefficient 



1. 3 A General Review of the Literature 
Research that has been performed that is pertinent to 
the subject problem can be conveniently grouped into two 
principal areas: 

(a) research concerned primarily with the investigation of 
the detailed nature of the flow of a gas over a thin liquid 
film; and 

(b) research concerned with an actual application of liquid- 
film cooling. 

The first area of research includes the investigation 

of such fundamental phenomena as interfacial structure and 

film Instability* (5, 6, 7, 8, 9), entrainment (10, 11, 12), 

wetting characteristics (13, 14, 15), and pressure Crop in 

♦Film instability is defined as the occurrence of large 
scale disturbances oti the surface of the liquid film. 

These disturbances are often characterized as roll waves. 
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annular, two-phase flow (16, 17). In general, those investi- 
gations were characterized by ambient pressure and tempera- 
ture, essentially zero heat transfer rates, and relatively 
low interfacial shear forces. Such conditions are signifi- 
cantly different from those that characterize the typical 
application of liquid-film cooling. Furthermore, few data 
exist for the case where the liquid phase was not water. 
Although the experimental advantages of using water as the 
liquid phase are obvious, the fact remains that few other 
liquids have physical properties (particularly surface 
tension) that are similar to those of water. Considerable 
doubt thus exists regarding the application of results that 
were obtained with water as the liquid phase to other 
liquids. A detailed review of the literature in this first 
area of research is not presented here because of the ques- 
tion as to the relevance of the results--particularly in a 
quantitative sense--to the present investigation. The 
interested reader is referred to the excellent reviews pre- 
sented In references (5, 10, 18). However, several specific 
references are made to those experimental results throughout 
this report when that reference gives some insight into the 
interpretation of the experimental results of the present 
investigation. 

The research of the second area (namely, applications 
of liquid-film cooling) is characterized by relatively high 
rates of heat and mass transfer between the gas stream and 
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the liquid film (19, 20, 21, 22, 23, 24). The objectives 
of those research programs were basically (a.) to establish 
the feasibility of liquid-film cooling the thrust chamber 
and/or nozzle of a rocket motor; and (b) to determine the 
relationship between the rate of coolant consumption and the 
resultant wetted surface area as a function of the gas stream 
Reynolds number, the gas stream pressure and temperature, 
the properties of the liquid coolant, etc. Attempts to 
develop empirical correlations for the experimental data 
have met with only limited success. The apparent reason for 
the lack of success has been the failure to incorporate into 
those correlations the influence of the interfacial phe- 
nomena* on the rate of mass transfer from the liquid film 
to the gas stream (25, 26). Moreover, theoretical analyses 
developed for liquid-film cooling (e.g., 20, 21) have either 
lacked agreement with the experimental data or have resulted 
in correlations that were not justified physically; the pre- 
vious analyses for liquid-film cooling have been critically 
reviewed in Ref. (26). The failure of the theoretical 
analyses is apparently due to the fact that they have as- 
sumed a direct correspondence between the problem of liquid- 
film cooling and that of simpl e-mass- transfer cooling. Ana- 
lytical models which had resulted in satisfactory ^clyses 
for the case of simple-mass transfer were also assumed to 
describe the heat and mass transfer characteristics for the 

| *That is, the effect of entrainment and of film surface 

{ roughness. 
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case of liquid-film cooling. Figure 3 aids in an explana- 
tion of the motivation for such an approach. Illustrated 
in the figure is a control surface "S" placed at the surface 
of an evaporating liquid film and at the surface of a 
smooth porous wall through which mass is being transferred 
into the gas stream. Assume that for each case the flow 
parameters, the boundary conditions, and the rate of mass 
transfer across the control surface into the gas stream are 
the same. It is then reasonable to argue that the follow- 
ing procedure can be employed to analyze the problem of 
liquid-film cooling. Thus 

(a) The simple theory for mass transfer is employed to 
calculate the rate at which energy is transferred across 
the control surface "S" to the liquid film. The simple 
theory suggests that the rate of heat transfer to the liquid 
film decreases as the rate of mass transfer from the liquid 
film into the gas stream increases. That phenomenon is 
sometimes referred to as thermal "blockage". 

(b) The rate of mass transfer from the liquid film into the 

gas stream is determined by writing an appropriate energy 

balance across the gas-liquid interface*. This step Is 

coupled to step (a) through the "blockage" phenomenon so 

that there is only one value for both the rate of heat 

transfer to the liquid film and the rate of mass transfer 

Into the gas phase that Is predicted for any particular 

*See Subsection 3.1.2 for a discussion on the Interfacial 
energy balance. 
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liquid-film cooling 



smooth porous wall 


simple-mass- transfer 
cooling 


Fig. 3 Comparison of liquid-Film Cooling 
and Simple-Mass-Transfer Cooling 
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situation. The following paragraph discusses how the ex- 
perimental data of the primary references agree with the 
foregoing argument. 

Emmons (20) concluded from his experimental investi- 
gation that the analogy between simple-mass-transfer cooling 
and liquid-film cooling was valid. Emmons suggested that 
the rate of heat transfer to an evaporating liquid film was 
always less than the rate of heat transfer to a dry wall, 
all other conditions being the same. He attributed this, of 
course, to thermal blockage. However, that conclusion drawn 
by Emmons is contradicted by the experimental data obtained 
by other investigators*. Kinney, Abramson, and Sloop (19) 
have reported experimental data showing that the rate of 
heat transfer to a liquid film can be substantially greater 
(by a factor of 3 or more) than the rate of heat transfer to 
a dry wall**. Lulkov (27) obtained similar results from 
well controlled experiments in which extreme care was taken 
to ensure Identical gas flow conditions and surface tempera- 
tures for (a) dry surface and (b) liquid surface test 
sections. The coefficient of heat transfer to the liquid 
film was consistently 20 to 60 percent greater than the co- 
efficient of heat transfer to the dry wall. Hermann, 
Leltlnger and Melnik (24) considered the problem of film- 

* Evidence Is presented In Section 4.4 of this report that 
the data of Emmons do not. In fact, support the analogy 
between simple-mass-transfer and liquid-film cooling. 

**The data of Ref. ( 19 ) are analyzed In Section 4.3. 
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cooling a convergent-divergent nozzle with water*. They 
found that the amount of coolant required to completely wet 
the convergent portion of the nozzle was roughly 2 times 
that predicted by the theory for simple mass transfer. They 
also estimated that the amount of coolant required for the 
divergent portion of the nozzle was 3.3 times that predicted 
by the simple theory**. 

The foregoing discussion raises a serious question 
regarding the adequacy of a simple-mass-transfer model for 
the analysis of liquid-film cooling. The difficulty of ob- 
taining an adequate empirical correlation or developing a 
satisfactory theoretical analysis for liquid-film cooling on 
the basis of the previously reported data for the rate of 
mass transfer from a thin liquid film to a high velocity gas 
stream is due primarily to: 

(a) the 1 imi ted' amount of data that has been obtained for 
conditions similar to those encountered in an application of 
liquid-film cooling; 

(b) the contradictory trends exhibited by the data of the 
primary references; 

(c) the possibility of appreciable experimental error due 
to the severity of the flow parameters Investigated In many 
cases; and 

* The data of Ref. (24) are analyzed in Section 4.5. 

**The discussion of Section 4.5 (In particular, that of 
Subsection' 4.5.5) suggests that the estimates by Hermann, 
et al. (24) of the mass transfer rates for both the sub- 
sonic and supersonic portion of a nozzle cannot be gener- 
ally applied. 
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(d) the scarcity of data for film coolants other than water. 

The present research program was initiated because of 
the incomplete understanding of the phenomena that charac- 
terize liquid-film cooling, and because of the limited in- 
formation that could be obtained from the existing literature. 

1 • 4 The S cope of th e I nves tigati on 

The subject investigation had the primary objective of 
experimentally determining the net rate of mass transfer 
from a thin liquid film to a high velocity gas stream under 
conditions of relatively high pressure, temperature, and 
interfacial shear forces with significant variations in the 
physical properties of the liquid phase. Two cases were 
investigated: (a) the case where the velocity of the gas 

stream was essentially constant over the length of the 
film-cooled region, and (b) the case where the gas stream 
was strongly accelerated over the wetted surface. The latter 
case was investigated because of the interest in the appli- 
cation of liquid-film cooling to cool the exhaust nozzles of 
high-energy propulsion devices. 

A limited number of experiments were conducted for 
Case (a) in which the gas stream was not heated. These 
experiments served to (1) isolate the mass transfer that 
was due to entrair.::.v.:t from the mass transfer that was due 
to the evaporation of the liquid In the film, and to (2) 
significantly extend the range of values for the pertinent 
flow parameters that was investigated in the present 
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Investigation. 

A correlation was developed for the mass transfer data*, 
isolating the contributions to the total mass transfer rate 
that were due to both (a) the simple mass transfer, and (b) 
the interfacial phenomena. The data correlation gives an 
explicit expression for the net rate of mass transfer from 
a thin liquid film to a high velocity gas stream in terms 
of the pertinent variables for both the gas and the liquid 
phases . 

The maximum temperature that the liquid film obtains 
when it 1$ introduced onto the surface of an adiabatic wall 
was measured for each experimental test that was conducted. 

A simple theory was developed for the prediction of the 
maximum film temperature that was in good agreement with the 
experimental data. The theory can be employed to estimate 
the maximum liquid-film temperature for thy general case 
where the film temperature is not a known parameter. 

A number of still photographs were obtained in the in- 
vestigation for the purpose of making some qualitative 
observations about how the pertinent liquid and gas stream 
parameters Influence the character of the Interfacial 
structure. Those observations were utilized In the develop- 
ment of the aforementioned correlation for the experimental 
mass transfer data. 

♦Approximately 600 datum points are reported herein for the 
rate of mass transfer from a thin liquid film to a high 
velocity gas stream. 


« 
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The utility of the design formula (i.e., tho corre- 
lation for mass transfer) that was developed for liquid-film 
cooling was investigated by attempting to correlate the 
published data of the primary references (19, 20, 24). In 
many cases, the values of the pertinent flow parameters 
(e.g., gas stream temperature, pressure, velocity and film- 
cooled length) investigated by those references differed 
substantially from those of the subject investigation. 
Analysis of those data, therefore, gave considerable insight 
into the validity of the design formula proposed herein for 
liquid-film cooling. 

A subordinate objective of this investigation was to 
present an analysis for the gas-vapor cooled region immedi- 
ately downstream of the terminus of the liquid film. The 
analysis presented is essentially the same as that previ- 
ously reported by Gater, L'Ecuyer, and Warner (26). Since 
that original presentation, however, the analysis has been 
further evaluated and additional comments regarding its 
usefulness can be made. 
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2. THE EXPERIMENTAL INVESTIGATION 

2 . 1 The Basic Experimental Apparatus* 

Figure 4 illustrates schematically the experimental 
apparatus that was employed in the investigation. The basic 
components of the apparatus were a hot gas generator and a 
two-dimensional mixing section, ar. adiabatic approach sec- 
tion, a film-cooled test section, and a variable area ex- 
haust nozzle. 

The hot gas generator employed methyl alcohol and air 
as the reactants. Because so much diluent air was added to 
the combustion gases before they entered the experimental 
tunnel, the properties of the gas stream were assumed to be 
the same as those for dry air. 

Prior to entering the test section, the gas stream 

passed through an assortment of fine mesh screens in the 

mixing section, thereby appreciably reducing the turbulence 

level and the non-homogeni ties of the flow.' The gas then 

passed through a 36-inch long adiabatic approach section so 

that the remaining irregularities in the gas stream had a 

chance tu dissipate and the velocity boundary layer had an 

*The experimental apparatus is discussed in greater detail 
in Appendix B. 
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opportunity to develop before the liquid film was introduced. 

(The thickness of the velocity boundary layer at the point T 

f 

of liquid injection was on the order of 0.2 to 0.4 inches.) 

The inside surface of the bottom wall of the approach see- « 

tion was surface ground so that a hydrodynami cal ly smooth 
boundary was presented to the developing boundary layer. 

Figure 5 shows schematically the test section employed 

\ 

in the investigation. The nominal internal dimensions of 
the test section are 2 in. by 5.5. in. cross section by 18 
in. long. Mounted in each side wall of the test section were 
6 circular Pyrex windows of 1 in. O.D. A single oblong 
Pyrex window was mounted in the top wall of the test section. 

These windows permitted visual and photographic observation 
of the liquid film. The photographic technique that was I 

employed is illustrated in Fig. 5. The test plate to which 
the liquid film was confined was recessed into the bottom 
wall of the test section. The liquid coolant was introduced ? 

onto the test plate through an 1/8-in. thick piece of 
"regimesh" that was 1/2 in. long and 2 in. wide. This in- 
jector design provided a uniform distribution of coolant and 
also kept the normal component of velocity for the liquid 
at the injector face to something less than one ft per sec. ; 

The liquid film was confined to the central portion of the 
test plate by means of a trough that was machined into the j 

plate. The trough was 2 in. wide and 0.030 in. deep, and 
it extended from the point of liquid injection to a point 1 
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Fig. 5 Schematic Diagram of Experimental Test Section 
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10 inches downstream. At that point the liquid film was 
mechanically terminated by means of a knife-edge capture 
slot. The projected area for the capture of the liquid film 
was 2.25 in. wide* by 0.060 in. high**. In preliminary 
experimentation with the effectiveness of the capture slot 
technique, it was found that it was necessary to allow some 
of the gas stream to blee^ through the slot with the liquid. 
Attempts to withdraw only the liquid film resulted in spill- 
age of the liquid film over the capture slot. The signifi- 
cance of this bleed gas flow is discussed in Section 2.5. 

In addition, it was visually observed that a small amount of 
liquid frequently spilled around the outside corners of the 
capture slot. This was particularly true if the rate of 
liquid injection was very high. To eliminate that spillage, 
the streamlines of the gas flow were made to converge 
slightly in the immediate neighborhood of the capture slot. 
Then, because the flow behavior of the liquid film was 
foun^. to be strongly affected by the flow pattern of the 
gas stream, the liquid film was made to converge slightly 
toward the center of the capture slot thereby preventing 
spillage of liquid around the corners of the slot. For most 
of the experimental tests in which the gas stream velocity 
was essentially constant over the length of the wetted test 

* 0.25 inch wider than the liquid film 

^approximately an order of magnitude greater than the 
maximum thickness that was investigated for the liquid 
film 
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plate, the convergence of the streamlines for the gas flow 
was accomplished by replacing the last circular Pyrex window 
in each side wall of the test section with an 1/2-in. thick 
by 1-in. O.D. stainless steel plug to which was welded an 
1/8-in. thick by 1-in. wide by 1-1/2-in. long stainless 
steel plate. The degree to which these plates deflected the 
gas stream was regulated by rotating the plates, thereby 
permitting the projected area of each of the plates to be 
varied from 1/8 in. by 1-1/2 in. to 1 in. by 1-1/2 in. 

These plates were normally adjusted so that the liquid film 
would converge from a width of 2 iv' ts at a point approxi- 
mately 1/2 inch upstream of the leading edge of the capture 
slot to a width of 1-1/2 inches at the capture slot. The 
behavior of the liquid film for the remainder of the tests 
that employed a gas stream of constant velocity and for all 
of those tests where the gas flow was accelerated over the 
test plate was controlled by mounting two pieces of 1/32-in. 
thick stainless steel sheet directly on the test plate on 
each side of the capture slot. These sheets were formed so 
that they protruded 1/4 in. into the gas stream and caused 
the gas flow (and, therefore, the liquid film) to converge 
toward the center of the capture slot in the manner as pre- 
viously described. Both of the aforementioned control 
devices were equally effective in the control of the be- 
havior of the liquid film In the Immediate neighborhood of 
the capture slot. 


The gas and liquid mixture that was withdrawn from the 
test plate was passed through a cyclone separator. The gas 
phase was vented off and the liquid phase was collected in 
a beaker. Thus, by measuring the rate of liquid injection 
and the rate of liquid withdrawal, the net rate of mass 
transfer from the liquid film to the qas stream was deter- 
mined. Also determined during the course of each mass 
transfer measurement were the temperatures for the liquid 
at the point of injection and at the point of withdrawal. 
The liquid temperatures were measured with thermocouples 
and they were employed to accurately specify the energy 
balance for the liquid film and to determine the values for 
the pertinent physical properties for the liquid phase. 

For those experiments where the gas stream was accel- 
erated over the length of the wetted test plate, the accel- 
eration of the gas stream was accomplished by inserting an 
aluminum wedge into the top of the test section (see the 
dashed line in Fig. 5). Five static pressure taps were 
installed in the aluminum wedge so that the static pressure 
(and, therefore, the gas stream Mach number, velocity, and 
static termperature) could be determined as a function of 
the coordinate x. In those experiments the nozzle was re- 
moved from the test section and the flow was choked at the 
minimum area formed by the wedge. 

In the subject investigation, the experimental appara- 
tus was designed so that extraneous heat transfer and fluid 
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dynamic effects could be safely neglected, rather than try- 
ing to correct for them analytically. The liquid film was 
injected so as to wet only the center 2 inches of the 5.5-in. 
wide lower surface of the test section, thereby permitting 
the influence of the complex flow phenomena for the corners 
of the rectangular duct on the rate of convective heat 
transfer to the liquid film to be neglected. In establish- 
ing an energy balance for the liquid film, it is difficult 
to determine the amount of energy transferred through the 
wetted test plate to the liquid film in a lateral direction. 
Therefore, in order to substantially reduce the lateral con- 
duction of heat through the test plate to the liquid film, 
two insulating grooves (0.060 in. wide) were cut into the 
underside of the test plate. The grooves extended to within 
0.035 in. of the surface that was wetted by the hot gas, and 
were placed just outside the region that was wetted by the 
liquid film. The insulation grooves, extending from the 
poi it of liquid injection to the point of liquid withdrawal, 
increased the resistance of the plate to lateral heat con- 
duction by roughly a factor of 50. 

2.2 The Experimental Model 
Figure 6 illustrates schematically the experimental 
mcdpl for the subject investigation. The figure presents a 
high velocity, hot gas stream flowing over a horizontally 
mounted flat plate. The film coolant was Injected at a 


Gas Stream 
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rate m 1 (lb/sec-ft)* at x 1 = 40 inches and was withdrawn 
at a rate (lb/'sec-ft) at = 50 inches. The film-cooled 
length, denoted by 1 **, was equal to 10 inches throughout 
the investigation. 

The gas stream flow conditions that were investigated 
were such that the velocity boundary layer developing on the 
flat plate could be treated as being fully turbulent from 
the leading edge of the plate. The wall surface upstream of 
the point of liquid injection was essentially adiabatic so 
that a thermal entrance lengti x t c X 1 = inches was 
established. The thickness of the velocity boundary layer 
at the point of liquid injection was on the order of 0.2 in. 
to 0.4 in. Due to the fact that a velocity boundary layer 
developed on each wall of the experimental tunnel, a correc- 
tion had to be made in the calculation of the gas stream 
velocity to account for displacement effects. The calcu- 
lations showed that the gas stream velocity generally in- 
creased by approximately 5 percent In going from x * 0 to 
x x t - For those investigations in which the gas stream was 
accelerated, the acceleration began at a point 6 Inches up- 
stream of the point of liquid Injection and continued to a 

* The units (lb/sec-ft) follow from dividing the rate of 
liquid Injection (lb/scc) by the width of the Injector 
(ft). 

**In general, the film-cooled length will be denoted by 
"l". The term 1 Q refers to the special case of the 

present Investigation where the film-cooled length was 
fixed at 10 inches. 
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point 2 inches downstream or the end of the film cooled plate. 

As indicated in Section 2.3, the liquid film was ter- 
minated by means of a knife-edge capture slot. Thus, unlike 
previous investigations of liquid-film cooling, the film- 
cooled length for the present investigation was a fixed pa- 
rameter. That characteristic of the experimental model 
allowed the influence of the rate of liquid injection on the 
cooling process to be investigated independently of the other 
flow oarameters. Moreover, it enabled the wetted surface 
area to te accurately determined so that the experimental 
data could be properly evaluated. In previous experiments 
(e.g., Refs. 20-24), the wetted surface area was determined 
by injecting a film of known width and then estimating the 
average film cooled length by means of thermocouples that 
were embedded into the test section wall. However, in such 
experiments it has been found that it is very difficult to 
obtain a uniform termination of the liquid film (19, 20, 28). 
That is, when the film becomes progressively thinner as a 
result of the transfer of mass from the film to the gas 
stream, it reaches a point where surface tension effects 
and/or flow non-uniformities in the gas stream cause the 
lir/.iid film to br^ak down locally. It is, therefore, gener- 
ally very difficult to properly interpret such external wall 
temperature measurements. 

An experimental technique that was introduced into the 
subject investigation was that of preheating the liquid 
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before injection to the maximum temperature that the film 
assumes when it is in equilibrium with the gas stream. 

(When the wetted surface is diabatic, that maximum film 
temperature is essentially the wet-bulb temperature for the 
liquid.) The advanatages of this experimental technique are 
essentially twofold: First, because the physical properties 

for the liquid phase are rather strong functions of the 
liquid temperature, it is much easier to accurately specify 
the values for the physical properties of the liquid if the 
liquid film is isothermal. The second advantage to pre- 
heating the liquid before injection is that the energy 
balance at the qas-liquid interface can then be accurately 
specified. In previous experimental investigations only an 
'’average 11 energy balance for the entire film could be pre- 
scribed because the liquid was introduced at a relatively 
low temperature, thereby making it impossible to ascertain 
v/ hai. portion of the local heat transfer resulted in an in- 
crease in the sensible enthalpy of the liquid and what 
portion resulted in the evaporation of the liquid. To de- 
termine the local rate of heat transfer to the liquid film, 
however, it is necessary that the local energy balance be 
prescribed more accurately. The problem of specifying that 
energy balance is simplified considerably by injecting the 
liquid at the maximum temperature it will assume when it is 
in thermodynamic equilibrium with the hot gas stream. 


Due to limitations in the heating capacity of the heat 
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exchanger that was employed, the liquid could not always be 
preheated completely to the maximum film temperature. But 
the liquid was always injected at a temperature substantially 
greater than the ambient temperature so that the afore- 
mentioned advantages to the preheating of the liquid coolant 
were largely realized. 

2 . 3 The Experimental Parameters 
In order to minimize the experimental error and to ob- 
tain experimental data for mass transfer of sufficient accu- 
racy to elucidate the phenomena and the parameters of impor- 
tance, some sacrifice had to be made in the severity of the 
flow parameters that were investigated, especially the gas 
stream temperature. A number of datum points exist from 
previous investigations for the case of relatively high gas 
stream temperatures (to 4100R), but the accuracy of those 
datum points has necessarily suffered because of the mea- 
surement problems involved. Moreover, such severe gas 
stream temperatures place a physical limitation on the num- 
ber of experimental datum points that can be obtained with 
a reasonable degree of effort. For those reasons, therefore, 
since both accuracy and quantity were badly needed for 
experimental data on liquid-film cooling, the gas stream 
temperature for the subject research program was limited to 
600F. However, this temperature still resulted in an appre- 
ciable rate of evaporation of the liquids investigated 
because of the relatively high gas stream Reynolds numbers 
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that were employed. 

As mentioned previously, experiments were conducted 
for both the case where the gas stream velocity was essen- 
tially constant over the length of the film-cooled test 
section and for the case where the gas stream was accelera- 
ted over the length of the wetted test section. Moreover, 
a limited number of experiments were conducted for the 
first case for which the gas stream was not heated. Those 
experimental tests are termed herein the cold-flow experi- 
ments . 

The liquids employed in the investigation were water, 
methanol, butanol, and RP-1 (a hydrocarbon fuel that is 
essentially kerosene). These liquids resulted in a rather 
wide variation in the pertinent physical properties*. For 
example, the liquid viscosity was varied by roughly a 
factor of 15:1, the surface tension by 3:1, and the heat of 
vaporization by 4:1. The significant variation in the 
physical properties for the liquid phase enabled an accurate 
determination of the influence that each property exerted 
on the heat and mass transfer process. 

Table 1 summarizes the nominal flow parameters for 

those experiments in which the velocity for the gas stream 

was essentially constant over the length of the wetted test 

plate. The table also serves to introduce some pertinent 

♦Appendix D presents the heat of vaporization, the liquid 
viscosity, the surface tension, and the vapor pressure as 
a function of the liquid temperature for the liquids that 
were investigated. 
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Table 1 


Nominal* Flow Parameters--Constant 

Velocity Gas Stream 

T , Gas Stream Static Temperature 


Hot Flow 

400 F, 600 F 

Cold Flow 

40 F 

p, Gas Stream Static Pressure 

75 psia, 150 psia 

Ug, Gas Stream Velocity 


Hot Flow 

50-400 fps 

Cold Flow 

40-200 fps 

Re, Gas Stream Reynolds Number 

10 6 -10 7 per ft 

1 , Film Cooled Lenqth 
o 

10 inches 

Liquid Film Coolants 

Water 


Methanol 


Butanol 


RP-1 

(Hydrocarbon Fuel ) 


♦See Table 6.1 of Appendix G for a Tabulation of the Specific 
Flow Parameters for Each Experimental Test With a Constant 
Velocity Gas Stream. 
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notation that is employed throughout the remainder of this 
report . 

Table 2 summarizes the nominal values of the pertinent 
flow parameters for the experimental tests in which the gas 
flow was accelerated over the wetted test plate. 

2.4 Photographic Data for the Film Surface Characteristics 

A number of photographs of the surface characteristics 
of the liquid film were obtained for a significant range of 
values of the pertinent liquid and gas flow parameters. A 
35 mm camera with a close-copy lens, together with a 3- 
microsecond strobe unit for the light source, was employed 
to photograph the surface of the liquid film. The camera 
was focused on the surface of the liquid film through the 
Pyrex window that was mounted in the top wall of the test 
section (refer to Fig. 5). The strobe unit was focused on 
the surface of the liquid film through one of the side 
windows in the test section. All of the photographs pre- 
sented herein were taken at a point 6 inches downstream of 
the point of liquid injection. 

The surface of a thin liquid film over which flows a 
high velocity gas stream is generally characterized by a 
large number of small scale disturbances (or waves) which 
are essentially three dimensional in character These small 
scale disturbances give the surface of the liquid film a 
"pebbled" or "cross-hatched" appearance. If large scale 
disturbances (instability waves) develop, they travel over 
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Table 2 

Nominal* Flow Parameters--Accel erating Gas Stream 

Tg, Gas Stream Stagnation Temperature 400 F 

p°, Gas Stream Stagnation Pressure 50, 75, 100, 125, 

150 psia 


Ug, Gas Stream Velocity** 
p, Gas Stream Static Pressure** 


200-660 fps 
50 - 40 psia 
75 - 60 psia 
100 - 82 psia 
125 - 105 psia 
150 - 123 psia 


M a , Gas Stream Mach Number** 

a 

1 Q , Film Cooled Length 
Liquid Film Coolants 


0.14-0.51 
10 inches 
Water 
Methanol 


* See Table G.2 of Appendix G for a Tabulation of the 
Specific Flow Parameters for Each Experimental Test With 
an Accelerating Gas Flow. 

**The Numbers Indicate the Variation in the Indicated 
Parameter Over the Length of the Film Cooled Plate. 
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this pebbled substructure at a velocity which is greater than 
the velocity of the smaller scale disturbances (10). A num- 
ber of investigators (e.g., 19, 29-31) have proposed a cri- 
terion for the onset of these large scale disturbance waves 
on the surface of the liquid film. They have proposed that 
if the Reynolds number for the liquid film exceeds a certain 
critical value, then the large scale disturbance waves will 
develop on the surface of the liquid film. Furthermore, it 
has been suggested (e.g., 9, 20, 2l, 22, 29, 30) that appre- 
ciable entrainment of liquid from a liquid film into a high 
velocity gas stream should occur only if the liquid-film 
Reynolds number exceeds that critical value. First it should 
be commented that the photographic data of the present in- 
vestigation did not allow the suggested criterion for film 
instability to be evaluated. It should also be noted, how- 
ever, that the data for mass transfer (discussed later in 
Section 2.5) did not confirm the existence of a critical 
liquid-film Reynolds number with reference to the phenomenon 
of entrainment. In fact, it was observed that the liquid 
viscosity (which enters into the definition of the liquid- 
film Reynolds number) had no measurable influence on the 
rale of mass transfer from a thin liquid film to a high 
velocity gas stream. Indeed, it is significant that the 
experimental data for mass transfer could be readily corre- 
lated without making any direct reference to the phenomenon 
of film instability. 
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The remainder of this section discusses how the 
pertinent flow parameters for the liquid and gas phases 
influence the gross structural characteristics of the sur- 
face of the liquid film. A limited number of photographs of 
the surface of the liquid film are presented to illustrate 
several important points. It should be noted, however, that 
the observations that are made with reference to the influ- 
ence of the pertinent flow parameters on the film surface 
characteristics are the result of .. analysis of nearly 200 
such photographs, and not just the analysis of those photo- 
graphs that are presented herein. These observations will 
be pertinent to the correlation thatis developed in Section 
3.1 for the expe imental data on mass transfer. 

Subsection 2.4.1, which follows, discusses how the sur- 
face characteristics of a liquid film are influenced by the 
parameters for the liquid phase; namely, the rate of liquid 
injection m^, the liquid viscosity , and the surface ten- 
sion o. Subsection 2.4.2 discusses how the surface charac- 
teristics are dependent on the gas stream velocity u , the 
gas stream mass velocity 6, and the momentum parameter M Q , 
where the parameters 6 and M Q are defined by 

6 = p g u g Ob/ft^-sec) 


( 2 - 1 ) 
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and 


M 0 • V# * Ob f /ft 2 ) (2-2) 

respectively. 

2.4.1 The Influence of the Liquid Parameters 
on the Film Surface Characteristics 

The purpose of the present section is to discuss how 
the following parameters influence the characteristics of 
the liquid film surface: (a) the rate of liquid injection 

m^, (b) the liquid surface tension a, and (c) the liquid 
viscosity . 

Rate of Liquid Injection m^. Figures 7 and 8 present 
photographs of a 1-inch-square surface element of the liquid 
film. The pertinent flow parameters are listed in the fig- 
ures for each photograph. The photographs shown in Mg. 7 
were obtained for a case where the gas stream velocity had ♦. 
the relatively low value of 35 ft/sec. Figure 7 (a) illus- 
trates that for the low liquid flow rate of 0.0144 Ib/sec-ft, 
the film surface is characterized by a uniform distribution 
of essentially two dimensional waves. Figure 7 (b) demon- 
strates that as the liquid flow rate is increased to 0.0822 
Ib/sec-ft, with the other parameters remaining constant, the 
interfacial structure becomes less ordered. However, that 
particular change in the character of the liquid film surface 
with increasing rate of liquid injection was observed only 
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for low gas stream velocities, less than, say, 50 ft/sec. 
Figure 8 presents two photographs for the more important case 
where the gas stream velocity is relatively h ■* g h (115 ft/sec 
in this case). Comparing these photographs with those shown 
in Fig. 7, it is seen that the size and the wave length of 
the disturbance waves on the surface of the liquid film has 
decreased significantly with the increase in the gas stream 
velocity. Moreover, comparison of Fig. 8 (a) and Fig. 8 0) 
shows that the scale o* + he disturbance waves is ^ot notice- 
ably altered by increasing the rate of liquid injection m^ 
from 0.0144 to 0.0016 lb/sec-ft. The study of a large num- 
ber of such photographs Indicated that, with the exception 
of those few cases that were investigated where the gas 
stream velocity was very low, the conclusion could be drawn 
that the basic scale of the small disturbance waves on the 
surface of the liquid film was essentially independent of 
the rate of liquid injection. 

Surface Tension o . Figure 9 presents two photographs 
that were obtained under essentially the same flow conditions 
with two different liquids. The liquid shown In Fig. 9 (a) 

Is RP-1 and the liquid shown in Fig. 9 (b) is water. The 
values of surface tension at the corresponding liquid tem- 
peratures in these photographs are approximately 19 dynes/cm 
for the RP-1 and 61 dynes/cm for the water. These photo- 
graphs demonstrate that the scale of the interfacial distur- 
bances is generally more ordered *or the case where RP-1 Is 
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the liquid. The surface of the water film is characterized 
by two types of disturbance waves, one having a larger scale 
and being more two-dimensional in appearance than the other. 
In fact, water was the only liquid that clearly showed two 
types of disturbance waves. The other liquids that were 
photographed, namely, methanol, butanol, and RP-1, did not 
exhibit that characteristic. Since the surface tension for 
water was significantly greater in all cases than the surface 
tension for the other liquids investigated*, it is suggested 
that this phenomenon of the existence of two distinct types 
of disturbance waves for the water films was the result of 
the relatively high values of surface tension for that 
liquid. 

Liquid Viscosity . The photographs of the subject 
investigation did not permit a direct observation of the 
influence of the liquid viscosity on the surface character- 
istics of the liquid film. For any two photographic experi- 
ments for which the liquid viscosity differed substantially, 
there also occurred a substantial variation in one or more 
of the remaining pertinent flow parameters. However, by 
crudely interpolating between the results of several photo- 
graphs, it was apparent that if the liquid viscosity had an 
influence on the surface characteristics for the liquid film, 
it was not pronounced. Relative to the surface tension for 
the liquid and to the pertinent gas stream parameters (which 


*See Fig. D.l in Appendix 0. 
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are discussed in the following subsection), the liquid vis- 
cosity appeared to be of secondary importance in determining 
the structural characteristics of the gas-liquid interface. 

2.4.2 The Influence of the Gas Stream Parameters 
on the Film Surface Characteristics 

The present subsection discusses how the following pa- 
rameters for the gas stream alter the surface characteristics 
of the liquid film: (a) the gas stream velocity u^, (b) the 

gas mass velocity G (defined by Eq. (2-1)), and (c) the mo- 
mentum parameter M q (defined by Eq. (2-2)). Because these 
three parameters are interrelated, it was not possible to 
conduct a series of experiments wherein any two of these pa- 
rameters were held constant while the remaining parameter 
was varied. It was possible, however, to obtain photographic 
data where one of these parameters was held constant while 
the remaining two parameters were varied. It is thus possi- 
ble, in an indirect way, to arrive at some conclusions about 
how the subject gas stream parameters influenced the film 
surface characteristics. 

Constant G, Variable u and M . Figure 10 presents two 

— — — — - 9 , ,, 

photographs of a liquid film where water is the liquid phase. 
The parameters listed in the figure show that the gas mass 
velocity G is essentially the same for both photographs while 
the momentum parameter M and the gas stream velocity u vary 
substantially. Figure 10 shows that the scale of the inter- 
facial disturbances decreases as the velocity u„ (or the 

9 

momentum parameter M Q ) increases. . It follows, therefore, 
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that the value for the gas mass velocity G (or, the gas 
stream Reynolds number) does not give sufficient information 
to characterize the scale of the interfacial waves. 

Constant u , Variable M Q and G. Figure 11 presents two 
— — — — ■ 

photographs of a liquid film for which methanol is the liquid 
phase. The parameters listed in the figure show that the 
gas stream velocity u g is essentially the same for both pho- 
tographs while the parameters M Q and G vary significantly. 
Figure 11 shows that the liquid film for which the gas mass 
velocity G (and the momentum parameter M Q ) has the largest 
value exhibits the smaller interfacial scale. It follows, 
therefore, that the value of the gas stream velocity u g does 
not give sufficient information to characterize the inter- 
facial film structure. 

Constant M Q , Variable u and G. Figure 12 presents two 
photographs for a methanol liquid film for which the momen- 
tum parameter M Q has essentially the same value while the 
gas stream velocity u^ and gas mass velocity G vary sub- 
stantially. The scale of the interfacial wave pattern in 
each photograph is virtually the same. The photographic 
data suggest, therefore, that the momentum parameter M Q is 
the variable for the gas stream that most governs the resul- 
tant film surface characteristics. This observation has not 
been made previously in the published literature, due in 
part to the fact that very few (if any) detailed investiga- 
tions into film surface characteristics have been made for 
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gas stream pressures and temperatures significantly dif- 
ferent from ambient. That fact has precluded the indepen- 
dent investigation of the influence of each of the parameters 
Ug, M Q , and G on the characteristics of the surface of the 
liquid film. (That is, for a constant value of the gas 
stream density, it is not possible to hold any two of the 
parameters u^, G and M Q constant while the remaining param- 
eter is varied.) The foregoing observation, however, is 
related in a fundamental way to the correlation for the mass 
transfer data that is developed in Section 3.1 of this report. 

2.4.3 Summary 

The observations made in Subsection 2.4.1 and in Sub- 
section 2.4.2 regarding the influence of the pertinent liq- 
uid and gas stream variables on the character of the surface 
of a liquid film can be summarized as follows: 

(a) The interfacial structure is not noticeably altered by 
substantial variations in the rate of liquid injection m^, 
except for those cases where the gas stream velocity is very 
low (less than approximately 50 ft/sec). 

(b) The liquid viscosity apparently does not Influence 
the surface characteristics to any significant degree. 

(c) The surface tension o apparently alters the surface 
wlmi acterlstlcs in that two distinct wave form* wore gener 
ally observed when water was the liquid phase (relatively 
high values for a) that were not observed for either butanol, 
methanol, or RP-1 as the liquid phase (relatively low values 
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for 0 ). 

(d) The scale of the small interfacial disturbances is well 
characterized by the value for the momentum parameter M 
(and not G or u^) for the gas stream; the scale of the inter- 
facial disturbances decreases as the value for the momentum 
parameter M Q increases. 

2.5 The Data for Mass Transfer 
for a Constant Velocity Gas Stream 

The following procedure was employed to obtain experi- 
mental data on the rate of mass transfer from the liquid film 
to the gas stream. The hot gas generator was started and 
preselected values for the gas stream temperature, pressure, 
and air flow rate (gas stream Reynolds number or gas mass 
velocity) were established. The liquid film flow rate was 
established and the electric heat exchanger for preheating 
the liquid was turned on. Sufficient time was allowed to 
attain thermal equilibrium (normally 10 minutes) before the 
mass transfer measurements were made. 

Prior to making quantitative mass transfer measurements 
it was necessary to insure that the liquid film was properly 
wetting the test surface and that the capture slot was func- 
tioning so as to remove all the liquid remaining on the test 
surface. First a visual observation was made through the 
Pyrex window in the roof of the test section to determine 
the wetting characteristics of the liquid film. During the 
course of the investigation it was found that if any 
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non-urn fornsi ties existed in the flow field, the liquid film 
would not remain in the trough that was machined into the 
test plate. Frequently one of the mixing screens w;vjlrf fail 
(due to the high gas stream temperature and the large pres- 
sure drop across the screen), resulting in an irregular gas 
flow field and a corresponding poorly wetted test plate 
The only experimental data that are reported herein are 
those for which the liquid film remained in the trough and 
completely wetted the 2-in. by 10-in. test surface. 

Since the rate of mass transfer from the liquid film 
was determined as the difference between the rate of liquid 
injection and the rate of liquid withdrawal, the necessity 
of an accurate determination of these flow rates is apparent. 
The rate of liquid injection was determined by means of a 
turbine flowmeter. The flowmeter was normally calibrated 
before each experimental test by the time-weight technique. 
During the course of an experimental test, the output of the 
turbine flowmeter was continuously recorded on a stripchart 
recorder. A cyclone separator was employed to measure .he 
rate at which liquid was withdrawn from the test plate at 
the capture slot*. As mentioned previously, the effective 
operation of the capture slot required that a certain amount 
of the gas flow be withdrawn with the liquid, thus producing 
a two-phase withdrawal flow. The gas and liquid phases were 
separated by means of the cyclone separator, thereby 


*See Appendix B.4 for the design of the separator. 
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enabling the measurement of the liquid withdrawal rate. The 
cyclone separator was found to provide effective separation 
of the gas and liquid phases. In preliminary experimentation, 
it was found that under conditions of negligible mass trans- 
fer from the liquid film (cold, low velocity gas flow) the 
rates of liquid injection and withdrawal generally agreed to 
within less than 5 percent. Even with the effective per- 
formance of the cyclone separator, it was recognized that 
the direct contact of the gas and liquid phases during the 
flow through the separator could result in some error in the 
measurement of the rate of liquid withdrawal due to satu- 
ration of the gas flow with liquid vapor. Thus, it was 
desired to keep the rate of gas flow through the separator 
as low as possible while maintaining effective operation of 
the capture slot. The procedure employed is described below. 

Once it was ascertained that the wetting character- 
istics of the film were correct, the liquid flow rate was 
set at the maximum value to be investigated. The valve on 
the separator which controlled the rate at which the gas 

phase was vented to the surroundings (and thus controlled 

♦ 

the rate at which gas was withdrawn at the capture slot) was 
opened a small amount. A measurement was made of the rate 
at which the liquid was being captured at the withdrawal 
slot. The vent valve was then opened further and a second 
measurement of the rate of liquid withdrawal was made. This 
procedure was repeated until a further opening of the vent 
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valve did not increase the rate of liquid withdrawal. 

Figure 13 presents a typical plot of the rate of liquid 
withdrawal as a function of the rate at which air was vented 
from the cyclone separator. The figure illustrates that the 
rate of liquid withdrawal increases rapidly at first with an 
increasing vent flow rate and then reaches a maximum, or 
plateau, such that a further increase in the vent flow rate 
does not appreciably affect the rate of liquid withdrawal. 
The increase in m 2 with vent gas flow rate reflects an in- 
crease in effectiveness of the capture slot in withdrawing 
the liquid from the test surface. When the vented qas flow- 
rate is too low, some of the liquid spills over the capture 
slot and then m 2 is not the total film flowrate at the point 
of withdrawal. As the vented gas flowrate is increased, the 
spillage is reduced until with sufficient vent flow the 
capture slot collects and withdraws all of the liquid on the 
test surface. In all of the experimental tests, the vent 
flow rate was always set equal to or greater than the value 
required to realize no change in the rate of liquid with- 
drawal with an increase in the vent flow rate. For the 
magnitude of the vent flow rate normally required, rough 
calculations showed that even if the gas phase that entered 
the separator was completely dry air, and the gas phase 
that was vented from the separator was completely saturated 
with vapor at the corresponding pressure within the separa- 
tor, the amount of mass transfer that would have occurred 




Fig. 13. Rate of Liquid Withdrawal as a Function 
ot Vent Flow Rate for Separator 
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within the separator would be negligible in comparison to 
the total rate of mass transfer that was normally realized 
from the liquid film to the gas stream (not more than 5 per- 
cent). 

After the vent flow rate was set, the quantitative mass 
transfer data were obtained by decreasing the rate of liquid 
injection in incremental stages and making the corresponding 
measurements of the rate of liquid withdrawal. Thus, during 
any one experimental test, the only parameter that was varied 
was the rate of liquid injection. 

Figure 14 illustrates the typical nature of the experi- 
mental data for mass transfer that were obtained in the 
present investigation. The data shown in Fig. 14 are for 
methanol as the liquid phase and they are for a nominal gas 
stream temperature and pressure of 400F and 75 psia, respec- 
tively. Plotted in the figure is the rate of liquid with- 
drawal, mg, as a function of the rate of liquid injection, 
p^, with the gas stream .elocity, u^, as a parameter. The 
solid curve shown at an inclination of 45° represents the 
condition of no mass transfer; at any point on that curve, 
the rate of liquid injection is equal to the rate of liquid 
withdrawal. The vertical distance between the solid curve 
and any one datum point represents the total rate of mass 
transfer from the liquid film for that particular com- 
bination of flow conditions. Figure 14 illustrates the 
important fact that the rate of mass transfer from the 
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liquid film to the gas stream was found to be a function of 
the rate at which liquid was injected onto the test plate ; 
the theory for simple mass transfer, however, does not pre- 
dict any such dependence on the rate of liquid injection. 
Figure 14 shows further that, to within the accuracy of the 
experimental data, the total rate of mass transfer from the 
liquid film to the gas stream increased linearly with the 
rate of liquid injection. This fact was not foreseen before 
the experimental investigation was begun, but it greatly 
simplified the correlation of the experimental data, as is 
discussed in detail in Section 3.1 of this report. 

The remainder of this section discusses how the perti- 
nent flow parameters for both the liquid and gas phases 
influenced the net rate of mass transfer from the liquid 
film to the gas stream. The mass transfer data of some of 
the experimental tests is considered in detail. Although 
it was not feasible to discuss the data of each experimental 
test of the present investigation, it should be noted that 
the illustrative cases that are considered are representative 
of all of the tests of the present investigation. That in- 
eludes not only those tests for which a constant velocity, 
hot gas stream was employed, but also those tests where the 
hot gas stream was accelerated over the wetted test plate 
and those tests for which the gas flow was not heated (cold 
flow, constant velocity gas stream). An inclusive graphical 
presentation of the data for mass transfer is given in 
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Appendix G. 

Subsection 2.5.1, which follows below, discusses the 
influence of the liquid surface tension and the liquid vis- 
cosity on the total rate of mass transfer. Subsection 2.5.2 
discusses the influence of the gas stream flow parameters on 
the measured rate of mass transfer from the liquid film to 
the gas stream. 

2.5.1 The Influence of the Surface Tension and 
the Liquid Viscosity on the Data for Mass Transfer 

Figure 15 serves to illustrate the influence of both 
the surface tension a and the liquid viscosity p upon the 
experimental data for mass transfer. Presented in the 
figure are the experimental data for mass transfer that were 
obtained for the experimental tests 113, 205, and 37. The 
liquids employed in these three tests were methanol, butanol, 
and water, respectively. The pertinent flow parameters for 
each test are listed in the figure. The dashed curves pre- 
sented in Fig. 15 represent the condition of simple mass 
transfer as computed for the pertinent gas flow conditions*. 
If there had been no entrainment of liquid into the gas 
stream and if the liquid film had presented a smooth bound- 
ary to the developing boundary layer, the experimental data 
would lie along the dashed curves. Thus, the vertical 
distance between any one of the dashed curves and the datum 

♦Subsection 3.1.2 discusses how the simple mass transfer 
rates were evaluated. Table G.l of Appendix G lists the 
average value for the simple mass transfer rate for each 
experimental test that is reported. 
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points for tnat curve represents the contribution to the 
total rate of mass transfer that is due to the interfacial 
phenomena (film surface roughness and entrainment). 

The three cases presented in Fig. 15 were chosen for 
purposes of illustration because the pertinent flow param- 
eters for the gas stream, with the exception of the gas 
stream temperature, T^, are essentially the same while the 
viscosity and surface tension for the liquid phase are sig- 
nificantly different*. 

Compare first the data that were obtained for water 
(test 37) and the data that were obtained for methanol (test 
113). Figure 15 shows that the' mass transfer due to the 
interfacial phenomena is much greater for the case where, 
methanol is the liquid phase. Referring to the values for 
the pertinent parameters listed in Fig. 15, it is seen that 
only the surface tension a is significantly different for 
these two cases. These data suggest, therefore, that the 
surface tension influences to a considerable degree the rate 
of mass transfer from a thin liquid film to a high velocity 
gas stream. 

Comparing the data for methanol (test 113) and the data 

for butanol (test 205) in Fig. 15, we note that the measured 

rates of mass transfer are essentially the same for both 

cases. Moreover, the simple mass transfer contribution Is 

*The temperature employed to evaluate the liquid viscosity 
Wl and the surface tension a was the average of the liquid 

temperature at the point of liquid injection and at the 
point of liquid withdrawal. 
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virtually the same for each case so that the contribution 
due to interfacial phenomena is the same. In comparing 
these two cases, it is important to note that all of the 
parameters listed in the figure are nearly the same, except 
for the liquid viscosity Hj. The liquid viscosity for 
butanol for test 205 is roughly 2.5 times that for methanol 
for test 113. These data suggest, therefore, that the 
liquid viscosity has no appreciable influence on the rate of 
mass transfer from a liquid film to a high velocity gas 
stream. The significance of that result was previously 
commented on in Section 2.4. It was pointed out that a num- 
ber of investigators have proposed that a liquid-film 
Reynolds number (which is inversely proportional to the 
liquid viscosity ) should be of primary importance in 
determining the rate at which liquid is entrained from the 
surface of a thin liquid film by a high velocity gas stream. 
The results of the present investigation, however, do not 
support that proposal. 

It should be emphasized that the aforementioned de- 
pendence of the net mass transfer rate on the surface 
tension a was exhibited by all of the experimental data of 
the present investigation. Including those obtained under 
cold-flow conditions, and those realized for the case where 
the gas stream was accelerated over the wetted test plate. 
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2.5.2 The Influence of the Gas Stream 
Parameters on the Data for Mass Transfer 

Figure 16 serves to illustrate how the pertinent param- 
eters for the gas stream (u g , G, and M Q ) influenced the 
experimental data for mass transfer. The figure presents 
the data for mass transfer that were obtained for the experi- 
mental tests 210, 205, and 213. The flow conditions for 
each test are indicated in the figure. The liquid phase 
employed in each test was butanol. The dashed lines in 
Fig. 16 indicate the condition of simple mass transfer. 

Compare first the experimental tests 210 and 205. 

Figure 16 shows that the total .rate of mass transfer for 
test 205 was significantly greater than for test 2i0. The 
contribution due to simple mass transfer, however, was 
essentially the same for both cases. It follows, therefore, 
that the mass transfer due to the interfacial phenomena was 
greater for test 205 than it was for test 210. The param- 
eters listed in Fig. 16 indicate that the gas stream veloc- 
ity, u , and the momentum parameter for the gas stream, M , 
differed substantially for the two cases, while the gas 
mass velocity, G, was essentially the same for each case. 
Consider next the data shown in Fig. 16 for test 210 and 
for test 213. The parameters listed In the figure indicate 
that the momentum parameter M Q was the same for both cases, 
while the remaining parameters show a varying degree of 
difference. The realized rates of mass transfer for the 









63 


two cases are essentially the same. The data of Fig. 16 
suggest, therefore, that the momentum parameter M Q for the 
gas stream characterizes, to a large degree, the extent of 
the contribution to the net rate of mass transfer that is 
due to the interfacial phenomena. Again it should be em- 
phasized that this observation was supported by all of the 
experimental data of the present investigation, and not by 
just those that are presented in Fig. 16. 

2.5.3 Summary 

The foregoing discussion suggests that the contribution 
to the net rate of mass transfer that is due to the inter- 
facial phenomena: 

(a) increases in a linear manner with the rate of liquid 
injection, m^; 

(b) increases with decreasing values for the surface 
tension a ; 

(c) is apparently not a function of the viscosity for the 
liquid, y-j ; and 

(d) is characterized by the momentum parameter for the gas 
stream, M , with the rate of mass transfer increasing with 
increasing values for M 0 - 

These observations form the basis for the mass transfer 
correlation that is developed in Section o.i or this report. 


64 


2 • 6 T he D at a for Hass Tran sf er for an 
Accelera ti ng Gas Stream 

The stagnation temperature for the gas stream was main- 
tained constant at 400F for that phase of the investigation 
where the gas stream was accelerated over the wetted test 
surface. Moreover, the Mach number profile was fixed, such 
that the Mach number increased from 0.14 at the point of 
liquid injection to 0.51 at the point of liquid withdrawal*. 
The stagnation pressure was varied from 50 to 150 psia in 
increments of 25 psia, and both water and methanol were 
employed as the liquid phase. The experimental procedure 
employed in these investigations was basically the same as 
that outlined in the previous section. 

Figure 17 presents all of the experimental data that 
was obtained in this portion of the investigation. Plotted 
in the figure is the rate of liquid withdrawal as a function 
of the rate of liquid injection with the stagnation pressure, 
p°, as a parameter. The experimental data obtained for both 
water and methanol as the liquid phase are shown in Fig. 17. 
Again, like the data that was presented in Figs. 14, 15, and 
16. the total rate of mass transfer (the vertical distance 
between the line for no mass transfer and any one datum 
point) increases linearly with the rate of liquid injection. 
The data in Fig. 17 demonstrate that for the same flow con- 
ditions (same total pressure, total temperature, and degree 

*i.e., only one area-contraction ratio (one aluminum wedge 
in Fig. 5) was emoloyed in the investigation 



Fig. 17. Experimental Data for Mass Transfer. Accelerating 
Gas Stream. Methanol and Water Film Coolant 
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of acceleration of the gas stream) the rate of mass transfer 
v/as much greater for the case where methanol was the liquid 
phase than for the case where water was the liquid phase. 

That is due, in part, to the fact that the heat of vapori- 
zation for water is roughly 3 times that for methanol, so 
that for a given rate of heat transfer to the surface of the 
liquid film, the rate of vaporization of methanol is greater 
than is the rate of vaporization of water. But, more impor- 
tantly, the rate of entrainment of methanol into the gas 
stream was much greater than was the rate of entrainment of 
water because the surface tension for methanol is roughly 
one-third that of water. The rather strong dependence of the 
rate of mass transfer from the liquid film into the gas 
stream on the surface tension was discussed in Subsection 
2.5.1 of the present report. 

A correlation is developed in Section 3.3 for the mass 
transfer data presented in Fig. 17. It will be shown that 
the data correlation is a natural extension of the corre- 
lation that is developed in Section 3.1 for the experimental 
data for mass transfer for the case where the gas stream 
velocity was essentially constant over the length of the 
wetted test plate. 


i 
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2 . 7 The Data for the Max i mum_ Liquid 
Film Temperatu re 

In general, T^ , the temperature of the liquid at any 
point within the liquid film, is less than T £ , the tempera- 
ture at the surface of the liquid film, and is greater than 

T , the temperature at the surface of the wetted wall. How- 
w 

ever, for the special case where the wetted wall is adiabatic 
(as was the case for the present investigation), it is 
possible to obtain the condition 


That condition occurs when all of the energy transfer to the 

liquid film goes into the vaporization of the liquid in the 

film*. Moreover, at the point downstream from the point of 

liquid injection where that condition is obtained, the 

temperature Tj (=T = T g ) is a maximum. That maximum liquid 

temperature is denoted herein by T-j 

The maximum liquid temperature, T^ m> was determined 

experimentally for the present investigation by measuring 

Tj 2 * the temperature for the liquid coolant at the point of 

liquid withdrawal. Recall that the attempt was made in the 

subject investigation to preheat the liquid before injection 

so that the sensible enthalpy rise that would occur for the 

liquid after it was introduced onto the test plate would 

*That is, none of the energy transfer to the surface of the 
liquid film goes to increase the sensible enthalpy for the 
liquid in the film. 
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be minimized. Therefore, for those tests where it was possi- 
ble to completely preheat the liquid before injection, it 
follows that the temperatures m and T-j ^were equivalent. 
Moreover, for those experimental tests wherein the liquid 
could not be completely preheated before injection (due to 
limitations in the heating capacity of the heat exchanger 
that was employed), it was found that the measured with- 
drawal temperature was not influenced significantly by the 
rate at which liquid was injected onto the test plate, there- 
by suggesting that the rate of energy transfer to the liquid 
film was sufficiently rapid to cause the allowable increase 
in the sensible enthalpy for the liquid to occur before the 
liquid film reached the capture slot. It follows, therefore, 
that the temperaturess 2 an{ * Tj m could be reasonably 
equated for all of the experimental tests of the subject 
i nvestigati on . 

The liquid temperature T^ 2 was measured with a chromel- 
alumel thermocouple*, and the indicated temperature was con- 
tinuously recorded on a stripchart recorder. Experimental 
values for T<| 2 were obtained for each of the four liquids 
that was employed, and for the entire range of values for 
the gas stream parameters that was investigated**. Table 
G.l of Appendix G lists the average value of Tj 2 for each 
of the experimental tests that are reported herein. Also, 

* See Fig. B.6 of Appendix B for the physical location of 
the thermocouple. 

**Refer to TabU 1 in Section 2.3. 
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there is a graphical way in which the data can be presented 
that is of interest. In most previous investigations of 
liquid-film cooling (e.g., 20, 22, 23), the maximum tempera- 
ture for the liquid film was taken to be the boiling temper- 
ature for the liquid at the prevailing pressure. However, 
it is demonstrated theoretically in Refs. (21) and (26), and 
in Section 3.4 of this report, that the temperature of the 
liquid film can never be equal to that boiling temperature. 
Figure 18 graphically illustrates this point for the data 
obtained in the subject investigation. Presented in Fig. 18 
is a plot of the boiling temperature for the liquid at the 
prevailing pressure versus the measured liquid-film 

temperature. Data are presented for water, methanol, butanol, 
and RP-1. Figure 18 shows that the maximum film temperature 
can be substantially less than the boiling temperature at 
the prevailing pressure. The maximum difference between 
these two temperatures encountered in the present investi- 
gation is approximately 355F, and it was realized for RP-1 
at a pressure of 150 psia and a gas stream temperature of 
400F. 


v • 



measured (F) 
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3. CORRELATION OF THE EXPERIMENTAL DATA 

The present section is concerned with the correlation 
of the experimental data that were discussed previously in 
Section 2. It should be stated at the outset that it was 
the intention of the present investigation to develop for- 
mulae that could be employed to accurately design a liquid- 
film cooled system but which would not involve excessive 
computations.. Therefore, the present investigation sought 
to consider the details of the complex phenomena that char- 
acterize the interaction of a high velocity gas stream and 
a thin liquid film only to the extent necessary to develop 
those design formulae. These interfacial phenomena have 
been the object of a large number of experimental and/or 
analytical investigations reported in the literature. To 
illustrate, references (5, 32) were concerned with the ampli- 
tude and wave length spectra that the gas-liquid interface 
exhibits for various liquid and gas stream flow conditions; 
references (33, 34) were concerned, in part, with applying 
stability theory to analyze the phenomenon of film Insta- 
bility; and references (28, 35) investigated the structure 
of the boundary layer which develops over the surface of a 
liquid film. It is virtually impossible; however, to 
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utilize these very interesting data--particular1y in a 
quantitative sense--to answer the basic engineering ques- 
tion: .what is the net rate of mass transfer from a thin 

liquid film to a high velocity, hot gas stream? The rela- 
tively simple formulae developed herein answer that question 
expl i ci tly . 

Section 3.1 presents a correlation of the data for mass 
transfer with a constant velocity gas stream. Correlation 
of the data for mass transfer where the gas stream was 
accelerating over the wetted test surface is discussed in 
Section 3.3 Analysis of the latter data requires first the 
consideration of the problem of how the net rate of mass 
transfer varies with the liquid-film cooleu lenght 1; Section 
3.2 considers that problem. Section 3.4 presents an analysis 
of the data for the maximum liquid temperature. 

3 . 1 Correlation of the Data for Mass Transfer 
for a Constant Velocity Gas Stream 

3.1.1. General Discussion 

The experimental data for mass transfer for the subject 
investigation showed that for fixed gas stream flow con- 
ditions, the rate of liquid withdrawal, mg. could be approxi- 
mately expressed as a linear function of the rate of liquid 
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i njecti on , m^*. Thus 

m 2 = - A + B m 1 (3-1) 

If m^ is subtracted from both sides of Eq. (3-1), the resul- 
tant expression can be written in the form 

m l~ m 2 = ^ + (1-B) m^ (3-2) 

Noting that the quantity (m^ - m 2 ) is simply m', the net 
rate of mass transfer from the liquid film to the gas stream, 
we obtain from Eq. (3-2) the expression 

m' * A + (1-B) m 1 (3-3) 

Equation (3-3) expresses the net rate of mass transfer from 
the liquid film as a linear function of the rate of liquid 
injection. 

In the original statement of the problem in Section 1.2, 
it was indicated that tn' was made up of two fundamental con- 
tributions: (a) the contribution due to the evaporation of 

the liquid in the liquid film, and (b) the contribution due 
to the entrainment of bulk liquid into the gas stream. The 
problem, therefore. Is to determine how each of these two 
contributions can be evaluated from Eq. (3-3). With 
♦Refer to Figs. 14, 15, and 16 and to Appendix 6. 
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reference to that problem, consider the following two points: 

(a) The phenomen of liquid entrainment generally has been 
associated with the shearing off of liquid droplets from 
large scale disturbance waves (instability waves) on the sur- 
face of the liquid film. Therefore, it is reasonable to 
assume that the rate at which liquid is entrained fro the 
surface of these large scale disturbance waves into the gas 
stream would increase with the frequency of occurrence of 
these waves. It has been demonstrated (10) that the fre- 
quency of occurrence of these waves increases with m^, the 
rate at which liquid is introduced onto the wetted wall 
surface. Therefore, the rate of entrainment of liquid into 
the gas stream should be proportional to m^. 

(b) It can be argued that the evaporation contribution to 
the total rate of mass transfer from the liquid film to the 
gas stream should be largely independent of the flow param- 
eters for the liquid phase. The extent to which the tur- 
bulent transport of heat to the surface of the liquid film 
is intensified because of the effective roughness of the 
surface of the liquid film should be primarily a function of 
the physical scale of the substructure on the surface of the 
liquid film. The photographic data discussed in Section 2.4 
suqgested, however, that while the scale of that substructure 
Is influenced to a considerable degree by the flow parameters 
for the gas stream (e.g., u g , 6, and M Q ), it Is not notice- 
ably altered either by sirnificant variations in the rate 
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at which liquid is injected onto the test plate or by 
significant variations in the physical properties for the 
liquid phase. It thus seems reasonable to assume, at least 
as a f^rst approximation, that the evaporative contribution 
is not dependent on the parameters for the liquid phase. 

The foregoing two points and the form of the Eq. (3-3) 
suggest the simple hypothesis that (a) the first term on 
the right-hand side of Eq. (3-3) represents the contribution 
to the total mass transfer rate that is due to the evapo- 
ration of the liquid in the film, and (b) the second term 
on the right-hand side of Eq. (3-3) represents the contri- 
bution that is due to the entrainment of bulk liquid from 
the surface of the liquid film into the gas stream. Thus, 
it is hypothesized herein that 

A = the evaporation contribution 

(1-B) m^ = the entrainment contribution 
The ultimate justification for this simple hypothesis i c 
the fact that it resulted in a satisfactory correlation of 
the experimental data. This will be demonstrated in Sub- 
section 3.1.5 of the present report. 

Referring again to the statement of the subject problem 
presented in Section 1.2, it was indicated that the evapora- 
tive contribution to the total rate of mass transfer m' 
could be subdivided into (a) that due to the simple mass 
transfer and (b) that due to the interfacial phenomena (the 
effective surface roughness). If the term A in Eq. (3-3) 
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gives the total contribution to m' that is due to evapo- 
ration, and if rn^ denotes the rate of simple mass transfer, 
then a parameter r, termed the roughness p a rameter , can be 
defined such that 

A = (l + r) (lb/sec-ft) (3-4) 

The parameter r is a measure of the degree to which the 
interfacial structure intensifies the turbulent transport of 
heat from the gas stream to the surface of the liquid film. 
For the limiting case where the surface of the liquid film 
becomes hydrodynami cal ly smooth, the parameter r goes to 
zero, and the problem reduces to the classical problem of 
simple mass transfer. 

In addition to the roughness parameter r, it is con- 
venient to define an entrainment parameter e such that 

e 0 = 1-B (3-5) 

where B is defined by Eq. (3-1). The discussion of Section 
2.5 suggests that e Q is primarily a function of the surface 
tension a for the liquid phase and of the momentum parameter 
M 0 for the gas stream. The subscript "o" on indicates 
that e 0 was determined from experimental data for which the 
liquid-film cooled length was maintained constant at 10 
inches (1=1= 10 inches). Section 3.2 of the subject 
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report shows how the experimental results obtained for that 
fixed film cooled length can be utilised to analyze the 
general case of liquid-film cooling where the liquid-film 
cooled length 1 is arbitrary. 

Substitution of Eqs. (3-4) and 3-5) into Eq. (3-3) 
results in the following expression. Thus 

m' * (1 + r) m^ + e Q (3-6) 

Equation (3-6) is the primary equation that was utilized in 
the present investigation to correlate the net rate of mass 
transfer from the liquid film to the gas stream. The sig- 
nificance of each of the terms in Eq. (3-6) is shown graphi- 
cally in Fig. 19. The lower solid curve in the figure 
represents a least-squares fit of the data for mass transfer 
that were realized for a typical experimental test. The 
upper solid curve in Fig. 19 represents the condition of no 
mass transfer, and the upper dashed curve is the curve that 
is predicted by the simple theory for mass transfer. The 
lower dashed curve represents the curve that would be pre- 
dicted by the present theory, accounting for the roughness 
of the film, if there was no entrainment of liquid into the 
gas stream. Figure 19 shows how each of the terms in Eq. 
(3-6), m', rm ' , and e m n , can be conveniently interpreted 
graphically. 

Figure 19 and Eq . (3-6) show that the subject problem 
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of correlating the experimental data for m' is reduced to 
that of calculating the simple mass transfer rate m^ and to 
the correlation of the parameters r and e Q on the basis of 
the experimental data. Subsection 3.1.2 discusses the 
method that was employed to calculate m' in the present in- 
vestigation. Subsection 3.1.4 and 3.1.5 discuss the tech- 
niques for correlating the parameters e Q and r, respectively, 
in terms of pertinent flow parameters. Subsection 3.1.6 
then illustrates the agreement of the resultant data corre- 
lation with the experimental mass transfer data of the 
present investigation. 


3.1.2 The Evaluation of the Rate of 
Simple Mass Transfer m^* 

The net rate of simple mass transfer from the liquid 
film to the gas stream is related to the local rate of mass 
transfer by the following expression. Thus 



dx 


(3-7) 


* <*2 - * ’o 


(3-8) 


*The primary reference for this section is the textbook by 
Kays (36) on convective heat and mass transfer. 
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where 

nig = the net rate of simple mass transfer 
= the local rate of simple mass transfer 

1 = the film cooled length ( = 10 inches) 

in" - the average value for m" over the film cooled 
5 length 1 0 s 

A method is outlined in the present section for the 
evaluation of the local rate of simple mass transfer m^. 

The method consists of two basic steps: (a) writing an 

energy balance across the gas-liquid interface so that m^ 
can be directly related to q^, the local rate of heat trans- 
fer from the hot gas stream to the surface of the liquid 
film, and (b) determining q $ from the comprehensive theory 
that is presented by Kays (36) for that problem. The 
solution for m^ is then completed by developing a closed- 
form expression for the average rate of simple mass trans- 
fer, in", that appears in Eq. (3-8). 

The Energy Balance at the Gas-Liquid Interface . Figure 
20 illustrates a general interfacial element placed at the 
surface of a stable, non-reactive liquid film with the dotted 
lines representing the surfaces considered in writing an 
energy balance for the element. The different energy fluxes 
involved are defined as follows: 

(i; q s = (k 3T/3y) s = the rate at which energy is trans- 

3 ferred across the S-surface by 
conduction. 

(2) q, = (k 3T/ ay ) , = the rate at which energy is trans-. 

L ferred across the L-surface by 
conduction. 
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Fig. 20. Control Volume for Interfacial Energy Balance 


w - w 1 



82 


(3) m" h s = the rate at which energy is convected across 

s v * the S-surface by evaporated vapor. 

(4) m" h, = the rate at which energy is convected across 

the L-surface by the evaporating liquid. 

(5) q <. = the rate at which energy is transferred across 

r * the S-surface by radiation. 

(6) q - = the rate at which energy is transferred across 

r ’ the L-surface by radiation. 

The energy balance for the interface element illustrated 
in Fig. 20 can be written as follows: 


q s = q L + ■; 0> v ,s - l> L ) + (q r ,L - Vs> (3 ' 9) 

I II III IV 


Term III in Eq. (3-9) can be rewritten in the form 
■; H v 

where is the latent heat of vaporization for the liquid 
coolant which corresponds to the saturation condition that 
is prescribed by the film surface temperature T . Further- 
more, for the temperature levels encountered in the subject 
investigation, term IV in Eq. (3-9) can be considered neg- 
ligible in comparison to a term such as term III. Thus, 
neglecting that term, Eq. (3-9) can be rewritten to yield 

q$ = q L + H y (3-10) 

The term q L in Eq. (3-10) is that portion of the heat trans- 
fer to the surface of the liqui-d film which crosses the 
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interface and goes into raising the sensible enthalpy of the 
liquid and/or is transferred through the wetted plate to the 
surroundings. In the subject investigation, the test plate 
was essentially adiabatic so that the heat transfer across 
the test plate to the surroundings was negligible. Moreover, 
an attempt was made in the investigation to preheat the 
liquid before injection so that no rise in sensible enthalpy 
would occur for the liquid once it was introduced onto 
the test plate. The elimination of the rise in sensible 
enthalpy for the liquid would make the quantity q^ in Eq. 
(3-10) equal to zero and thus simplify the energy balance 
as mentioned in Section 2.2. However, due to limitations in 
the heating capacity of the heat exchanger employed to pre- 
heat the liquid, some rise in the sensible enthalpy for the 
liquid usually occurred after the liquid was injected. To 
correct for that fact, the term q^ in Eq. (3-10) was approxi- 
mated by the expression* 

q L " m l * h l,2 * h l,l) (3-H) 


where 

h, « * sensible enthalpy for the liquid at the point 
’ of withdrawal , x * x 2 

hi i * sensioie enthalpy for the llquio at the point 
* of injection, x * 


*This approximation is the same as that employed in 
Refs. (19, 29). 
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Substitution of Eq. (3-11) into Eq. (3-10) results in 
the following basic expression for the interfacial energy 
balance. Thus 



H v + h 1 , 2 



(3-12) 


It is convenient to introduce a parameter 4> that is termed 
herein the energy balance parameter , where <{> is defined as 
the group of terms on the right-hand side of Eq. (3-12). 
Thus 


* * H v ♦ h l,2 



(3-13) 



(3-14) 


Since the energy balance parameter <p is expressed solely 
in terms of quantities that can readily be evaluated, the 
problem of determining m" thus reduces to that of evaluating 

q S* 

The Evaluation of q $ . It Is convenient to express q<. 
in terms of St , the local Stanton number for simple heat 
transfer, where St $ is defined by 


St s “ 5”c — (T-r ) 
pg g s' 


(3-15) 
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where 

G - the mass velocity for the gas stream 

c = the specific heat at constant pressure for the 

gas stream 

Tg = the temperature for the gas stream 

T s = the temperature at the surface of the liquid film 

The complete evaluation of the local Stanton number St $ 
requires that the effect of the following on the transport 
phenomena be accounted for: (a) the thermal entrance length 

x t s Xj, (see Fig. 6); (b) the dependence of the physical 
properties for the gas phase on the temperature T and on the 
local concentration of the injected vapor C; (c) the rate of 
transfer of vapor from the surface of the liquid film into 
the gas phase; and (d) non-unity values for the Prandtl 
number Pr for the gas phase. If St $0 denotes the local 
value for the Stanton number for the case where none of 
these effects are present, then the Stanton number St^ can 
be related to St $0 through a number of correction factors 
that are denoted by F. Thus 


St 


s 




Pr 


St 


so 


(3-16) 


wherre 


F * the correction factor for the***"*! entrance 
t length effects (x t * x A ) 

F, . * the correction factor for the dependence of the 

physical properties of the gas on the temperature 
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F c = the correction factor for the dependence of the 
physical properties for the gas on the concen- 
tration of the injected vapor C 

F „ = the correction factor for the injection of the 
mass (i.e., vapor) into the developing boundary 
layer* 

Fp - the correction factor for a non-unity Prandtl 
number for the gas stream 

Kays (36) suggests that for the case where the gas stream 
velocity is constant**, the correction factors in Eq. (3-16) 
can be evaluated by the following empirical equations. Thus 


F x ■ 

x t 

< l -3T> 

(3-17) 

F T * 

j 0.25 

(A 

s 

(3-18) 

F c s 

H °- 4 

(3-19) 

F m" * 

ln(l+B h ) 

(3-20) 

F Pr = 

"U 

1 

o 

• 

(3-21) 


where 

Mg * the molecular weight for the gas stream 

M = the molecular weight for the. gas-vapor mixture 

at the surface of the liquid film 

* F m „ accounts for the thermal “blockage" effect that was 
discussed In Section 1.3. 

**This restriction actually applies to only F . 

x t 
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= a dimensionless heat-mass transfer parameter 


(T„ - TJ 


(3-22)* 


To complete the solution to q<., both the term St $() in 
Eq. (3-16) and the term M g in Eq. (3-19) must be evaluated. 
The following classical expression for turbulent heat trans- 
fer to a flat plate was employed to evaluate St SQ (36). Thus 


St = 0.0295 Re " 
so X 


(3-23) 


where 


= the Reynolds number for the gas stream based 
on the characteristic length x 


(3-24) 


Vi = the dynamic viscosity for the gas stream evalu- 
9 ated at T g 

The molecular weight M is related to M , the molecular 

j V 

weight for the vapor, and M g , the molecular weight for the 
free-stream gas, by 


: S - C s> 


(3-25) 


*If compressibility effects are important (Mach number 
greater than, say, 0.5), then the temperature difference 
i T g - T $ ) should be replaced by (T g>aw - T $ ), where T g>aw 

is the adiabatic-wall recovery temperature for the gas 
stream. 


where 


C = the concentration fo v ' the vapor at the surface 
of the liquid film 

The vapor concentration C $ can be determined by the follow- 
ing approximate expression (36). Thus 


C 


s 



1 ♦ B h 


(3-26)* 


The Evaluation of m^ . Equations (3-14), (3-15), (3-16) 
and (3-22) can be combined to obtain the following basic 
expression for m£ as a function of the coordinate x. Thus 


■! = G B h F x t f T F C F m" F Pr st ,o 


(3-27) 


For the subject case where the gas stream velocity is con- 
stant over the film cooled length 1 , it can be reasonably 
assumed that only the terms F and St in Eq. (3-27) are 

A ^ 5 0 

a function of the coordinate x. Therefore, for that special 
case, the result of the substitution of Eq. (3-27) into 
Eq. (3-7) can be written as 

x 2 

"i ■ G B h F T F C F m" F Pr f F x, St so dx < 3 - 28 > 

X 1 

iu is not possible to evaluate the integral in tq. (3-28) !n 

*See the footnote following Eq. (3-87a) of the present 
report for an explanation of the origin of this equation. 


w 
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closed form if Eq. (3-17) is substituted for F . However, 

x t 

Reynolds, Kays, and Kline (37) have developed the following 
approximate expression for that integral term: 


8 



where St „ « is the value for St_ that is obtained from 
Eq. (3-23) with x = x 2 . For a value of Xj, = 40 inches and 
a value of x 2 = 50 inches, the non-dimensional term inside 
the brackets { } in Eq. (3-29) has a value of 1.374. Sub- 
stituting that value into Eq. (3-29) and the resultant 
expression into Eq. (3-28), one obtains the basic expression 
that was employed in the present investigation to evaluate 
the rate of simple mass transfer from the liquid film into 
the constant velocity gas stream . Thus 

• 1.374 G B h F t F c F n „ F pr St |#>2 (3-30) 

Table G.l of Appendix G lists the average, value for the 
simple rate of mass transfer, m^, that was calculated for 
each experimental test. The range of value* chat was com- 
puted for each of the dimensionless terms in Eq. (3-30) and 
also for the interfacial concentration C & and the heat-mass . 
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transfer parameter is summarized as follows: 

F t = 1.015 to 1.150 

F c = 0.976 to 1.096 

F m „ = 0.86 to 0.96 

Fp r = 1.168 (for all cases) 

St_ 0 = 0.0021 to 0.0030 
r , = 0.081 to 0.247 

B h = 0.083 to 0.328 

3.1.3 The Correlation of the Entrainment 
Parameter e Q 

The present section is concerned with the correlation 
of the entrainment parameter e Q that occurs in the primary 
equation (3-6). The procedure that was employed herein to 
develop the correlation was purely inductive. That is, the 
procedure that was employed was to first note how each of 
the pertinent flow parameters influenced the rate of en- 
trainment of liquid from the liquid film into the gas 
stream, and to then develop a general correlation that would 
reflect each of those observations. Thus an analytical 
model is not postulated in this report for the approximate 
analysis of the entrainment phenomenon. Moreover, the 
method of dimensional analysis was not directly employed to 
develop the subject correlation for e 0 . The ultimate justi- 
fication for the simple procedure that was employed is the 
fact that it resulted in a satisfactory correlation of the 
parameter e Q in terms of the pertinent variables for both 
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the gas and liquid phases. The next subsection of this 
report, Subsection 3.1.4, does, however, consider in some 
detail the possible physical significance of the resultant 
correlation that was developed for e Q . The discussion in 
Subsection 3.1.4 suggests a simple procedure for correlating 
the roughness parameter "r" that was defined by Eqs. (3-4) 
and (3-6). Subsection 3.1:5 presents the correlation for 
the parameter r. 

The discussion of the experimental data for mass trans- 
fer that was presented in Section 2.5, and summarized in 
Subsection 2.5.3, showed that the two parameters which most 
influenced the net rate of entrainment of liquid from the 
liquid film into the gas stream were the momentum parameter 
M q for the gas stream and the surface tension a for the 
liquid phase. The rate of entrainment of liquid increased 
with increasing values for M Q and with decreasing values for 
o. In view of those observations, the attempt was made to 
correlate the parameter e Q in terms of a dimensional entrain- 
ment group "X e " which is defined as* 

. a 

M a T 7 

X- * (t 1 ) ' (3-31) 

S 

where "a" is a constant that was determined experimentally. 

The term (T g /T s ) a/ ^ in Eq. (3-31) arises from the substitution 
for the gas stream density p in the parameter M Q of the 

51 

— 111 ■ ” - ■ " ’ » ■ ■■■■ — 1 ■ 11 ' —" ■■■■ ■■ - ■ - 

♦Evidence is presented in Subsection 3.1.4 that X fi is equiva- 
lent to a dimensionless Weber number. 
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corrected density p (T /T )* 5 . This is a technique that is 

9 9 s 

frequently employed in the correlation of heat transfer data 

to account for the fact that the density for the gas within 

the boundary layer region is inversely proportional to the 

temperature for the gas (36)*. It should be noted, however, 

that because of the low gas stream temperatures that were 

employed in the subject investigation, the term (T /T ) a ^ 

9 * 

was not crucial to the correlation of the experimental data 
for the entrainment parameter e Q (discussed immediately 
below). It is included in the definition of the entrainment 
group X g primarily to be consistent with the accepted prac- 
tice for accounting for variable property effects. 

A lengthy trial -and-error procedure employing all of the 
experimental data presented herein showed that the form of 
the entrainment group presented in Eq. (3-31) resulted in 
the best correlation of the data over the complete range of 
parameters investigated, and that the parameter e Q could be 
•“orrelated satisfactorily if the constant "a" in Eq. (3-31) 
was set equal to H. Figure 21 presents the resultant corre- 
lation of the data for a - h. Presented in the figure is 
(l-e Q ) as a function of X fi , where (l-e o ) was determined for 
each experimental test from the slope of the least-squares 
fit curve for the experimental data of that test (e.g., as 
shown in Fig. 14). The nominal values for the gas stream 
temperature and pressure are indicated in Fig. 21, together 

♦This technique was employed, in part, to derive the ex- 
pression for the aforementioned correction factor F T (see 
Eq. (3-18)). 1 




Entrainment Parameter e„ as a Function of the Entrainment Group 
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with the four liquids that are represented by the data. Al- 
though there is admittedly some scatter in the experimental 
results, some of the scatter can be attributed to the fact 
that ( 1- e 0 ) was determined from a least-squares fit of the 
"raw" experimental data for mass transfer. Indeed, because 
the data presented in Fig. 21 do represent such a substan- 
tial variation in both the pertinent parameters for the liq- 
uid phase and for the gas phase, it is somewhat surprising 
that e Q can be correlated so well in terms of the relatively 
simple parameter X . Perhaps most significant is the fact 

V 

that e Q was not found to show any dependence on the liquid 

viscosity y-j. The data presented in Fig. 21 represent a 

variation in y-j from 1.7x10“^ lb/ft-sec for water at a tem- 

-4 

perature of 265F to 30x10 lb/ft-sec for butanol at a tem- 
perature of 40F. The results of the subject investigation 
are thus contrary to the suggestion made by some investi- 
gators that a liquid-film Reynolds number (which is inversely 
proportional to y^) is the parameter that most characterizes 
the entrainment phenomenon*. 

The curve shown in Fig. 21 is given by the following 
expression. Thus 

e Q = 1. 0-exp[-5 x 10 _5 (X e - 1000)] (3-32) 

where X fi is substituted with the dimensions lb^” 1 *. Equation 

*Refcr to the previous discussion on this point in Sub- 
section 2.4.1. 
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(3-32) suggests that a critical value for X g of 1000 lb^” ,s 

exists, belov/ which no entrainment of liquid occurs. The 

-k 

critical or limiting value for X g of 1000 lb^, was deter- 
mined partly from the experimental data of the present in- 
vestigation, but also from the experimental data due to 
Steen and Wallis (37, 38). They considered in some detail 
the problem of determining the critical flow conditions at 
which entrainment from a liquid film first occurs. As a 
result of their experiments, it was suggested that a critical 
gas stream velocity could be defined by the condition that 

ir c = 2.5 X 10" 4 (3-33) 
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The data reported by Steen et al. were for cold air flow, 

with most of the data obtained for water as the liquid phase. 

Assuming that the temperature of both the gas stream and the 

liquid film was approximately 60°F, and substituting the 

corresponding values for p for air and p, for water at that 

9 1 

temperature, Eq. (3-26) gives the following relationship 
between X g and it for the data of Steen et al. Thus 

X e = 3.8 x 10 6 tt (1 b f ~ is ) (3-37) 

Substitution of the critical value for tt that was proposed 
by Steen and Wallis into Eq. (3-37) results in the corre- 
sponding critical value for X fi of 950 lb^"* 5 . The critical 
value of 1000 lb^ ? ’’ ,s for X g in Eq. (3-32) was selected as a 
convenient compromise between the data of the present inves- 
tigation and those for the investigation of Steen et al. 

The experimental data of Steen and Wallis (37, 38) for 
the onset of entrainment, though obtained under much less 
severe flow conditions than those that characterized the 
present investigation, are in good agreement with the data 
of the subject investigation. It is Important to note, how- 
ever, that it was found that it was not possible to corre- 
late the experimental data of the present Investigation for 

e„ in terms of tt Instead of In terms of X„. Unlike the 
o e 

Investigation due to Steen et cl., the gas stream viscosity 
Pg-for the present investigation varied substantially. 


97 


f 

I 

f 

r 

i 

i 

i 

r 

f 

{ 

i 


Although the parameter n does have the advantage over the 
parameter X g of being dimensionless, the successful corre- 
lation of the data of the subject investigation suggests that 
the parameter X e better characterizes the phenomena. It will 
be argued in the following subsection that although the 
parameter X g is not itself a dimensionless variable, the 
correlation of e Q in terms of that parameter is equivalent 
to the correlation of e 0 in terms of a 'Mmensionless Weber 
number. 

It is interestir. . note from Fig.' 21 that for the 

present Investigation the entrainment parameter e Q varied 

from essentially zero to a maximum value of approximately 

0.6. The maximum value for e A that is shown in Fig. 21 

o 

corresponds to the experimental test 117 where methanol was 
employed as the film coolant. For that test, the computed 
simple rate of mass transfer was 0.022 Ib/sec-ft and the 
maximum mass transfer rate that was realized experimentally 
was 0.154 lb/sec-ft, a difference of 700 percent. Almost 
all of that difference (about 690 percent) was due to the 
entrainment phenomenon. In fact, it is Important to note 
that with the exception of a very small number of tests 
(most of which employed water as the film coolant which has 
a relatively high vaJue of surface tension), the contribution 
to the total rate of mass transfer that was due to the en- 
trainment of liquid into the gas stream was substantially 
larger than the contribution that was due to evaporation. 
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This was true considering both the simple rate of evaporative 
mass transfer and the evaporative mass trunsfer that was due 
to the phenomenon of interfacial roughness. 

3.1.4 The Physical Significance 
of the Correlation for e Q 

The physical significance of the satisfactory corre- 
lation of e Q in terms of the parameter X g is discussed in 
the present section. The discussion will suggest a procedure 
that is employed in Subsection 3.1. a to correlate the rough- 
ness parameter r. 

Since it is generally accepted that the phenomenon of 
entrainment is primarily the result of the shearing off of 
liquid droplets from disturbance waves on the liquid film 
surface, it is instructive to consider the forces that act 
on a wave that develops on the surface of the liquid film. 

If 1,. characterizes the scale of those interfacial waves, 
then the aerodynamic drag force that is exerted by the gas 
stream on the wave that acts to destroy the wave Is propor- 
tional to and the force that Is due to the surface 

tension that acts to keep the wave Intact Is proportional 
to ol . The balance between these two forces should be the 
controlling factor In tin breakdown of a disturbance wave on 
the surface of the liquid film. Therefore, that balance of 
forces should be Important In determining the rate at which 
liquid Is entrained from the liquid film Into the gas stream. 
Thus, If a dimensionless Weber number, denoted by We, 1st 
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defined as the ratio of those tv/o forces such that 


We 


Ml 
o s 


(3-38) 


it is reasonable to assume that the entrainment parameter e 0 
could be correlated as 


e 0 = e o (We) (3-39) 

The results of the foregoing subsection showed that 

e o ■ e 0 < X e> (3 - 4C) 

Upon comparing the defining equations for X and W , and 

v C 

ignoring for the sake of argument the term ( T q /T s ) ^* , it is 
seen that the apparent contradiction between Eq. (3-39) and 
Eq. (3-40) could be reconciled if a dimensional constant 
"K" existed such that 

K = 1 $ M^ *= constant (3-41) 

There is, in fact, some physical evidence presented herein 

that supports the assumption that the dimensional constant 

K does exist. Recall from the qualitative analysis of the 

photographic data in Section 2.4 of this report that the 

*The maximum value of this term for the subject investigation 
was 1.15. 


urn mm 
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parameter M Q (and not G or u g or any parameter for the 
liquid phase) was the parameter that most characterized the 
scale of the waves on the surface of the liquid film. More- 
over, the photographs showed that the scale decreased with 
increasing values for M q . Therefore, although the validity 
of the exponent of h in Eq. (3-41) cannot be argued on the 
basis of the qjalitative photographic data, there is some 
physical evidence to support the form of the equation. 

Another interesting observation can be made with ref- 
erence to Eq. (3-41) that gives some insight into the phe- 
nomena. It is well documented (18, p671) that the shatter- 
ing (or atomization) of liquid drops that are suspended in 
a gas stream can be characterized by a dimensionless Weber 
number, where the characteristic length (1 $ in Eq. (3-38)) 
is chosen as the radius of the drop. If that Weber number 
is greater than about 6, the liquid drop will break down 
into smaller drops. Now, if an analogy is drawn between the 
breakdown of such liquid drops and the breakdown of the 
liquid waves on the surface of the liquid film, then a crude 
calculation can be made for the constant K in Eq. (3-41). 
Assuming that the critical value for We for the case where 
1 characterizes the scale of the waves on the surface of 
the liquid film is still 6, and recalling that the critical 

valfue for the parameter X g below which no entrainment of 

_}* 

liquid occurred was approximately 1000 lb f - , then K can be 
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approximated by the following calculation. Thus 


6 

1000 


0.006 lb, 2 
f 


(3-42) 


Continuing this heuristic argument, it is interesting to use 
the above value for K, together with Eq. (3-41), to estimate 
the range of values for the scale 1 $ for the subject inves- 
tigation. From Table 6.1 of Appendix 6, the maximum and the 
minimum experimental value for the momentum parameter M Q was 
1087 lb^/ft^ and 20 lb^/ft^, respectively. Thus, employing 

j, 

the value for K of 0.006 lb^ , the maximum and the minimum 
value for the scale 1 for the subject investigation can be 
estimated from Eq. (3-41). Thus 


s ,max 


0.00135 ft 


1 . = 0.00018 ft 

s , m l n 

While these calculations are admittedly very approximate, 

the values obtained for the scale 1 are surprisingly similar 

to the dimensions that characterize the structure of a liquid 

film that is subjected to relatively large interfacial shear 

forces (as is the case in the present investigation). For 

example, Ref. (19) shows that the mean thickness of such a 

-4 

liquid film is of the order of 10 ft. Moreover, the wave 
length of the small scale disturbances can be roughly approxi- 
mated as 10 times the mean thickness of the liquid film (8, 5) 

_ 3 

or the orde^ of 10 ft, and the amplitude of those distur- 
bance waves is of the same order as the thickness of the 
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liquid film (5), or 10”^ ft. The preceding calculations, 
therefore, add some additional support to the assumption 
that the dimensional constant K exists. Moreover, they also 
suggest that the breakdown of the interfacial waves and the 
resultant entrainment of liquid into the gas stream may be 
somewhat analogous to the atomization of liquid drops that 
are suspended in a gas stream. 

To summarize briefly, the above discussion suggests that 
(a) the dimensional parameter X g and the dimensionless Weber 
number We are physically equivalent parameters, and (b) a 

j* 

constant K = 1 M „ 2 exists such that the characteristics 
s o 

(really, the gross characteristics) of the interfacial film 
structure can be characterized in terms of the flow param- 
eters for the gas stream. The second of these two points, 
which follows directly from the first, is most pertinent to 
the correlation of the roughness parameter r that is pre- 
sented in the following section. 

3.1.5 The Correlation of the 
Roughness Parameter r 

The present section is concerned with the correlation 
of the roughness parameter r that was defined by Eqs. (3-4) 
and (3-6). T he primary assumption that was employed to 
develop the correlation was that, at least to a first ap- 
proximation, the degree to which the effective roughness of 
.the surface of the liquid film intensifies the turbulent 
transfer of energy from the hot gas stream to the surface 




W 
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of the liquid film is primarily a function of the scale of 
the small disturbance waves that develop on the surface of 
the film. That is, it was assumed that the roughness param- 
eter r could be correlated in the form 

r = r ( 1 s ) (3-43) 

Unf ortunately , the utility of Eq. (3-43) is limited because 
of the problems inherent in trying to measure the scale 1 . 
Indeed, there is considerable uncertainty as to what physical 
dimension for the waves on the surface of the liquid film 
corresponds to the length 1 $ . However, an argument was pre- 
sented in the foregoing subsection that the scale 1 $ could be 
approximately related to the momentum parameter M by the 
following expression. Thus 

*A (3 - 44) 

H o 

Therefore, if a roughness group "X r " is defined by 

T ^ 

x r = (yfl.) (3-45)* 

Equation (3-43) can be rewritten in the following equivalent 
form. Thus 

*Again, as for Eq. (3-33) in Subsection 3.1.3, the term 
(T /T,)' is an empirical correction for the dependence of • 
the density of the gas in the boundary layer on temperature. 
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r = r(X r ) (3-46) 

Equation (3-46) is the basic expression that was employed in 
the subject investigation to correlate the roughness param- 
eter r. The value of the roughness parameter r was estimated 
for each experimental test as follows: 

(a) The average value for the simple rate of mass transfer, 
m^, was calculated following the procedure outlined in Sub- 
section 3.1.2. 

(b) The value of the entrainment parameter e Q was determined 
from the average value of the entrainment group X g for that 
test, together with the analytical expression given by Eq. 
(3-32). 

(c) The value of the intercept "A" (see Fig. 19) was esti- 
mated by determining which curve having the slope (l-e o ) 
resulted in the best "eyeball" fit of the experimental data. 

(d) The roughness parameter r was then determined from the 
expression 

r ■ j-j-r - 1.0 (3-47) 

5 

Figure 22 presents the roughness parameter r--determined 
by the aforementioned procedure--as a function of the param- 
eter X r> The various conditions for the gas stream pressure 
and temperature are listed in the figure. Although the data 


H 

B 
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in Fig. 22 do exhibit some scatter, much of the scatter can 
be attributed to the fact that the value for the intercept 
A was very sensitive to relatively small errors in the mea- 
surement of the net rate of mass transfer from the liquid 
film to the gas stream. Furthermore, any error in determin- 
ing the value of the intercept A was compounded in the calcu- 
lation of the roughness parameter r. To illustrate, assume 
that the true value for the parameter A/m^ is 1.2 and that 
the experimentally determinec >.\Jue is 1.4. From Eq. (3-47), 
the corresponding values for the roughness parameter r are 
0.2 and 0.4, respectively. Thus, for this example, an error 
of approximately 17 percent in the intercept A results in an 
error of 5* percent in the evaluation of r. In addition to 
that problem, some of the scatter of the data in Fig. 22 is 
probably due to the fact that the rather complex phenomenon 
of interfacial roughness cannot be characterized completely 
in terms of the relatively simple parameter X f . .However, 
considering the wide range of values for the liquid and the 
gas flow parameters that is represented by the data in Fig. 
22, the correlation of r in terms of the group X r would 
appear to have some physical justification. 

The solid curve that Is shown in Fig. 22 is given by 
the following analytical expression. Thus 


3.0 

0.8 



r 


(3-48) 
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\c 

where X p has dimensions of 1 b ^ / f t . While the experimental 
data in Fig. 22 are in general agreement with this analytical 
expression, the final form of Eq. (3-48) was arrived at by 
considering the following three points, which are listed in 
their order of importance: 

(a) The correlation of all of the experimental data of the 
present investigation, including those for which the gas 
stream was accelerated over the wetted test plate; 

(b) The correlation of the mass transfer data that were 
reported by Kinney, Abramson, and Sloop (19) (these data are 
analyzed in Section 4.3); and 

(c) The correlation of the experimental data for mass trans- 
fer that are reported by Emmons (20) (these data are ana- 
lyzed in Section 4.4). 

It is important to note that Eq.' (3-48) and the data in 
Fig. 22 suggest that the interfacial roughness can contribute 
substantially to the convective transport of energy from the 
hot gas stream to the surface of the liquid film. Indeed, 
Fig. 22 suggests that the simple mass transfer rate m^ could 
be increased by 100 percent (r B 1), or more, due to the 
interfacial roughness. In comparison, for the case of heat 
transfer to a dry rough wall. Fig. 7-56 of Ref. (39) shows 
that the value for r can be as great as 2.5. It would seem, 
therefore, that a value of r of 1.0 for the problem of 
liquid-film cooling is not unreasonable. It should be em- 
phasized, however, that in the analysis of the experimental 
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data of the present investigation and in the analysis of the 
data of the primary references (19) and (20), that the con- 
tribution to the net mass transfer rate that was due to the 
phenomenon of interfacial roughness was, with very few excep- 
tions, substantially less than the contribution that was due 
to the entrainment phenomenon. Frequently, the mass transfer 
due to entrainment was 10 to 20 times the evaporative mass 
transfer that was attributed to the effective roughness of 
the surface of the liquid film. 

3.1.6 The comparison of the Predicted and 
the Measured Rates of Mass Transfer 

A comparison is presented in this subsection of the 
rates of mass transfer from the liquid film to a constant ve- 
locity gas stream that were predicted from the correlation de- 
veloped herein to those that were determined experimentally. 
The comparison is presented for both the datum points that 
were obtained for the hot-flow experiments (approximately 
500 datum points) and those that were obtained for the cold- 
flow experiments (76 datum points). In addition, to empha- 
size the utility of the present theory, the experimental 
data are also compared with the results that were predicted 
from the simple theory of mass transfer. The latter com- 
parison serves to Illustrate graphically the general inade*- 
quacy of the simple theory tor the subject problem. 

The primary equation that was employed to calculate the 
rates of mass transfer from the liquid film to the gas 
stream is repeated here for convenient. Thus 
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m* = (1 + r) nig + e Q m (3-6) 

The simple mass transfer m^ was calculated following the de- 
velopments of Subsection 3.1.2. The roughness parameter r 
was calculated from Eq. (3-48) and the e.itrainment parameter 
e Q was determined from Eq. (3-32). The surface tension o 
that appears in the entrainment group X g in Eq. (3-32) was 
based on the average of the temperature of the liquid at the 
point of injection and at the point of withdrawal*. All of 
the numerical computations were performed on an IBM 7094 
digital computer. 

Figures 23 through 26 summarize graphically the results 
of the aforementioned calculations for those experimental da- 
ta where the gas stream was heated. Each figure represents 
one of the four liquids that was investigated (Fig. 23: 
methanol; Fig. 24: butanol; Fig. 25: water; Fig. 26: RP-l). 
The nominal value , for the gas stream pressure and the gas 
stream temperature for each set of experimental data are 
indicated in the figures. Part (a) of each figure is a 
plot of the simple mass transfer rate, m^, that was pre- 
dicted from the theory outlined ! n Subsection 3.1.2 ver- 
sus the measured rate of mass transfer. Part (b) of 
each figure compares the measured rate of mass transfer 
with the rate of mass transfer that was predicted from 

♦See Table G.l of Appendix G for a tabulation of the average 
inlet and outlet liquid temperatures for. each experimental ' 
test. 
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Eq. (3-6) for the present theory. Each figure graphically 
demonstrates that, in general, the simple theory for’mass 
transfer is a poor approximation to the complex phenomena 
that characterize liquid-film cooling. The maximum differ- 
ence between the experimental result and the result that was 
predicted from the simple theory is shown in Fig. 23(a) for 
methanol as the film coolant. For that figure, the maximum 
value for the measured mass transfer rate is 0.154 lb/sec-ft 
and the corresponding simple rate of mass transfer is 0.022 
lb/sec-ft, a difference of 700 percent. The present theory, 
however, accounts for the existence of the interfacial phe- 
nomena, and the resultant agreement between the predicted and 
the measured mass transfer rates that is shown in part (b) of 
each figure is quite satisfactory. The maximum difference 
between the result that is predicted from the present theory 
and the experimental result is shown in Fig. 26(b) for RP-1 
as the film coolant. For a measured rate of mass transfer of 
0.0944 lb/sec-ft in that figure, the corresponding mass 
transfer rate that is predicted by the present theory is 
0.0765 lb/sec-ft, a difference of about 20 percent. 

Figure 27 presents the results of the computations that 
were performed for the case where the gas stream was not 
heated. Presented in the figure is a plot of the mass 
transfer rate that was predicted by the present theory 
versus the measured rate of mass transfer. The experi- 
mental data in the figure are for one nominal gas stream 
pressure of 75 psia a.nd for the liquids methanol, butanol, 
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and water*. It was assumed in the computations that the 
rate of simple mass transfer, mj , was essentially zero be- 
cause of the low gas stream temperatures. Thus, the realized 
rates of mass transfer were assumed to be due to entrainment 
only. The maximum value that was computed for the entrain- 
ment parameter e Q was 0.32, as compared to a maximum value 
of 0.63 for the hot-flow tests. The agreement between the 
measured and the predicted rates of mass transfer is quite 
satisfactory, particularly in view of the fact that the 
pertinent physical properties for both the liquid and the 
gas phases differed so significantly between the case where 
the gas stream was heated and the case where it was not. 

For example, the viscosity for the liquid varied by as much 
as a factor of 15:1, the surface tension for the liquid by 
as much as a factor of 1.5:1, and the viscosity for the gas 
stream by as much as a factor of 2:1 between the hot and 
cold gas stream flow conditions. It would appear, therefore, 
that the relatively simple theory that is proposed herein 
for determining the rate of mass transfer from a thin liquid 
film to a high velocity gas stream can be applied quite 
generally. The utility of the present theory is demonstrated 
further in Section 4 of this report where the correlation of 

*The number of datum points that could be obtained with water 
as the film coolant wc: limited because the gas stream was 
generally cooled sufficiently by the expansion from the 
storage tank pressure of 1000 to 2500 psia to the test sec- 
tion pressure cf 75 psia that a portion of the water that 
was Injected onto the test plate froze before it reached 
the capture slot. Moreover, no e xoeri mental tests were 
performed with RP-1 as the liquid phase because of the 
limited supply of th.is liquid. 
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the experimental data for liquid-film cooling ot the pri- 
mary references (19, 20, 24) is considered. The gas stream 
flow conditions that characterized those investigations will 
be shown to be, in most instances, much more severe than 
those that characterized the present investigation. 

The analyses of the experimental data of the afore- 
mentioned primary references requires first the consider- 
ation of the problem of how the basic correlation for mass 
transfer can be extended to the more general case where the 
film cooled length 1 is arbitrary. That is, the case where 1 
does not necessarily equal 1 Q = 10 inches. The following 
section presents an analysis of that problem. The result of 
the analysis will also suggest a method for correlating the 
experimental data for mass transfer that were presented in 
Section 2.6 of the present report for the case where the gas 
stream was accelerated over the film cooled length 1 Q . The 
correlation of those data is then discussed in Section 3.3. 

3.2 The Modification of the Basic Correlation 
for Mass Transfer to Account 
for an Arbitrary film Cooled Length 

The correlation that was developed for the entrainment 
parameter e Q in Subsection 3.1.3 is necessarily dependent on 
1 0 , the length of the film cooled test section for the sub- 
ject investigation. That is. if the film cooled length 1 0 

i 

had been, say, greater than the 10-1n. length that was em- 
ployed, then, for the same flow conditions, there would have 
been more entrainment of liquid into the gas stream because 


V 
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of the prolonged exposure of the surface of the liquid film 
to the gas stream shear forces. In order to be able to pre- 
dict the mass transfer that occurs for an arbitrary film 
cooled length 1, the basic correlation for e Q has to be mod- 
ified. A rather simple model and analysis is developed in 
iubsection 3.2.1 to solve that problem. Subsection 3.2.2 
then suggests an empirical approximation to the result of 
that analysis that is generally more convenient to employ. 

3.2.1 Model and Analysis 

Figure 28 presents the simple model that is employed 
in the analysis. The pertinent notation for the figure is 
defined as follows: 

x' * x - 

* the rate of liquid injection at x* = o 

1 = the liquid-film cooled length 

1 * the length of each film element* 

n ■ 1/1 (assumed to be an integer) 

mj* = the net rate of mass transfer from the 
K element of the liquid film, k * 1, 2, 3,...n 

m' . = the net rate of simple mass transfer for the 

S ' K k th film element, k ® 1, 2, 3,...n 

*For the purposes of the. present analysis. It is not nec- 
essary to specify the* 1 * 10 Inches. However *« "ti- 

lize the result of the analysis to perform any numerical 
calculations, it will be necessary to make that speci- 
fication. 
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The basic assumptions employed in the analysis are as 
f ol 1 ow: : 

(a) The flow parameters for the gas stream are independent 
of x' . 

(b) The roughness parameter r is independent of x'. 

(c) The local rate of simple mass transfer varies with the 
coordinate x' as 

m" = c x' p (3-50) 

where c and p are appropriate constants. 

(d) The primary equation (3-6) applies to each of the 
elements of the liquid film. 

(e) The liquid film terminates when 1 ■ n 1 0 , where n is an 
1 nteger. 

Thus, for the first element of the liquid film, the net 
rate of mass transfer Is given by 

c (1 + r) ir^ j + e Q (m 1 ) (3-51) 

Simllarily, the mass transfer for the second film element is 
given by 

■ (1 + r) m^ 2 + e () (m 1 - mp (3-52) 

In general, for the k th element of the liquid film, the net 
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rate of mass transfer m£ is expressed as 

m£ = (1 + r) m^ + e^n^-mj -mj -...-mj^) (3-53) 

Substitution of Eq. (3-51) into Eq. (3-52) gives the 
result 


"i ' * r > »;,2 ♦ e 0 t"i 

Substitution of Eqs. (3-51) and 
k = 3 yields 

*3 - (1 * r )*5 ( 3 e„ {« x 

Vl - (1 * r) "s,2 ‘ e o C " 


- (i * r) in* e^] (3-54) 
(3-54) into Eq. (3-53) with 

‘ (1 + r > m i,l “ 


e^^]} . (3-55) 

This substitution procedure can be continued in a stepwise 
fashion so that for the general kth term m£, the terms mj, 
m^, . .. m k-x can be eliminated from the resultant expression. 
Moreover, after considerable algebra and a fortuitous com- 
bining of terms, it is possible to show that 

"i ■ "a - ♦ <1 ♦ r> «; (k (l-e 0 )"- k 


(3-56) 
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But, if it is assumed that the iiquid in the film is consumed 
at a point where 1 = n 1 Q , then it follows from the conser- 
vation of the liquid that 

n 

l ni = m, (3-57) 

k=l K 1 

Equation (3-56), therefore, can be rewritten in the form 

n t 

m l = (1 + r ) I m i,k (l ‘ e o ) (3-58) 

From assumption (c) of the model it follows that the 
rate of simple mass transfer from the k th film element is 
given by 




(k-Di 0 




dx' 


<(kl 0 ) ptl -[(k-l)l 0 ] pn } 
(P + 1) 


(3-59) 


Moreover, the total rate of simple mass transfer from the 
liquid film is equal to 


■; 



(nlj p+1 

(p+i) 


(3-60) 
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Finally, substituting Eq. (3-59) into Eq. (3-58) and dividing 
through both sides of the resultant expression by Eq. (3-60) 
yields the following non-dimensional equation. Thus 



(1 + r) I [(£) 


k=l 


p+1 , , p+1 _. 

- <*=*> 3U-e 0 ) k 


(3-61)* 


Equation (3-61) is the desired analytical result. For 
given gas stream and liquid flow conditions (and, therefore, 
given values for the parameters p, r and e Q ), Eq. (3-61) can 
be utilized to calculate the net rate of mass transfer from 
the liquid film into a constant velocity gas stream if the 
film cooled length 1 is equal to n 1 0 , where n is an integer. 
The obvious disadvantage of Eq. (3-61) is the fact that it 
can be utilized for only those special cases where 1 = n 1 Q . 
Therefore, it is convenient to replace that, analytical re- 
sult by an approximate empirical expression that is a con- 
tinuous function of 1, where 1 can be either greater than or 
less than 1 Q . Such an empirical expression is developed in 
the following subsection. 

♦Note that for the special case where the Interfacial 
phenomena are not Important, 1 .e. , r » e Q * o, then 

the RHS of Eq. (3-611 reduces to unity as required. 
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3.2.2 An Empirical Approximation 
to the Analytical Result 

An empirical approximation to the analytical equation 
(3-51) can be derived as follows. First note that for the 
special case where 1 = 1 Q , Eq. (3-61) reduces to the simple 
expressi on 



(1 + r) (l-e 0 )" 1 


(3-62) 


Now assume that Eq. (3-62) can be modified as follows for 
the case where the film cooled length is arbitrary. Thus 



= (1 + r) (l-e 0 ) 


-0/1/ 


(3-63) 


The parameter b in Eq. (3-63) is determined by first speci- 
fying a value for the parameter p that appears in the ana- 
lytical expression (3-61), and then "matching" the empirical 
expression (3-63) and the analytical expression (3-61) at 
those points where ( 1 / 1 Q ) = an integer. 

To illustrate the aforementioned procedure for the 
evaluation of the parameter b, consider the classical case 
of the turbulent transport of heat and mass from a flat 
plate where the parameter p in Eq. (3-50) is equal to -1/5. 
Figure 29 presents a plot of the quantity (m^/m^)/ (1 + r) 
as a function of the ratio (1/1 Q ) with the entrainment pa- 
rameter e Q as a parameter. The solid curves shown In Fig. 29 
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were calculated from Eq. (3-63) with a value for b of 0.6, 
and the circles are the "match" points that were calculated 
from the analytical expression (3-61). 

Figure 30 presents another illustrative example for the 
case where the parameter p in Eq. (3-50) is equal to 1.0; 
that is, the case where the rate of simple mass transfer 
increases linearly with the coordinate x'*. The solid curves 
were calculated from Eq. (3-63) where the value of 0.9 for 
the parameter b was employed, and the "match" points were 
calculated from Eq. (3-61). 

Figures 29 and 30 show that the analytical result of 
Subsection 3.2.1 can be approximated with sufficient accuracy 
by the relatively simple empirical formula of the present 
section. Moreover, these figures show that a substantial 
change in the parameter p in Eq. (3-61) requires a relatively 
small change in the parameter b in the empirical expression 
(3-63) in order to obtain an adequate matching of the ana- 
lytical and empirical results. 

The significance of the foregoing developments can be 

further illustrated by writing Eq. (3-63) in a more explicit 

form. Substituting the analytical expression (3-48) for 

the roughness parameter r and the analytical expression 

(3-32) for the entrainment parameter e 0 , Eq. (3-63) can be 

♦This example is presented solely for the purposes of 
illustration. It may or may not haye any physical 
significance. 
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written in the form 



(1 + -^g) exp [5 x IQ' 5 (X e - 1000) (l/1 0 ) b ] 

X r 

(3-64)* 


Thus, if the value for the parameter b is approximated from 
the aforementioned "matching" procedure, Eq. (3-64) gives 
an explicit expression that can be employed to calculate the 
net rate of mass transfer from a liquid film of arbitrary 
length 1 to a constant velocity gas stream. 

The validity of the foregoing analysis (and, therefore, 

Eq. (3-64)) will be investigated in Section 4 of the present 
report. Considered in that section Is the correlation of 
the experimental data for mass transfer that were reported 
in the primary references (19, 20). The range of values for 
the film cooled length 1 that Is represented by those data 
is 3 inches to 36 inches long. 

Equation (3-64) also serves to suggest a method for the 

correlation of the experimental data for mass transfer that 

were obtained for the case where the gas stream was accelerated 

*So that the reader Is not confused, the substitution of Eq. 
(3-32) for the entrainment parameter e 0 requires that the 

length 1 0 in Eq. (3-64) be taken as 10 Inches. Moreover, 

tne numerical constants In Eq. (3-64) require that the 

roughness group X f be substituted with the units (Ib^/ft) 

. and the entrainment group X. be substituted with the units 

Obf**). 


129 


over the film cooled length 1 Q . Those data were presented 
in Section 2.6 of the present report. The correlation of 
the subject data Is considered In the following section. 

3.3 Correlation of the Data for Mass Transfer 
for ah Accelerating Gas Stream 

The purpose of the present section Is to correlate the 
experimental data for mass transfer that were obtained for 
the case where the gas stream was accelerated over the length 
of the film cooled test plate. These data were shown graphi- 
cally In Fig. 17 of Section 2.6. Subsection 3.3.1, which 
follows, shows how the simple rate of mass transfer was com- 
puted for the subject case. Subsection 3.3.2 then demon- 
strates how that simple rate of mass transfer was employed, 
together with a logical extension of the developments that 
were presented in Sections 3.1 and 3.2, to correlate the 
experimental data for an accelerating gas stream. The com- 
parison between the experimental and the predicted net mass 
transfer rates Is presented In Subsection 3.3.3. 

3.3.1 The Evaluation of the Simple Rate of 
Mass Transfer for an Accelerating Gas Flow 

Kays (36) suggests that the following expression can 
be employed to evaluate the local rate of simple mass trans- 
fer for the case where the gas stream is accelerating. Thus 
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■; * G Wc F nT F Pr C°- 0295 Re i"° 8 Re i t ‘°' 12 3 

(3-65)* 


where 






x 


I 


0 


I 


G dx 


G dx 


(3-66) 


(3-67) 


The "F" correction factors In Eq. (3-65) are given by Eqs. 
(3-18) - (3-21) of Subsection 3.1.2. For the purposes of 
the present calculations, the energy balance parameter <f> 
that enters Into the definition of the parameter B^** was 
approximated from Eq. (3-13) of Subsection 3.1.2. Moreover, 
because of the limited Mach number range that was Investi- 
gated, the variation In the gas stream temperature, T g , over 
the film cooled length 1 Q was so small that the average of 


* Kays (36) Introduces an additional empirical correction 
factor Into Eq. (3-64) to correct for acceleration effects 
that are not predicted analytically. However, preliminary 
calculation of that correction factor showed that It was 
wf the order of 0.98 for the range of flow cunniilons for 
the subject Investigation. Therefore, that correction 
factor was not considered In the present calculations for 


**See Eq. (3-22). 
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the values of at x = and at x s was employed to 
calculate the parameters , B^, and Fj in Eq. (3-65). To 
complete the solution for m£ , it is necessary to determine 
the local values for the static pressure p, the g«s stream 
Mach number M a , the gas stream velocity u^, and the mass 
velocity for the gas stream 6 . The procedure employed to 
evaluate those parameters Is outlined as follows. The static 
pressure was measured at the location x = 20 Inches in the 
approach section and at 5 axial locations (equlspaced at 
2 Inches) in the test section by wall taps that were located 
In the roof of the experimental tunnel. The stagnation 
pressure for the gas. flow was calculated from the static 
pressure that was measured In the approach section, together 
with the value that was computed for the gas stream velocity 
at the entrance to the test section*. The local value for 
the Mach number within the experimental test section was 
computed by forming the ratio of the measured static pressure 
and the calculated stagnation pressure, and then employing 
the following relationship for Isentroplc flow with a 
specific heat ratio y ■ 1.385. Thus 

*5 * I 1 + ¥ (3-68) 


*See Appendix F.2 for the procedure that was employed to 
calculate the. gas stream velocity at the entrance to the 
test section. 
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The complete profile for the Mach number throughout the 
test section was obtained by fitting a third-order poly- 
nominal to the local values for the Mach number that were 
computed from the experimental measurements for the static 
pressure. The resultant form or the approximate expression 
for the Mach number as a function of x', the distance from 
the point of liquid injection, was 

M a = 0.145 + 0. 02x' - 0.0014X' 2 + 0.00028x' 3 (3-69) 

The local values for the gas stream velocity, u^, were com- 
puted from Eq. (3-69) and the equation for the acoustic 
speed In a perfect gas: 

u g ■ h . y 0-5 < 3 - 70 > 

Finally, the local values for the gas mass velocity, 6, were 
calculated from the defining equation 


Pg Ug 


( 2 - 1 ) 


where the gas stream density, p g(< was evaluated from the 
equation of state for a perfect gas: 


p g ETT" 
9 9 


(3-71) 
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Figure 31 shows graphically the values that were computed 
from the foregoing procedure for the parameters p/p°, M . 

a 

and u„ as a function of x'. The solid curves that are pre- 
9 

sented in the figure are based on the analytical expression 
(3-69) for the Mach number. The circles presented in Fig. 31 
represent the static-stagnation pressure ratios that were 
realized experimentally. Each circle presented in the fig- 
ure represents the approximate range of values for (p/p°) 
that was obtained at that particular location for the entire 
spectrum of gas flow conditions that was investigated. 

The foregoing equations were programmed for an IBM 7094 
digital computer and the integral expressions for the terms 
Re^ and Re^ t were evaluated numerically. 

3.3.2 The Model and Analysis 
Figure 32 illustrates the model that was employed in 
the analysis of the subject data. The figure shows that the 
liquid film was subdivided into a number of incremental ele- 
ments of length <5x. The pertinent notation that applies to 
the general kth film element is as follows: 

(a) x fc , = (k-l)6x = the distance from the point of 

liquid injection to the beginning of 
the kth film element, k=l,2,3, . . .n. 

(b) u.,p. .... * the value for the indicated flow 

K K parameter for the kt_h film element 

evaluated at x^ - <sVx, k*l,2,3,. . .n. 

(c) X X . * the value for the indicated dimen- 

* ,K sional group for the kth film ele- 

ment evaluated at x k -“^5x, 
k-1,2,3, . . .n. 
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Fig 32. Model for Analysis of Experimental Data for Mass 
Transfer. Accelerating Gas Stream 
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(d) m k = the rate of liquid flow into the kth 

K film element, k=l , 2 ,3 , . . . n . 

(e) m k = the net rate of mass transfer from the 

K k th film element, k=l,2,3, . . .n. 

(f) m' . = the rate of simple mass transfer from 

s * the k th film element, k=l,2,3 , . . . n . 

The primary assumption that was employed in the analysis 
is that the net rate of mass transfer from the k th film 
element, m‘, , is given by the expression 

m k = U + r |<) m s,k + e k m k (3-72)* 

where 


r 


k 


3.0 

y 0 • 8 

A r,k 


(3-73) 


and 


e k = 1.0 - exp[- 5 x 10‘ 5 (X e#k - 1000) (|^) ] 


(3-74)** 


The value for the parameter d in Eq. (3-74) was determined 
from the experimental data. 

The procedure that was employed to calculate the net 
rate of mass transfer from the liquid film to the acceler- 
ating gas stream for each datum point can be outlined as 

* Note that Eq. (3-72) is analogous to Eq. (3-6). 

**flote that Eq. (3-74) is motivated by the developments of 
Subsection 3.2.2. 
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follows. 

(a) The pertinent gas stream flow parameters were calculated 
for each location following the procedure outlined in 
Subsection 3.3.1. 

(b) The entrainment group, X g k , and the roughness group, 

X r k> were calculated for each location from the flow 
parameters that were determined in step (a). 

(c) The rate of simple mass transfer, m^ k , for each film 
element was calculated from the equation 


m 


s.k 


= m 


s , k 


6x 


(3-75) 


where m” k was computed at the location x k - *s6x from the 
procedure that was outlined in Subsection 3.3.1. 

(d) A value for the parameter d in Eq. (3-74) was selected. 

(e) The net rate of mass transfer for the first film ele- 
ment (i.e., k = 1) was calculated from 


m i “ t 1 + r i) m s,l + e l m l (3-76) 

where m 1 was the experimentally measured rate of liquid 
1 njection. 

(f) The net rate of mass transfer from the second film 
element (i.e., k » 2) was calculated from 

= (l + r^) Rig 2 * e 2 m 2 


( 3 - 77 ) 
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where 


m 2 8 (3-78) 

(g) Steps (e) and (f) were continued in a stepwise fashion 
so that for the general k th film element the net rate of 
mass transfer was calculated from 

m k * U + r j ( ) + e k (3-79) 

where 


m k 




(3-80) 


(h) When the net rate of mass transfer for each film ele- 
ment had been calculated, the net rate of mass transfer from 
the entire liquid film was calculated from 





where 


n 



(3-81) 


(.3-82) 


(i) The result of step (h) was compared with the experi- 
mental result for m' ( = i, - I* those two results 


I 

I 

I 

I 

I 
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were not in reasonable agreement, then a new value for the 
parameter d was selected and steps (e) through (h) were 
repeated until a satisfactory agreement of those results was 
obtai ned. 

The aforementioned procedure was employed to analyze 
each of the 50 datum points that was shown in Fig. 17 of 
Section 2.6. The computations were performed for values of 
6x of 0.1, 0.2, and 0.4 inches to see if that parameter 
significantly influenced the results. The calculations 
showed that the simple value of 1.0 for the parameter p in 
Eq. (3-74) resulted in a good correlation of all of the 
datum points for each of the aforementioned values for the 
parameter 6x. 

3.3.3 The Comparison of the Predicted and the 
Measured Rates of Mass Transfer 

Figure 33 presents the comparison of the experimental 
data for the net rate of mass transfer from the liquid film 
to the accelerating gas stream with the results that were 
predicted from both (a) the simple theory for mass transfer 
that was outlined in Subsection 3.3.1, and (b) the present 
theory for mass transfer that was developed in Subsection 
3.3.2. The stagnation pressure, p°, and the film coolant 
that is represented by each set of experimental data is 
noted in the figure. 

Figure 33(a) again graphically Illustrates that the 
simple theory for mass transfer is, in general, a very poor 
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approximation to the complex phenomena that characterize 
liquid-film cooling. The maximum difference between the 
measured and the analytical result that is shown In the 
figure occurs for methanol as the liquid phase at a mea- 
sured rate of mass transfer of 0.1668 lb/sec-ft and a 
corresponding predicted value of 0.0309 lb/sec-ft, a dif- 
ference of more than 500 percent. 

Figure 33(b) shows that the present theory, with a 
value of 1.0 for the parameter d in Eq. (3-74), resulted in 
a satisfactory correlation of the experimental data. The 
maximum difference between the measured and the predicted 
result occurred for water as the film coolant at an experi- 
mental value for the net mass transfer rate of 0.0462 
lb/sec-ft and a corresponding predicted value of 0.0620 
lb/sec-ft, a difference of less than 40 percent. 

The utility of the present analysis is investigated in 
further detail in Section 4. Considered In that section is 
the correlation of the data that were reported in the pri- 
mary reference (24) for the liquid-film cooling of a Mach 
number 1.92, axi-symmetric nozzle. Also presented in that 
section Is an illustrative example of the coolant require- 
ments that are predicted by the present theory to cool an 
exhaust nozzle for which the gas stream has a stagnation 
temperature of 4000R and a stagnation pressure of 500 psla. 
These stagnation conditions were chosen as being typical of 
those that characterize the gas flow In the exhaust nozzle 
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of a high-energy propulsion device. 

3.4 Correlation of the Experimental Data 
for the Maximum Liquid Temperature 

The correlation of the experimental data that were 
reported in Section 2.7 for the maximum liquid temperature, 

T-j m » is considered in this section. Subsection 3.4.1 pre- 
sents a simple theory for calculating that liquid temperature. 
The comparison of the theoretical and the experimental re- 
sults is shown in Subsection 3.4.2. 

3.4.1 A Method for Predicting the Maximum 
Liquid Film Temperature* 

The maximum liquid-film temperature, T-j m> can be 

related to the bulk flow parameters for the gas stream by 

specifying (a) the appropriate energy balance across the 

gas-liquid interface, and (b) an analogy between the rates 

of convective** heat and mass transfer. The discussion of 

Section 2.7 showed that the maximum liquid temperature is 

obtained when all of the energy transfer to the liquid film 

results in vaporization of liquid in the film, with none of 

the energy transfer resulting in an increase In the sensible 

enthalpy for the liquid in the film. It follows from the 

discussion of the interfacial energy balance In Subsection 

3.1.2 that the maximum liquid temperature Is obtained when 

* This method is essentially the same as that previously 
presented by the author In Refs. (25) and (26). 

**Mass transfer due to entrainment Is not considered. 
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the energy balance parameter <J> at the surface of the liquid 
film reduces to 


4 > « H Vj|n (3-83) 

where H„ „ Is the heat of vaporization for the liquid eval- 
v ,m r 

uated at the maximum liquid temperature T^ m . 

The simplest analogy between convective heat and mass 
transfer Is that given by (36) 

St s «= St; (3-84) 


where St! is the Stanton number for simple mass transfer 
& 

defined by 


St 


i 

s 




(3-85)* 


The term p y m In Eq. (3-85) Is the partial pressure for the 
vapor at the surface of the liquid film evaluated at the 
temperature T, 

I 

It Is Important to note that this analogy, although 

very simple, should be approximately correct for quite 

general flow conditions. Including the case of accelerating 

*Equation (3-85) can also be written In terms of the con- 
centration of the vapor at the surface of the liquid film, 
C . Thus 



(3-86) 
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gas flow and the case where considerable entrainment of 
liquid into the gus stream is occurring. In other words, 
the acceleration of l .he gas stream and/or the entrainment of 
liquid should have a similar influence on the development of 
both the thermal and the vapor-concentration boundary layers 
over the surface of the liquid film. 

Combining Eq. (3-15) with Eqs. (3-83) - (3-85) yields 


P - P 


P r 1 

w v,s V 


v,s 




(3-87a )* 


H 

TC^nr 


Vjjn 


'pg 


1 ,m J 


( 3- 87 b ) 


Equation (3-87b) is the desired analytical expression. 
Since a value for the liquid temperature T-j m uniquely 
determines the value for both the vapor pressure p y m and 
the heat of vaporization H y m . It follows from Eq. (3-87b) 
that the maximum liquid temperature is expressed solely in 


♦Note that if Eqs. (3-85) and (3-86) are combined to elimi- 
nate the left-hand side of Eq. (3-87a), and if the param- 
eter 6^ Is Introduced, then the resultant expression can 

be written In the form 


B. 


1 + B, 


(3-26) 


Equation (3-26) was Introduced In Subsection 3.1.2 
without derivation. 
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terms of the flow parameters p, T , c , and M for the gas 

y r y y 

stream. 

Figure 34 presents, for purposes of illustration, 
calculated values of the maximum liquid temperature T-j m for 
an air-water system*. Shown in the figure is a plot of T-j m 
as a function of the static pressure p with the gas stream 
temperature T g as a parameter. The saturation curve for 
water that is presented in the figure illustrates the 
approximate error that is introduced by the assumption that 
the maximum film temperature is equal to the saturation 
temperature for the liquid coolant at the local static 
pressure, p. That assumption was employed, for example, 
by Refs. (20, 30, 23, 42). References (25), (26) and (28) 
suggest that the calculated results shown in Fig. 34 are in 
good agreement with the reported experimental data for the 
film temperature T-j m for an air-water system. 

3.4.2 Comparison of the Experimental and the Predicted 
Results for the Maximum Film Temperature 

Section 2.7 previously discussed the nature of the 

experimental measurements for the maximum liquid temperature 

T, **. Presented in Fig. 18 of that section was a comparison 
I *m 

* The calculations were performed by employing Eq. ( 3- 87 b ) 
in conjunction with the tabulated data of Refs. (40) and 
(41). A sample calculation Is presented Appendix F. 

♦♦Recall from Section 2.7 that the maximum liquid tempera- 
ture Tj was determined experimentally by measuring the 

temperature for the liquid at the point of liquid with- 
drawal , T, 9 . 



Tj « maxlmun liquid temperature (R) 
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Fig. 34. Computed Results for Maximum Liquid 
Temperature as a Function of T g and p 

for an Air-Water System 
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between the temperature m and the boiling temperature for 
the various liquids* of the subject investigation at the pre- 
vailing gas stream pressure. Figure 18 illustrated that the 
difference between those two temperatures was generally 
significant, with the maximum difference for the subject 
investigation being equal to approximately 350F. A similar 
comparison is made in this section between the experimental 
values for the maximum liquid temperature and those that 
are predicted by the theory outlined in the foregoing sub- 
section. 

The primary equation ( 3 - 87 b ) was employed, together 
with the thermodynamic data presented in Appendix D of the 
present report, to calculate the values for m for the gas 
stream pressure and temperature corresponding to the experi- 
mental tests conducted. Figure 35 illustrates the results 
of those computations. Presented in the figure is a plot 
of the measured value for the liquid temperature versus the 
predicted value for the temperature m . The cross-hatched 
areas in Fig. 35 represent the range of the experimental 
values for the maximum liquid temperature that were re- 
alized for each liquid at a particular value for the gas 
stream temperature and pressure**. The agreement between 
the experimental “nalytical results shown in Fig. 35 is 
generally quite satlsfactor*. The only substantial 

* water, methanol, butanol, and RP-1 

♦♦Several hundred datum points are represented by the 
cross-hatched areas in Fig. 35. 


measured (F) 



Tj m i present theory (F) 


Fig. 35. Comparison of Present Theory With Data 
for Maximum Liquid Temperature 
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difference between the measured and the predicted values for 
Tj m occurs for some of the experimental data that were ob- 
tained with RP-1 as the film coolant. However, part of that 
difference may be due to the fact that RP-1 was the liquid 
for which the maximum temperature rise occurred after it was 
introduced onto the test plate. Therefore, the liquid film 
may not have always achieved its maximum temperature by the 
time the film reached the point of withdrawal. However, even 
if that were not the case, the results of Fig. 35 suggest 
that the simple theory presented herein for the prediction 
of the maximum film temperature is sufficiently accurate for 
the purposes of most engi neeri ng- type calculations. 
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4. CORRELATION OF THE DATA OF 
THE PRIMARY REFERENCES 

4.1 Introduction 

The experimental data for mass transfer that are re- 
ported in the subject report are limited in two fundamental 
respects: (a) the liquid-film cooled length was held con- 

stant at 10 inches, and (b) the gas stream temperature was 
limited to oOOF and the gas stream pressure was limited to 
150 psia so that the experimental problems would be mini- 
mized. The purpose of the present section is to analyze 
those data that have been reported in the literature for the 
problem of liquid-film cooling that are not so restricted. 
The analysis of those data will give some insight into the 
general usefulness of the correlation for mass transfer that 
is presented in the subject report. 

The experimental data that are reported by the primary 
references (19) (20) and (24) are analyzed. The data 
reported by Kinney, et al. (19) and those reported by 
Emmons (20) were obtained under essentially constant gas 
» l< earn velocity conditions. The basic experimental config- 
uration for each of those two Investigations was the flow 
of a hot gas through a horizontally mounted tube with the 
liquid film Introduced around the circumference of the tube. 
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The experimental data reported by Hermann, et al . (24) are 
for the case of a strongly accelerating gas flow. They 
investigated the problem of liquid-film cooling an axi- 
symmetric convergent-divergent nozzle. For that investi- 
gation, the gas stream was accelerated from a Mach number 
of approximately 0.07 at the plane of liquid injection to 
a Mach number of 1.92 at the exit of the nozzle. 

Table 3 presents a comparison of the nominal flow pa- 
rameters for the experimental investigation of Kinney, et al. 
and of Emmons with those for that portion of the present 
investigation where the gas stream velocity was essentially 
constant over the length of the wetted test plate. Table 3 
shows that the experimental parameters for the subject 
references represent a substantial extension over those for 
the experimental program reported herein. For example, the 
maximum gas stream temperature, T^, investigated in both 
references is significantly greater than the maximum tem- 
perature of 600F for the present investigation. Moreover, 
the gas stream pressure that was employed by Emmons is more 
than three times that employed in the subject investigation; 
the maximum gas stream velocity that was investigated by 
Kinney is nearly four times the maximum that Is reported 
herein; and the maximum film cooled length that Is reported 
by Kinney Is more than three times longer than the 10-Inch 
film cooled length that was employed in the present investi- 
gation. 


Comparison of the Nominal Flow Parameters for the Primary References 
to Those for the Present Investigation--Constant Velocity Gas Stream 
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Table 4 presents a similar comparison of the nominal 
flow parameters for the experimental investigation of 
Hermann, et al. with those for that portion of the subject 
investigation where the gas stream was accelerated over the 
length of the film cooled test section. Again, the experi- 
mental parameters for the referenced investigation represent 
a considerable extension over those for the experimental 
program reported herein, particularly with reference to the 

Mach number, M_ , the gas stream velocity, u„, and the static- 
a g 

stagnation pressure ratio, p/p°. 

Obviously, any evaluation of the present theory in 
terms of the referenced experimental data is only as con- 
clusive as is the accuracy of those experimental data. In 
that regard, it should be mentioned that the wetted surface 
area (and therefore the net rate of mass transfer from the 
liquid film to the gas stream) was determined in each of the 
subject references by external wall temperature measurements. 
Moreover, each reference employed the assumption that the 
liquid film terminated at that point where the wall tem- 
perature just exceeded the boiling temperature for the 
liquid at the prevailing local static pressure. Section 3.4 
herein showed, however, that that assumption Is never ex- 
actly correct, with the error in the assumption Increasing 
with decreasing gas stream temperature and/or decreasing 
gas stream pressure. Thus, there Is some limitation in 
accuracy Introduced by the experimental technique employed 
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Table 4 

Comparison of the Nominal Flow Parameters of the 
Primary Reference (24) to Those of the Present 
Inve$tigation--Accelerati ng Gas Stream 



Parameter 

Present 

Investigation 

Reference (24) 1 

< 

T° 

T g ’ 

Gas Stream Stagnation 
Temperature, R 

860 

| 

1460-1940 1 

p° » 

Gas Stream Stagnation 
Pressure, psia 

50-150 

69-82 J 

p/p°. 

Stati c-Stagnation 
Pressure Ratio 

0.84-0.98 

0.14-0.99 | 

T/T° , 

Static-Stagnation 
Temperature Ratio 

0.95-0.99 

0.58-0.99 J 

M a * 

Gas Stream Mach Number 

0.15-0.51 

0.07-1.92 - 

u g * 

Gas Stream Velocity, fps 

200-660 

150-2750 * 

Liquid Film Coolants 

Water 

Methanol 

Water J 


] 

3 


1 

I 

I 
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In the aforementioned investigations. 

Furthermore, Kinney, et al. (19) comment that the ex- 
perimental data for mass transfer that they reported could 
be reproduced in only a few cases to within 5 percent, and 
more generally to within only 10 or 15 percent. They also 
state that the accuracy of the experimental data for the 
short film cooled 1 engths--! ess than, say, 12 inches--was 
poorer than that for the longer film cooled lengths that 
were investigated. The problem that they experienced for 
the shorter wetted lengths was that the liquid film seldom 
terminated uniformly around the circumference of the wetted 
tube. Thus, it was difficult to interpret the external wall 
temperature measurements and to accurately determine the 
actual wetted surface area. Of the three references con- 
sidered in the present section, the research program re- 
ported by Kinney, et al. is considered the most compre- 
hensive and was conducted under the most well controlled 
experimental conditions. Thus the confidence In the accu- 
racy of the experimental data for Refs. (20) and (24) is 
somewhat less than that for the experimental data of Ref. 
(19) because of the short film cooled lengths that were 
employed In those two Investigations. 

One additional problem occurred in the evaluation of 
the theory of the present report from the data that were 
presented In the subject references. Namely, how to 
properly determine the value for the surface tension o that 
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enters into the entrainment group X fi . In each of the three 
subject references, the liquid coolant was injected at 
essentially the ambient temperature, so that a substantial 
change occurred in the value for the temperature of the 
liquid (and therefore the value for o) over tha length of 
the wetted test section. This problem was particularly 
acute for the investigation due to Emmons because the final 
equilibrium film temperature at the end of the liquid s.lm 
was generally very high. For example, for the film cooling 
experiments reported by Emmons with water as the liquid 
phase, the temperature for the water increased from approxi- 
mately 70F at the point of liquid injection to something in 
excess of 400F at the terminus of the liquid film. The 
corresponding decrease in the surface tension for the water 
coolant was roughly 50 percent. Obviously, the change In 
the temperature of the liquid coolant over the length of the 
wetted test section creates some uncertainty about how t*e 
parameter X fi should be evaluated. It is also obvious, 
however, that such a variation in the liquid film tempera- 
ture would characterize the general application of liquid- 
film cooling. Because of the genera* Importance of the 
problem of determining a proper reference temperature for 
tl.c calculation of the surface tension o, this problem is 
discussed In greater detail In Section 4.2 that follows. 

Sections 4.3, 4.4, and 4.5 consider the analysis of the 
data reported In Refs. (19), (20), and (24), respectively. 
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The basic assumptions that were employed to analyze those 
experimental data are listed in each section, and a compari- 
son is made in each case between the reported rates of mass 
transfer and those that are predicted from the present 
theory. 


4.2 An Approximate Method for 
Evaluating the Surface Tension 

An approximate method for the evaluation of a reference 

liquid temperature, which will be denoted by T, , that can 

■ » ■ 

be employed to estimate the value for the surface tension in 
the group is presented in this section. The method 
strictly applies to only the case where the wetted surface 
Is adiabatic*, but it could be readily extended to the more 
general case where the film cooled region is also cooled 
externally. 

For the general case of liquid-film cooling, the bulk 
temperature for the liquid increases from a minimum value 
at the point of injection, T^ to a maximum value at (or 
near) the end of the liquid film. That maximum film tem- 
perature, denoted by T^ m , can be estimated from the theory 
outlined previously In Subsection 3.4.1. The basic problem, 
therefore, is to determine how the minimum and the maximum 
film temperatures should be "weighted" so as to determine 
the reference temperature Tj f . 

First it is convenient to Introduce a new parameter 4 


♦This case Includes the subject references (19), (20), (24). 
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which is defined as the average value for the energy balance 
parameter <J> over the length of the liquid film. The param- 
eter <fT can be reasonably approximated by the following ex- 
pression. Thus* 


<f> =H +h 1 -h, . 
Y v ,m 1 ,m 1,1 


(4-1) 


i, = the sensible enthalpy for the liqui 
* m at the temperature T^ m 


d evaluated 


The rate at which the bulk temperature for the liquid 

film approaches the maximum film temperature, T, m , once 

that it is injected onto the test plate is a function of 

how much of the total energy transfer to the surface of the 

liquid film results in an increase in the sensible enthalpy 

of the liquid and how much results in the evaporation of 

liquid. That is, it is dependent on the relative magnitudes 

of the terms H„ m and (h, m - h, ,) on the right of Eq. 
v ,m I ,m 1,1 

(4-1). Thus, a weighting parameter, denoted by "w", can be 

*It can be shown that in most cases, the group of terms 
(H y + h-j - h.| y) t where H y and h^ are evaluated at any 

temperature between Tj ^ and Tj m is essentially a con- 
stant over that range of liquid temperatures \ct). Thus 
the evaluation of H and h, at T, Is mostly a matter of 
convenience. v * m 
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6 


defined as follows: 

w ■ H Vjln / ♦ (4-2) 

To illustrate: if the value for w is large (i.e., near 

unity) so that most of the energy transfer to the liquid 
film results in evaporation of liquid, then the reference 
temperature T-j r should be chosen as nearly equal to the 
maximum film temperature m . However, if w is very small, 
so that most of the energy transfer to the liquid film re- 
sults in an increase in the sensible enthalpy of the liquid, 
then the average temperature for the liquid over the length 
of the wetted test surface will be much less than T, m , and 
the reference temperature Tj r should be weighted more in 
favor of the minimum liquid temperature Tj 1 . Therefore, in 
view of the lack of any definitive data for the problem, it 
is suggested herein that the reference temperature T^ m be 
approximated by the following expression. Thus 

T l,r 0-) T ,.l < 4 - 3 > l 

i 

! 

Equation (4-3) was employed to estimate the surface 1 

tension for the analysis of the data of the primary refer- 
ences (19) and (20). The values that were computed tor T, _ 
are given In the appropriate subsection. 

i 

■■ . ..I . ■■ i m m m i uww--- r '"*“ 
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4.3 The Data Due to Kinney, Abramson, and Sloop* 

4.3.1 Introductory Remarks 

Kinney, et al. (19) present experimental data for the 
liquid-film cooling of a 2-inch and two 4-inch inside diam- 
eter ducts. The 2-inch diameter duct and one of the 4-inch 
diameter ducts were seamless tubes having a honed inner sur- 
face. Tiiece two ducts are termed herein the smooth-surface 
ducts. The second 4-inch diameter duct was a 4-inch inside 
diameter rolled tube with a longitudinal weld. Although the 
inside surface of the tube was machined after welding, some 
surface roughness and waviness remained. This tube is 
termed herein the rough-surface tube. 

Subsection 4.3.2, which follows, describes the basic 
experimental apparatus for Ref. (19). The primary assump- 
tions that were employed in order to obtain a comparison of 
the subject data with the theory that is developed in the 
present report are outlined in Subsection 4.3.3. Subsection 
4.3.4 then presents the comparison between the reported 
experimental data and the calculations from the present 
theory. 


4.3.2 The Basic Apparatus 

Figure 36 illustrates schematically the basic experi- 
mental apparatus employed In the investigation reported in 
Ref. (19). The basic apparatus consisted of a horizontally 


♦Refer to Table 3 for the nominal flow parameters for these 
data. 
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mounted Inconel approach section that was 40 inches long and 
a film cooled test section that was 48 inches long. The 
walls of the experimental duct work were essentially adia- 
batic. Mixing baffles and a calming section were installed 
upstream of the approach section in order to obtain a uni- 
form gas stream having a relatively low level of turbulence. 

A gasoline-air combustor was employed to produce the hot gas 
stream. 

The liquid injector that was employed was fabricated 
from a porous wire cloth (36 mesh) that was mounted flush 
with the internal flow passage. The porous injector element 
was two inches long in the direction of the gas flow. The 
data reported were for a water film. 

4.3.3 The Assumptions 

The basic assumptions that were employed to analyze the 
experimental data for the subject reference are as follows: 

(a) The gas flow was treated as boundary-layer type flow 
over a flat plate. The velocity boundary layer was assumed 
to be fully turbulent and to develop from the leading edge 
of the 40-inch long approach section. 

(b) The approach section was assumed to be adiabatic so 
that a thermal entrance length x^ = 40 Inches existed. 

(c) The thermodynamic properties for the gas stream were 
assumed to be the same as those for dry air at the prevailing 
gas stream temperature. The value for the gas stream 
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Prandtl number, Pr ; was taken to be 0.7 for all computations. 

(d) Radiation heat transfer was assumed negligible. 

(e) The reference temperature T, , for the evaluation of 
the surface tension o in the entrainment group X e , was esti- 
mated from the procedure outlined in Section 4.2 to be 165F 
for the lowest gas stream temperature that was investigated 
(800F), and 210F for the highest gas stream temperature that 
was investigated (1600F). The values for the reference 
temperature for intermediate values of the gas stream tem- 
perature were determined by a linear interpolation between 
those limiting values. 

(f) The value for the energy balance parameter <|> was cal- 
culated from Eq. (4-1). 

4.3.4 Comparison of the Present Theory 
With the Experimental Data 

The primary equation of the present theory is repeated 
here for convenience. Thus 

^ * (1.0 + exp[5.0 x 10' 6 (x e - 1000) (l/l 0 ) b ] 

s X r 

(3-64) 

Equation (3-64) gives the desired relationship between the 
rate at which liquid is injected, n^, and the resultant film 
cooled length, 1. The simple rate of mass transfer, m‘, was 
computed from the procedure outlined In Subsection 3.1.2, 
together with the aforementioned assumptions. Figure 29 of 
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Subsection 3.2.2 predicts a value of approximately 0.6 for 
the parameter b in Eq. (3-64) for the subject flow situation. 
The present subsection will investigate the validity of that 
prediction (and, therefore, the validity of the model that 
was employed in the analysis of Subsection 3.2.1). 

Figure 37 presents a comparison of the experimental 
data for mass transfer in terms of the film cooled length 
versus the injected liquid flowrate for the 4-inch diameter, 
smooth-surface tube with the results that were computed for 
both (a) the simple theory for mass transfer (Fig. 37(a)), 
and (b) the present theory for mass transfer (Fig. 37(b)). 
Figure 38 presents a similar comparison for the 2-inch 
diameter, smooth-surface tube. The solid curves presented 
in these two figures represent the data that are presented 
in Fig. 6 of Ref. (19). Approximately 25 datum points are 
represented by each of the solid curves. The values for the 
gas stream temperature and the gas stream Reynolds number 
(based on the tube diameter D) are indicated in the figures 
for each set of experimental data. 

Figures 37(a) and 38(a) again illustrate how inadequate 
the simple theory is for the problem of liquid-film cooling. 
For example: Figure 37(a) shows that at a gas stream 

Reynolds number of 4.7 x 10® and a gas stream tpmne^ture of 
1600F, the value for required to establish a film cooled 
length of 2.2 ft is more than three times the value for m^ 
that is predicted from the simple theory for mass transfer. 
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The possible catastrophic consequences of designing a liquid- 
film cooled system on the basis of the simple theory are 
evident. 

Figures 37(b) and 38(b) present the resultant computa- 
tions for the present theory of mass transfer. Because the 
values for the gas stream pressure were not indicated explic- 
itly for the data that were reported in Ref. (19), a band is 
shown for the present theory for each set of experimental da- 
ta. For the two-inch diameter duct, p was assumed to vary 
from 30 to 35 psia. An., f the 4-inch diameter duct, p was 
assumed to vary from 25 to 30 psia. The larger values for m. 
at any particular value of 1 correspond to the lower values 
for pressure. (For the same value of gas stream Reynolds 
number, or gas mass velocity G, the value for the momentum 
parameter M Q increases with decreasing pressure. And the 
amount of entrained liquid increases with increasing M Q .) 

A value of 0.5 for the parameter b in Eq. (3-64) was employed 
to compute the results that are shown in Figs. 37(b) and 
38(b) for the smooth-surface ducts. Recall that the simple 
theory presented in Subsection 3.2.2 for the extension of 
the results of the present investigation to an arbitrary 
film cooled length showed that for the classical case of 
turbulent boundary layer flow over a flat plate, a value of 
approximately 0.6 for the parameter b was predicted. That 
value for b resulted in predicted mass transfer rates that 
were generally 10 percent greater than the calculated 
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results that are shown in Figs. 37(b) and 38(b). However, 
as discussed earlier in Section 4.1, that agreement between 
the predicted and the experimental results is probably with- 
in the accuracy of the experimental data. 

Comparison of Figs. 37(b) and 38(b) shows that the 
agreement between the experimental results and the results 
calculated from the present theory is somewhat better for 
the 4-inch diameter, smooth-surface duct than it is for the 
2-inch diameter duct. However, at least two reasons can be 
sugg*?sted as to why the experimental data for the 4-1:.ch 
djameter duct should be weighted more heavily in an eval- 
uation of the present theory: (a) the gas flow in the 

4-1nch diameter duct should more closely correspond to the 
case of boundary layer flow over a flat plate*, and more 
Importantly, (b) it was stated In the subject reference that 
the agreement between the data for any two experiments con- 
ducted at the same flow conditions was generally much poorer 
for the experiments with the 2-inch diameter duct than it 
was for the 4-inch diameter ducts. 

Figure 39 presents a similar comparison between the 
calculated results and the experimental results for the 
4-inch diameter, rough-surface tube. The solid curves 
shown in the figure represent the data that are presented 
in Fig. 6(c) of Ref. (19). Approximately 15 datum points 
are represented by each of the solid curves. 


*See assumption (a) in Subsection 4.3.3. 
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Figure 39(a) again shows the poor agreement between the 
experimental results and the results that were computed from 
the simple theory for mass transfer. The mass transfer 
rates that were realized experimentally are shown in Fig. 
39(a) to be as much as 400 percent greater than those pre- 
dicted by the simple theory. 

The computed results shown in Fig. 39(b) for the pres- 
ent theory were calculated from Eq. (3-64) with a value of 
0.6 for the parameter b. The band that is presented for each 
set of experimental data represents a variation in the static 
pressure of 25 to 30 psia. In general, the agreement be- 
tween the experimental data and the results that are calcu- 
lated from the present theory is satisfactory. The greatest 
difference occurs for the experimental data at the larger 
rates of liquid injection, m^. However, even for those data, 
the present theory is in agreement to within the reported 
accuracy of the experimental data. 

The following conclusion is drawn from the foregoing 
comparison of the computed results of the present theory 
with the experimental data of Ref. (19) for both the smooth- 
surface ducts and the rough-surface duct: to perform con- 

servative engineering calculations so as to allow for un- 
known factors such as the physical condition o* *be wetted 
surface, the value of 0.6 for the parameter b in Eq. (3-64) 
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should be employed*. That conclusion is in good agreement 
with the result that was predicted by the analysis of Sub- 
section 3.2.1 and Subsection 3.2.2. 

4.3.5 A Comment on the Presentation of Data for 
Liquid-Film Cooling in Terms of Dimensionless Parameters 

Consider the presentation of the experimental data of 
Ref. (19) in terms of a Nusselt number-Reynol ds number type 
of plot. The effective** average Nusselt number, denoted 
by Nu g , can be calculated from the following expression. 
Thus 


Nu. 


m^yi Re D Pr 


(4-4) 


Figure 12 of Ref. (19) presents plots of Nu g as a function 

of Rep showing some of the experimental data for both the 

smooth-surface tubes and the rough-surface tube. Figures 

40 and 41 of the present report show a comparison between 

those data and the results that are calculated from the 

* Note that this conclusion applies to only that case 
where the gas stream velocity is essentially constant 
over the length of the liquid film. 

**The term "effective 1 ’ is introduced because Nu g includes 

the mass transfer due to entrainment in addition to that 
due to evaporation. 


effective average NusseU number x 10 





effective average Nusselt number 
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present theory*. Figure 40 shows the results for the 
smooth-surface ducts, and Figure 41 presents the results for 
the rough-surface duct. Note that each set of experimental 
data is characterized by a certain value for the rate of 
liquid injection, m^, the values for m^ are listed in the 
figure. The interrupted curve shown in both Fig. 40 and 
Fig. 41 corresponds to the case of simple mass transfer, 
where the various "F" correction factors (see Subsection 
3.1.2) were employed in the evaluation of the average 
Nusselt number. Because these correction factors are de- 
pendent on parameters other than the gas stream Reynolds 
number, the curves for simple-mass transfer are labeled 
"approximate". In evaluating the net mass transfer rate, 
n^, a value of 0.5 for the parameter b in Eq. (3-64) was 
utilized in the calculations for the smooth-surface ducts 
and a value of 0.6 was employed for the rough-surface duct. 

*In the calculation of the Nusselt number in Pef. (19), 
the thermal conductivity for the gas stream, kg, was 

introduced. However, the values for k g that were employed 

in those calculations were generally 10 percent to 15 per- 
cent greater than those listed, for example, in Ref. (40). 
To illustrate. Appendix C of reference (19) gives a value 

for kg = 15.6x10*® (Btu/sec-f t-F ) for a gas stream tem- 
perature T 0 = 2000F. But Ref. (40) gives a value for 
9 

kg = 13.1x10“ (Btu/sec-f t-F) at the same temperature. 

Therefore, to make some correction for this fact, the 
values for Nu e presented in Figs. 40 and 41 are greater 

than those presented in the original reference by a factor 
of 1.10. . 




7T 
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The bands shown in Fig. 40 and Fig. 41 for the theory that 
was developed in the present report correspond to the gas 
stream pressure variation indicated in the figures. Figures 
40 and 41 demonstrate again the poor agreement of the simple 
mass theory and the much more satisfactory agreement of the 
present theory with the experimental results. However, a 
more important point should be noted from the comparisons 
illustrated in the figures. Investigators of the problem 
of liquid-film cooling have commonly presented their experi- 
mental data in the form of such a Nusselt-Reynolds number 
type of plot (19, 20, 24). The experimental data of the 
present investigation, however, together with the data due 
to Kinney, et al. (19), show that this is generally a very 
poor way to present such experimental data because it does 
not account for the considerable importance that the param- 
eters M , a, and m, can have on the net rate of mass trans- 
fer from a thin liquid film to a high velocity gas stream. 
This type of data presentation can be very misleading for 
the practicing engineer who is concerned with the problem 
of designing a liquid-film cooled system, but who is not 
cognizant of the importance of the interfac’al phenomena, 
particularly the entrainment phenomenon. It is hoped that 
the results of the experimental investigation reported here- 
in clearly demonstrate the general importance of the inter- 
facial phenomena in the calculation of the net rate of mass 
transfer from a thin liquid film to a hot gas. stream. 
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4 . 4 The Data Due t o F mm on s * 

4.4.1 Introductory Remarks 

Emmons (20) presents experimental data for liquid-film 
cooling of a cylindrical rocket motor combustion chamber. 
Experimental determinations were made of the effects of 
changes in the gas stream temperature, pressure, and Reynolds 
number upon the required flow rate of a single liquid (water) 
needed to establish a range of :m cooled lengths (1 = 4 to 
8 inches). Experiments were also conducted with other 
coolants**, at one gas stream flow condition, to study the 
influence of the physical property values for the liquid on 
the rate of lass transfer from the liquid film to the gas 
stream. Most of the data presented in Ref. (20) were ob- 
tained for a gas stream temperature of 4100 R and a chamber 
pressure of 500 psia. Only those data are reviewed herein. 

Subsection 4.4.2, which follows, describes the basic 
experimental apparatus that was employed by Emmons. The 
primary assumptions employed in the calculation of the mass 
transfer rates from the theory of the present report are 
discussed in Subsection 4.4.3, and the comparison is shown 
in Subsection 4.4.4 between the measured and the predicted 
results. 

* See table 3 for the nominal flow parameters for these 
data. 

**anhydrous ammonia, ethyl alcohol, and Freon-113 
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4.4.2 The Basic Apparatus 

Figure 42 illustrates schematically the basic experi- 
mental apparatus employed by Emmons in his study of liquid- 
film cooling. The experimental rocket motor burned gaseous 
hydrogen and air as the propellants. The inside diameter of 
the chamber was 3.14 inches. Two convecti vely cooled cylin- 
drical sections were installed upstream of the adiabatic 
film cooled test section. The first convectively cooled 
chamber was 9.42 inches long and the second convectively 
cooled chamber was 1.48 inches long. Emmons utilized the 
1.48-inch long section to calorimetrical ly estimate the heat 
transfer to the dry wall immediately upstream of the film 
cooled test section. An observation will be made in Sub- 
section 4.4.5 regarding the measurements that were obtained 
for that 1.48-inch long section. 

The liquid coolant was injected through 72 V-shaped 
grooves that were milled into the downstream face of the 
film coolant injector. The liquid was introduced onto the 
test surface with a component of velocity in the tangential, 
or circumferential, direction. 

4.4.3 The Assumptions 

The primary assumptions employed to analyze the experi- 
mental data of Ret. (20) are summarized as follows: 

(a) The gas flow was treated as boundary-laysr type flow 
over a flat plate. The boundary layer was assumed to be 



adiabatic filn-cooled 
test section 
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fully turbulent and to develop from the leading edge of the 
9.42-inch long convectively cooled section. 

(b) The temperature for the liquid film at the point of 
injection was much less than the wall temperature immediately 
upstream so that a thermal entrance length x t = 9.42 + 1.48 = 
10.90 inches existed. 

(c) No chemical reactions occurred between the vapor for 
the film coolant and the main stream gas within the boundary 
layer that developed over the liquid film. 

(d) The value for the energy balance parameter <j> was cal- 
culated from Eq. (4-1). The computed values for 4> are shown 
in Table 5 below. 

(e) Radiation heat transfer was neglected. 

(f) The properties of the gas stream were estimated from 
the results presented by Rannie in Ref. (43)*. For a gas 
stream temperature of 4100 R and a chamber pressure of 500 
psia, the pertinent thermodynamic properties that were 
estimated from Ref. (43) are as follows: 

M g - 23.0 lb/lb-mole 
c pg = 0.43 Btu/lb-R 

Y - 1.255 

u g e 4.75 x 10" 5 lb/ft-sec 

♦Emmons (20) docs not list the values that he assumed for 
the thermodynamic properties of the hot gas stream. 
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(g) The reference temperature* for the evaluation of the 
surface tension o was estimated from the procedure that was 
suggested in Se*ction 4.2. Table 5 presents the values for 
the pertinent parameters thav entered into that evaluation, 
together with the computed values for T-j r and o . Also 
listed in Table 5, for completeness, are the values for the 
heat transfer parameter B^ and the interfacial vapor con- 
centration C $ . was calculated from the right-hand side 
of Eq. (3-87b), and C $ was calculated from Eq. (3-26)*. 

Of the aforementioned assumptions. Assumption (3) that 
the heat transfer due to radiation is negligible is probably 
one of the leas* valid. The products of combustion of air 
and gaseous hydrogen contain a substantial quantity of water 
vapor which is a strong radiator, parti cul arly at the tem- 
perature level of 4100 F. Perhaps 25 percent or more of 
the net heat flux to the surface of the liquid film was due 
to radiation heat transfer. An attempt was not made, however, 
to estimate that radiation heat transfer. Significant con- 
clusions can be drawn about the usefulness of the theory 
developed in the present report for the prediction of the 
rate of mass transfer from a liquid film to a hot, high 
velocity gas stream without trying to solve that complex 
problem 

♦Appendix D.2 presents H , p , and o as a function of 

V V i s 

the liquid temperature for the liquids that were investi- 
gated by Emmons. 
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4.4.4 Comparison of the Present Theory 
With the Experimental Data 

The aforementioned assumptions were employed in con- 
junction with the primary equation (3-64) to analyze the 
data of Ref. (20). The results of those computations are 
summarized in the following paragraphs. 

Figure 43 shows the comparison of the computed results 
and the experimental data for the case where the gas stream 
Reynolds number was a variable and water was the film cool- 
ant. The gas stream Reynolds numbers are listed on the 
figure. Part (a) of the figure presents a plot of the mass 
transfer rate that was predicted from the simple theory for 
mass transfer versus the measured mass transfer rate. Part 
(b) of Fig. 43 then shows the comparison between the experi- 
mental data and the present theory. A value of 0.6 for +he 
parameter b in Eq. (3-64) was employed to compute the re- 
sults for the present theory. 

Figure 43(b) shows that the present theory for mass 
transfer predicts a mass transfer rate that is consistently 
less than that realized experimentally. However, Fig. 43 
also shows that the present theory is a substantial improve- 
ment over the simple-theory for mass transfer. Moreover, 
if the radiation heat transfer was accounted for, the agree- 
ment between the present theory and the experimental data 
would be improved. Indeed, if one also accounts for the 
likelihood of considerable experimental error in the 





183 


f 




to 

to 

TO 


** 

>» 

s: 

+- 

• 

U 


u 

•1 

O 

s> to 

«* 

a* 

O i- 

JS 

X! 

*4— <D 


f— 

JQ 

V 


TO B 

t 

-M 

4-> 3 

3 

¥» 

C 

to h: 

<• 

0) 

o 

£ 

to 

to 

* 

TO 

-C "O 

C* 

i- 

+* r- 


Cl. 

*f~ o 



3: c: 






>> to 


X> 

i~ CC 


"* — 1 

O 


<D E 
x: A3 
h- <U 


s- 

'M +> 
C CO 
TO 

to to 
TO TO 
i- CD 
a. 

10 

•a 3 

c: a 

tD *r* 
*- 
>> TO 
C- :> 
O 
TO 

SZ * 
f— CO 

c 
TO o 
r— E 
O, E 
B uj 


CO o 



>> 

+-> 


s- 


£ 

o 

o a> 

£ 

TO 

3 

* 

M 

x: 

C Q 

£ 

H- 

O 



to u 


0) 

•r~ TO 

? 

3 

r— • 

U 4- 

Cl 

TO V) 

%* 

« 

B 

D. C 

8f 

*r- 

E to 

• 

<S> 

o u 

r 


O H* 


<0 

• 



CO 







Ojoa«/<||) XUoim «|dui|« • *« 



184 


reported data*, it could probably be argued that the agree- 
ment between the present theory and the reported experi- 
mental data is quite satisfactory. 

Figure 44 presents a similar comparison for the data 
that were obtained for the film coolants water, ammonia, 
ethanol, and Freon-113 at a gas stream Reynolds number of 
0.68 x 10 5 . Again, part (a) of the figure compares the 
experimental data for mass transfer with the mass transfer 
that was predicted from the simple theory; and part (b) 
of Fig. 44 shows the comparison between the theory that was 
developed herein and the reported data. Figure 44 demon- 
strates further that the proposed correlation for mass trans- 
fer is an improvement over the simple theory. If the 
neglect of the radiation heat transfer and the probable 
experimental error is taken into consideration, the agree- 
ment between the present theory and the experimental data 
shown in Fig. 44 is quite satisfactory. 
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4.4.5 An Additional Comment of the 
Research Due to Emmons 

It was indicated in Subsection 4.4.2 that the 1.48-inch 

long convectively cooled section** (just upstream of the 

film cooled test section) was employed to determine the rate 

of heat transfer to a ury wall as a function of the pertinent 

* See the discussion of Section 4.1. 

**Refer to Fig. 42. 







Fig. 44. Comparison of Simple Theory and Present Theory with Dat 
for Mass Transfer Due to Emmons. Various Film Coolants 








gas stream flow conditions, Emmons compared those dry-wall 
measurements to the measurements ho obtained for the liquid- 
film cooled wall region and concluded that the rote of heat 
transfer to the liquid film was loss than the rote of heat 
transfer to the dry wall for the same gas stream conditions. 
In essence, he concluded that the analogy between simple- 
mass- transfer cooling and liquid-film cooling was valid. 
However, the experimental data reported herein, together 
with those reported by Kinney, et al. (19), show rather con- 
clusively that the analogy between the simple-mass transfer 
process and the process of liquid-film cooling is generally 
very poor. In fact, the computed results shown In fig. 43 
and Fig. 44 suggest that the experimental data reported by 
Emmons Ho not support the analogy. This apparent discrep 
oncy Is discussed In the following paragraph. 

Emmons found that the dry-wall heat transfer data could 
bo correlated by a Mussel t~ Rcynol ds number type of corre- 
lation that was suggested by Humble, et al. (44) for fully- 
developed turbulent tube flow. That particular corre- 
lation, however, specifies that the properties for the gas 
that enter Into both the Nusselt number and the Reynolds 
number be evaluated at the wall temperature. The evaluation 
- r the properties at the wal l temperature is goimron pro- 
cedure employed to correlate heat transfer data for tbo case 
where heat Is transferred from the wall to the gas stream 
(36) ( 1 . e . , the wall is at a higher temperature than the gas 
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cjas stream flow conditions# Emmons compared those dry-wall 
measurements to the measurements he obtained for the liquid* 
film cooled wall region and concluded that the rate of heat 
transfer to the liquid film was less than the rate of heat 
transfer to the dry wall for the same gas stream conditions, 

o 

In essence, he concluded that the analogy between simple* 
mass* transfer cooling and 1 1 qul d-f 1 lm leal Ing was valid. 
However, the experimental data reported herein, together 
with those reported by Kinney, et al. (10), show rather con- 
clusively that the analogy between the simple-mass transfer 
process and the process of liquid-film cooling Is generally 
very poor. In fact, the computed results shown In fig, 43 
and fig, 44 suggest that the experimental data reported by 
Emmons do not support the analogy, this apparent discrep- 
ancy is discussed in the following paragraph , 

Emmons found that , the dry-wall heat transfer data could 
be correlated by a Nusaolt- Reynolds number type of corre- 
1 a 1 1 o n t h a t w a s § u g g e s ted by Humble, et al, (44) for f u 1 1 y - 
developed turbulent tube flow, That particular corre- 
lation, however, specifies that the pfoperties for the gas 
that enter Into both th& Nussel t number and the Reynolds 
number^ be evaluated at the wall temperature , The oval uat|on 
* f the p roperti es at the wa 1 1 temperature i s « v*H*imon pro- 
cedure employed to correlate heat transfer data for the case 
where heat Is transferred from the wall to the gas stream 
(36) ( i . a , ; the wall is at a higher temperature than the gas 
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stream). That was the case that was investigated by Humble, 
et al . (44). However, for the case where the heat transfer 
is from the gas stream to the wall, as is the case for the 
data reported by Li ors, the standard procedure is to eval- 

uate the physical properties for the gas at some inter- 
mediate temperature between that for the gas stream and that 
for the wall. The choice of the reference temperature for 
the evaluation of the properties for the gas can be very 
important. To illustrate, Kays (36) suggests that the 
Nusselt number for the case of turbulent heat transfer from 
the wall to the gas stream is greater than the Nusselt 
number for the case of turbulent heat transfer from the 
gas stream to the wall, all other conditions being equal, 
by the factor (i /T $ ) ' . For the dry-wall experiments of 

Emmons, a characteristic value of (T /T g ) would be 4 with 
the resulting value for that factor being 1.41. It would 
appear, therefore, that the heat transfer rates that were 
realized for the dry-wall test section were approximately 
40 percent greater th*n those that would be predicted by a 
standard correlation for fully developed, turbulent gas 
flow with the direction of heat transfer being from the gas 
to the wall of the tube. Part of that excess heat transfer 
could probably be attributed to radiation heat transfer. 

More significantly, however. It can be argued that part 
(or perhaps' most) of that excess heat transfer was probably 
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due to the fact that the assumption of fully-developed tube 
flow was a poor approximation to the flow over the 1.48-inch 
convectively cooled test section, particularly as regards 
the thermal flow field. Reference to the schematic illus- 
tration of the experimental apparatus presented in Fig. 42 
shows that for the 9.42-inch long convectively cooled sec- 
tion, the external coolant (water) was introduced at the up- 
stream end. It follows, therefore, that because of the heat 
transfer from the combustion chamber to the external coolant, 
both the temperature for the wall and for the coolant would 
increase in the flow direction until the coolant was with- 
drawn at a point just upstream of the 1.48-inch long con- 
vectively cooled test section. Assuming that the temperature 
of the external coolant was the same at the inlet to the 
9.42-inch long section and at the inlet to the 1.48-inch 
long section, it then can be argued that there existed a 
decrease in wall temperature (that could have been substan- 
tial) in going from the downstream end of the 9.42-inch 
section to the 1.48-inch section. Such a change in wall 
temperature would cause the rate of heat transfer to the 
1.48-inch section to be significantly higher, and it makes 
the assumption of a fully-developed thermal flow field much 
poorer, than If the wall temperature was everywhere iso- 
thermal. Emmons did not attempt to determine what the 
pertinent thermal boundary conditions were. Thus the valid- 
ity of the above argument cannot be investigated further. 
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The argument is presented here solely as a possible expla- 
nation for the apparent high rates of heat transfer that 
were measured for the dry-wall test section. It is unfor- 
tunate, however, that Emmons arrived at such an important 
conclusion that the liquid-film cooling process and the 
process of simple-mass transfer were analogous without hav- 
ing considered the interpretation of the dry-wall heat trans- 
fer data in greater detail. Also it will be shown in Appendix 
A. 4 of the present report that the major conclusion that 
Emmons drew with reference to his experimental data for the 
gas-vapor cooled region immediately downstream of the ter- 
minus of the liquid film is open to some question because of 
the possible incomplete interpretation of the experimental 
data for the dry-wall heat transfer experiments. 

4 . 5 The Data Due to Hermann, Leitinger, and Melnik* 

4.5.1 Introductory Remarks 

Hermann, et al. (24) experimentally investigated the 
feasibility of liquid-film cooling the nozzle of a hyper- 
sonic wind tunnel. For the preliminary experiments reported 
in the subject reference, the film cooling of a Mach number 
1.92 axi-symmetric Laval nozzle was considered. They were 
not primarily Interested in obtaining quantitative data on 
the rate of mass transfer from a liquid film to an acceler- 
ating gas stream. In fact, only three datum points are 

♦Refer to Table 4 for the nominal flow parameters for 
these data. 
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presented that can be employed to evaluate the present 
theory. Moreover, the accuracy of these data points is 
questionable, primarily due to the fact that the instrumen- 
tation was limited. Specifically, only 15 thermocouples 
that were welded to the external surface of the film cooled 
wall were utilized to determine the average film cooled 
length. The thermocouples were distributed both axially 
along and circumferentially around the nozzle. Only 9 of 
the 15 thermocouples were installed downstream of the throat 
of the nozzle where the liquid film usually terminated. 

Thus, due to the non-uniform way in which the liquid film 
usually terminates, it would be impossible, in general, to 
accurately determine the actual wetted area (or average film 
cooled length) from only 9 thermocouple leadings. However, 
in spite of the foregoing limitations, it is still of some 
Interest to compare the reported data for mass transfer with 
the mass transfer that is predicted by the present theory 
for the case of an accelerated flow. 

Subsection 4.5.2 presents a discussion of the basic 
apparatus, and Subsection 4.5.3 outlines the primary as- 
sumptions that were employed to analyze the experimental 
data. Subsection 4.5.4 presents a comparison between the 
predicted and the experimental results. 

To obtain some Idea of the possible film coolant re- 
quirements for the cooling of the same axl -symmetric nozzle 
If It were Installed on a high-energy propulsion device. 
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Subsection 4.5.5 presents an illustrative calculation based 
on the present theory for the case where the stagnation 
temperature T° = 4000 R and the stagnation pressure p° = 

500 psia. 


4.5.2 The Basic Apparatus 

Figure 45 illustrates schematically the experimental 
apparatus of Ref. (24). The figure presents a Mach number 
1.92, axi -symmetri c Laval nozzle having a 3-inch long sub- 
sonic section and a 5-inch long supersonic section. The 
diameter at the inlet to the nozzle was 8 inches; the di- 
ameter at the throat of the nozzle was approximately 2.7 
inches; and the diameter at the exit of the nozzle was 
approximately 3.5 inches, the walls of the nozzle were 
essentially adiabatic. 

A u-47 jet engine combustor was employed to produce the 
hot gas stream. The maximum stagnation temperature that was 
investigated was 1940 R. 

The film coolant (water) was introduced at the inlet 
to the subsonic portion of the nozzle. The coolant was in- 
jected onto the walls of the nozzle through 48 holes that 
were evenly spaced around the circumference of the nozzle. 
For most of the Investigation, the diameter of each hole was 
0.011 inches. 
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4.5.3 The Assumptions 

The primary assumptions employed to analyze the experi- 
mental data of Ref. (24) are summarized as follows. 

(a) The theory outlined in Subsection 3.3.1 for predicting 
the rate of simple mass transfer from a liquid film to an 
accelerating gas stream was employed to determine the simple 
mass transfer rates for the subject case. The boundary layer 
was assumed to be fully turbulent and to develop from the 
leading edge of the 3-inch long subsonic portion of the 
nozzle. In addition, the adiabatic-wall recovery tempera- 
ture for the gas stream was assumed to be equal to the stag- 
nation temperature T°. 

9 

(b) Except for the thin boundary layer region next to the 
nozzle wall, the gas flow was assumed to be isentropic with 
a specific heat ratio of 1.33. With that assumption, the 
values for p, T g , and u^ were computed as a function of the 
position along the axis of the nozzle. No corrections were 
made in the calculations for displacement effects. 

(c) The thermodynamic properties for the gas stream were 
assumed to the same as those for dry air at the same 
temperature. 

(d) The average temperature for the liquid film throughout 
the nozzle was estimated from Fig. 28 of &»*. (24) to be 
220F for the case of the three datum points that are ana- 
lyzed herein. 

(e) A value of 1000 Btu/lb for the energy balance parameter 
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<j> was estimated from Eq. (4-1) with the substitution of the 
pertinent thermodynamic properties for water. 

(f) Radiation heat transfer was neglected. 

The foregoing assumptions were employed, together with 
the solution technique that was outlined in Subsection 3.3.2 
for the present theory, to analyze the subject experimental 
data. The value for the parameter d in Eq. (3-74) that was 
employed in the calculations was 1.0. That value for d was 
found to correlate the data for the present investigation 
for the case where the gas stream was accelerated over the 
wetted wall region. 

4.5.4 Comparison of the Present Theory 
With the Experimental Data 

Reference (24) presents, for three different gas stream 
conditions, the approximate film coolant flow rate required 
to wet the first 4 to 5 inches of the nozzle. The compari- 
son between the experimental results and the results pre- 
dicted from both the simple theory for mass transfer and 
from the present theory (assuming d = 1.0) is as follows: 
Point 1 : T° » 1460R, p° = 69 psia. 

experimental: m A = 0.037 lb/sec 

simple theory: m ^ = 0.017 to 0.021 lb/sec 

present theory: * 0.023 to 0.028 Ib/sec 
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Point 2 : T° = 1470R, 

p° . 

82 p s i a . 

experimental : 

m i = 

0.037 Ib/sec 

simple theory: 

m i = 

0.019 to 0.23 lb/sec 

present theory: 

m i = 

0.028 to 0.033 Ib/sec 

Point 3 : 1 ° = 1940R, 

. y 

p° - 

82 psia. 


experimental: m 1 = 0.097 Ib/sec 

simple theory: m^ = 0.029 to 0.036 lb/sec 

present theory: m^ = 0.043 to 0.053 lb/sec 

The foregoing results show that the agreement between 
the experimental and the predicted results is satisfactory 
for Point 2, is less satisfactory for Point 1, and is rather 
poor for Point 3. However, for each point, the result pre- 
dicted by the present theory is better than that predicted 
by the simple theory for mass transfer. 

It is unfortunate that additional experimental data 
were not presented in the subject reference in such a form 
that the present analysis could be further evaluated. More- 
over, it is unfortunate that the accuracy of the three datum 
points that are reported is uncertain. One question that 
remains to be answered is to what degree (if any) compress- 
ibility effects for Mach numbers greater than, say, 0.5 
alter the form of the correlation proposed herein for the 
rate of mass transfer from a thin liquid film to an acceler- 
ating, hot gas stream. 
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4.5.5 An Illustrative Example for 
T° = 4000 R and p° = 500 psia 

To illustrate further the present theory, an example is 
presented in this section to show the film coolant require- 
ments (predicted) for the cooling of the above axi -symmetri c 
nozzle for the case where the gas flow has a stagnation tem- 
perature of 4000 R and a stagnation pressure of 500 psia. 
These stagnation conditions are characteristic of the flow 
conditions that might be encountered in the application of 
film cooling to thermally protect an exhaust nozzle of a 
high-energy propulsion device. This example serves to indi- 
cate the type of deviation from the simple theory for mass 
transfer that can be expected in the liquid-film cooling of 
such an exhaust nozzle. 

The assumptions employed in the calculations are the 
same as those listed in Subsection 4.5.3. Again, the value 
of 1.0 for the parameter d in Eq. (3-74) was employed in the 
computati ons . 

Figure 46 presents a plot of the net rate of coolant 
injection as a function of the film cooled length of the 
exhaust nozzle (measured along the axis of the nozzle). The 
pertinent parameters that were employed in the calculations 
are summarized in the -figure. The interrupted curve shown 
in the figure represents the result that Is predicted from 
the simple theory for mass transfer. The two solid curves 
represent the result that is predicted from the present 
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theory for two different values of surface tension a. The 
larger value for c (50 dynes/cm) represents roughly the case 
where water is employed as the film coolant. Most of the 
other film coolants that might be of interest (e.g., alco- 
hols, hydrocarbon fuels) are approximately represented by 
the smaller value for a (25 dynes/cm). Figure 46 shows that 
the present theory predicts, for the case where a = 50 
dynes/cm, a coolant flow rate to cool the subsonic portion 
of the nozzle that is approximately 50 percent greater than 
that predicted from the simple theory. And, for the case 
where o = 25 dynes/cm, the present theory predicts an in- 
crease in the coolant requirements to cool the subsonic 
portion of the nozzle of 100 percent. More interestingly, 
however. Fig. 46 shows that to cool both the subsonic and 
the supersonic portions of the nozzle, by injecting the 
coolant at only the leading edge of the subsonic portion, 
the coolant requirements predicted from the present theory 
become substantial. For the case where a =‘50 dynes/cm, 
about 6 times more coolant is required to cool the entire 
nozzle than is predicted by the simple theory. For the 
case where o = 25 dynes/cm, the amount of coolant required 
to cool the entire nozzle is roughly 30 to 40 times that 
predicted by the simr* 1 o theory. Although thee'' p^'dicted 
coolant flow rates are considerably greater than those pre- 
dicted by the simple theory, the three datum points dis- 
cussed in Subsection 4.5.4 suggest that the estimates from. 
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the present theory may even be low. 

It is worthy of note that almost all of the difference 
between the results predicted by the simple theory and those 
predicted by the present theory is due to the entrainment 
phenomenon. That is, except for the very leading edge of 
the subsonic portion of the nozzle where the gas stream 
velocity is relatively low, the roughness parameter r is 
small, less than, say, 0.10. 

It also is interesting to note the local rates of mass 
transfer from the liquid film to the gas stream that are 
predicted by the present theory for the subject illustrative 
case. Figure 47 presents a plot of the predicted local rate 
of mass transfer from the liquid film into the gas stream as 
a function of the axial location throughout the nozzle. The 
pertinent parameters employed in the calculations are listed 
in the figure. Only the one value for surface tension of 
50 dynes/cm was employed in the computations. The result 
that is predicted by the simple theory is indicated in 
Fig. 47. In addition, three curves are presented in the 
figure to illustrate the result that is predicted by the 
present theory for three values of the net rate of coolant 
injection, m^. The coolant flow rates shown are 0.5, 1.2, 
and 3.0 lb/sec. From Fig. 46, those coolant flow rates 
correspond to predicted film cooled lengths of approximately 
4, 6, and 8 inches, respectively. Figure 47 Illustrates 
that for all cases, the maximum value for the local rate of 
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mass transfer occurs at (or very near) the throat of the 
nozzle. The figure also shows that the local rate of mass 
transfer predicted by the present theory can be substantially 
greater than that predicted by the simple theory. Of course, 
if the liquid coolant has a surface tension less than the 
50 dynes/cm that is represented by the calculated results 
in Fig. 47, then the local rates of mass transfer would be 
substantially greater than those shown. 

In summary, the illustrative example presented in this 
subsection suggests that the simple theory for mass trans- 
fer can lead to a catastrophic design of a liquid-film 
cooled exhaust nozzle. In most cases that are of practical 
Importance, the entrainment phenomenon would completely 
dominate the mass transfer process. 

It also follows from the results of the foregoing cal- 
culations that it is not practical, in general, to Introduce 
the film coolant at only one location. The net rate at 
which liquid Is entrained into the gas stream can be re- 
duced substantially if the coolant is Introduced at several 
axial locations. It would seem that this particular problem 
needs further experimental Investigation. 
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5. CONCLUSIONS 

The following principal conclusions are drawn from the 
results of the subject investigation. 

1. The scale of the small disturbance waves on the surface 

of the liquid film is principally characterized by the momentum 

parameter M Q for the gas stream, with the scale decreasing 

with increasing values of M . 

o 

2. The surface tension apparently affects the film surface 
characteristics because with water as the liquid phase 
(relatively high values of o), a second type of disturbance 
wave was frequently observed (two dimensional, long wave- 
length) that was never observed with either methanol, bu- 
tanol, or RP-1 as the liquid phase (relatively low values 
of surface tension). 

3. The basic Interfaclal structure Is not altered apprecia- 
bly by substantial variations In the liquid viscosity. Also 
It Is not noticeably Influenced by a significant variation 
In the rate of liquid flow In the liquid film so long as the 
gas stream velocity Is greater than about 50 fps. 

4. The rate of mass transfer from a thin liquid film Into 

a high velocity, hot gas stream can be substantially greater 
(typically 5 to 10 times greater) than that predicted by the 



simple theory for mass transfer, with the phenomenon of 
entrainment normally accounting for the bulk of the mass 
transfer. 

5. The mass transfer due to entrainment is a function of 
the rate of liquid injection m^, the momentum parameter M 0> 
and the surface tension o, with the amount of liquid en- 
trained into the gas stream increasing with increasing 

values for m. and fi and decreasing values for a. 
l o 

6. The liquid viscosity was not found to influence the rate 
of entrainment, thereby contradicting the proposal by many 
observers that the liquid-film Reynolds number is a primary 
parameter characterizing the phenomenon of entrainment. 

7. The "effective" roughness of the liquid film can result 
in an increase in the evaporative transfer of mass of as 
much as 80 percent or more over that predicted by the simple 
theory for mass transfer. 

8. The data for mass transfer with a constant velocity gas 
stream was satisfactorily correlated in terms of two dimen- 
sional parameters, one characterizing the entrainment phe- 
nomenon and one characterizing the phenomenon of film sur- 
face roughness. The comparison of the experimental data 
with the present theory showed that the agreement was gen- 
erally to within 5 to 10 percent. 

9. There is some physical evidence that the dimensional 
entrainment group is equivalent to a dimensionless 
Weber number if the characteristic length for the latter 
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is taken as the scale of the small disturbance waves on the 
surface of the liquid film. Moreover, there is some evi- 
dence that the phenomenon of entrainment of liquid from the 
surface of a thin liquid film is analogous to the phenomenon 
of atomization of liquid drops suspended in a gas stream. 

10. ^-Tlie data for mass transfer with an accelerating gas 
stream was satisfactorily correlated by locally applying 
the correlation that was developed for the mass transfer 
data with a constant velocity gas stream and then numeri- 
cally integrating over the length of the liquid film to 
obtain the total mass transfer. 

11. The analysis that was developed to extend the results 
of the present investigation to the general case of an 
arbitrary film cooled length was shown to agree favorably 
with the experimental data of Refs. (19) and (20). This 
favorable comparison also suggests that the relatively sim- 
ple mass transfer correlation developed herein has consider- 
able utility because the flow parameters that characterized 
these two investigations were substantially different from 
those that characterized the subject investigation. 

12. The liquid coolant requirements for cooling the nozzle 
of a high energy propulsion device can be reduced substan- 
tially by injecting the liquid coolant at several axial 
locations, and by selecting coolants which have a high sur- 
face tension. 
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13. Correlation of the data for the maximum liquid tempera- 
ture showed that the simple theory developed herein pre- 
dicted that temperature with sufficient accuracy for most 
engineer!* ng- type calculations. 
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NOTATION 


acoustic speed (Eq. (3-69)) 

constant (see Eq. (3-1)) 

constant (Eq. (3-63)) 

constant (see Eq. (3-2)) 

heat-mass transfer parameter (Eq. (3-22)) 

mass concentration of vapor 

skin friction coefficient (Eq. (3-20)) 

specific heat at constant pressure 

specific heat at constant volume 
constant (Eq. (3-73)) 
tube diameter 

entrainment parameter (see Eq. (3-5)) 
denotes a function 
correction factor (see Eq. (3-16)) 
mass velocity for gas stream (Eq. (2-1)) 
specific enthalpy 

heat transfer coefficient (Eq. (1-2)) 
latent heat of vaporization 
thermal conductivity 
dimensional constant (Eq. (3-41)) 
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1 liquid-film cooled length 

1 10-inch long film cooled length 

l e measure of scale of waves on surface of liquid 

s film 

m rate of coolant injection or withdrawal 

m' rate of mass transfer from liquid film (see Eqs. 

(3-2) and (3-6)) 

m" local rate of mass transfer from liquid film 

m" rate of gas flow through the boundary layer per 

unit spanwise length (Eq. (5-6)) 

M molecular weight 

M_ Mach number 

O 

M q Momentum parameter (Eq. (2-2)) 

n parameter in the power-law velocity profile 

(Eq. (5-3)) 

Nu Nusselt number (Eq. (4-1)) 

p static pressure 

p° stagnation pressure 

p y partial pressure for vapor 

V c d 

Pr Prandtl number = — ^ 

q energy transfer due to conduction 

q f energy transfer due to radiation 

r roughness parameter (see Eq. (3-4)) 

R perfect gas constant 

Re D Reynolds number based on D 

Re x Reynolds number based on x 

Re x see Eq. (3-65) 
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Re' see Eq. (3-66) 

x t 

Re Q Reynolds number based on Q (Eq. (5-21)) 

St corrected Stanton number for simple heat 

5 transfer (Eq. (3-15)) 

St uncorrected Stanton number for simple heat 

so transfer (Eq. (3-23)) 

St' Stanton number for simple mass transfer 

s (Eqs. (3-85) and (3-86)) 

T static temperature 

T° stagnation temperature 

u streamwise component of velocity 

w weighting parameter (Eq. (4-2)) 

We Weber number (Eq. (3-38)) 

x distance measured from effective starting 

point of velocity boundary layer 

x' distance measured from point of liquid injection 

x t thermal entrance length 

X g entrainment group (Eq. (3-31) with a ■ Jj) 

roughness group (Eq. (3-45)) 

y distance normal to the wall or the liquid 

film surface 

Greek 

Y specific heat ratio 

6 dynamic boundary layer thickness 

thermal boundary layer thickness 

*k 

6 displacement thickness (Eq. (5-14)) 

6 momentum thickness (Eq. (5-15)) 

p dynamic viscosity 
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dimensionless entrainment parameter (Eq. (3-34)) 
den s i ty 

surface tension 

energy balance parameter (Eqs. (3-13) and (3-14)) 


Subscripts 


evaluated at the plane of liquid coolant 
injection 

evaluated at the terminus of the liquid-film 
cooled length 

evaluated at the terminus of the gas-vapor 
cooled length 

adiabatic wall 

critical 

effective 

main-stream state (y = 6) 
main-stream gas 
gas-vapor cooling 
summation index 
liquid 

evaluated at the gas-liquid interface but 
in the liquid phase 

maximum 

reference 

evaluated at the gas-liquid Interface 
simple 

evaluated at the gas-liquid interface but In 
the gas phase 


v 


vapor 
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w wall 

Superscripts 


o 


stagnation conditions 
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APPENDIX A 

AN ANALYSIS FOR THE GAS- VAPOR-COOLED REGION 
A. 1 Introduction 

It is difficult, if not impossible, to make a complete 
mathematical analysis of the boundary layer and heat trans- 
fer characteristics of the region immediately downstream 
from the end of the liquid film. Moreover, an involved 
analysis would not be generally warranted because of the 
large number of assumptions that must be made concerning 
th¥ vel^oci ty and temperature profiles, the form factor or 
factors, the skin-friction law, the heat- transfer law, etc. 

As a consequence of this fact, the analysis presented here- 
in is based on an idealized and relatively simple flow model. 
It will be shown that the analysis cannot be completely 
developed until more definitive experimental data are ob- 
tained for the subject problem. However, the limited data 
that do exist show that the analysis presented herein has 
considerable promise. These data are discussed In Sub- 
section A. 2.4. 

Both Emmons (20) and Guinn (45) developed seml- 
emplrlcal analyses for the gas-vapor-cooled region; these 
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analyses were based on the flow model suggested by Hatch 
and Papell (46) for gas-film cooling. The analysis pre- 
sented here differs from those due to Emmons and to Guinn 
in the following fundamental respects: 

(a) boundary-layer type flow rather than fully-developed 
pipe low is considered; 

(b) the analysis is generalized so that it can be applied 
to a compressible gas flow in the presence of a streamwise 
pressure gradient; and 

(c) no direct similarity between gas-vapor and gas-film 
cooling is assumed. 

As discussed in Subsection 2.1.2, an analysis for the 
gas-vapor-cooled region requires the consideration of two 
cases: (a) the case where the wall is adiabatic; and (b) 

the case where the wall is externally cooled. These two 
cases are considered in Sections A. 2 and A. 3, respectively. 

A. 2 The Case Where the Mall is Adiabatic 

A. 2.1 The Flow Model 

The flow model employed for the analysis of the gas- 
vapor-cooled region is defined by the following conditions. 

(1) The wall is adiabatic. 

(2) The boundary layer flow is fully turbul**"*. 

(3) The thermal, the diffusional, and the dynamic boundary 
layer thicknesses are equal; the latter is denoted by y = 6. 


% 
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(4) The stagnation temperature for the main stream, T°, is 
a constant in the streamwise direction. 

(5) The specific heat for the gas, Cp, is a constant. 

(6) Radiation heat transfer and chemical reactions within 
the boundary layer are not important. 

(7) The velocity profile and the stagnation- temperature 
profile for the boundary layer are of the function forms: 

5~ ■ f i<|> 

9 


and 


• f 2<*> <*-*> 

‘g aw 

where T aui is the temperature of the fluid adjacent to the 

a W 

adiabatic wall. To illustrate, f^ could be assumed of the 
power-law form: 


f x = (£) 1/n ( A- 3 ) 

and fg could be assumed of the form 

fj 'At B(J) + C($) + 0(f) 3 (A-4) 


The constants in Eq. (A-4) can be evaluated from the 
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following boundary conditions that apply to f 2 . Thus 
f 2 (0) - 0 

f 2 (1) = 1 

3f / 


9y 

3f / 

9 ^ 


- (0) - 0 


( 1 ) = 0 


The resultant expression for f 2 is 


f? = 3(^) - 2(^) 


( A- 5 ) 


For the present analysis, however, the actual forms of f, 
and f 2 need not be assumed , only the assumption that f^ and 
f 2 exist for the gas-vapor-cool ed region is required. 

A. 2. 2 The Energy Balance 

Figure A.l presents, schematically, a control volume for 
the gas-vapor-cool ed region; the control surface of the con- 
trol volume is defined by x = x 2 , x = x, y = 5, and y = o. 

The following energy fluxes, written in accordance with the 
flow model defined in the foregoing section, are presented 
In the figure: 

(a) j pu Cp T° dy a the rate at which thermal and ki- 

0 netic energy is transferred Into the 

control volume per unit spanwlse 
length at the station x * x 2 . 


V T •«! >$' " ' V 4 


'«T 


'*■*"*** 



mainstream 



o <y 

5 i 

o o 

o > 


- i l 


Fig- A.l Flow Model for Analysis of Gas-Vapor Cooling. Adiabatic Wall 
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6 

(b) j pu Cp T° dy = the rate at which thermal and ki- 

0 netic energy is transferred out of 

the control volume per unit span- 
wise length at the station x = x. 

(c) (m" - m!i) c„ T° = the rate at which thermal and ki- 

c p g 

netic energy is transferred from 
the main stream into the control 
volume per unit spanwise length. 

The parameter m" in term (c) is defined as follows: 

m" = | pu dy = the rate at which mass is trans- 
0 ferred through the boundary layer 

per unit spanwise length at the 
station x * x 

( A- 6 ) 

The energy balance for the control volume in Fig. A.l 

can be written as follows. Thus 

6 2 S 

| pu c p T° dy + (m" - mj) c p T® - j pu c p T° dy 
o o 


(A-7 ) 
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Employing Eq. (A-6), and dividing through by Cp, Eq. (A-7) 
can be written in the form 

6 2 6 

| pu T° dy | pu T° dy 

m 2 r 2 + (m " - m V T g = m " I 

f pu dy j f> u 

b 0 

(A- 10) 

In order to realize the desired simplicity of the present 
analysis, the following approximation is introduced con- 
cerning the integral terms in Eq. (A-10). Thus* 


6 

j pu T° dy 

o 

6 

j pu dy 
o 


6 

j u T° dy 

*\/ o 
55 6 

j u dy 
o 


Thus, Eq. (A-10) can be written in the following form by 
introducing the above approximation and assumption (7) of 
the flow model : 

*Note that this approximation is not equivalent to and Is 
less restrictive than the assumption of constant density. 


SJT 


- 
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c. 

\ h << T g - T aw, 2 > f 2 + T aw,2> d * 

0 

6 2 

j f i d * 


+ (m" 


■i) T 9 


o 

| f l << T g - T aw> f 2 * T aw> d * 
m " “ 5 

j f i d * 

o 

The latter expression can be rewritten In the form 


"2 < T ? - T aw,2> 


j f l f 2 d <$> 


"2 < T g - T aw,2> I 


| f l d <«> 


J f l f 2 <> 


I" (Tq - T aw ) - -■ (Tj - T aw ) J. 


j »i d <S> 


Finally, the terms In Eq. (A- 11) can be combine- tu yield 
the following energy balance for the control volume pre- 
sented In Fig. A.l. Thus 
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(A -12) 


The relatively simple form of Eq. (A-12) is due principally 
to assumption (7) of the flow model. 

Equation (A-12) enables the determination of the 
adiabatic-wall temperature T as a function of x once T aui 9 
and T° are specified and m" as a function of x is known. 

The parameter T° is normally known from the statement of the 
problem, and the parameter T aw 2 can he estimated from the 
developments presented in Subsection 3.4.1. The evaluation 
of m" as a function of x is considered in Subsection A. 2. 3. 


A. 2. 3 The Evaluation of m" 

Equation (A-6), which defines the parameter m", can be 
rewritten in the form 


■’ * <P“)g <f - f-> 6 CA-13) 

where 

4 

4* * | (1 - ) <*y « the displacement thickness 

0 9 


(A-14) 
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6 

f 


ou 


j ‘vput; 


(i - 


r> d * 
g 


the momentum thickness 


(A- 15 ) 


Equation ( A- 13 ) shows that m" can be determined as a function 
* 

X X 

of x if and C are known as functions of x. The eval- 

uation of these quantities is considered in the following 
paragraphs . 

The Evaluation of 6/6 and 6*/8 . Persh and Lee (47) 
have tabulated the values of 6/6 and 6*/6 as functions of 
the main-stream Mach number and the wal 1 -to-main-stream 
temperature ratio. It is suggested, therefore, that these 
tabulated data be employed for evaluating 6/6 and 6*/8 when 
the variation in the density for the gas stream across the 
boundary layer Is significant. However, for the more general 
case where that density variation Is not significant, it is 
suggested that the subject parameters be calculated by 
applying the power-law velocity profile (see Eq. (A-5)). The 
resultant expressions for 6/8 and 6*/6 are (48, p. 536) 


6 

6 


(1 + n) (2 + n) 


(A-16) 


and 


6_ 

e 


[2 t *1 
n 


(A- 17) 
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The Evaluation of 9 . To determine the value of the 
momentum thickness 9 as a function of x, one must solve the 
momentum integral equation. That equation can be written in 
the form (49) 

(A-18) 


(A-19) 

where 




the skin-friction coefficient 


( A- 20 ) 


To solve* either Eq. (A-18) or Eq. (A-19) for 6, one 
must specify the form of the skin-friction law; the latter 
is denoted functionally by 



*A numerical solution is generally required for either 
Eq. (A-18) or (A-19). 
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where 


p u 0 

Re = -3 — 2 — = the Reynolds number for the main 

9 

stream based on the characteristic 
length 9 

Several empirical skin-friction laws are presented in 
Chapter 22 of Ref. (48) for the case where the gas flow is 
turbulent, incompressible, and over a flat plate. In the 
analysis of gas flows with favorable pressure gradients, 
it is normally assumed that the skin-friction law is 
the same as that established empirically for the flat- 
plate case (49, 50, 48). Furthermore, the case where 
heat transfer effects and/or compressibility effects are 
important can be analyzed by evaluating the thermodynamic 
properties for the gas stream at an appropriate reference 
temperature. For example, the reference temperature, denoted 
by T r , can be determined from (36): 

T r * 0.5(T + T ) + (Pr) 1/3 J (A-21) 

g w rc“ ~ 

P9 

It is suggested that in the evaluation of the physical 
properties that the composition of the gas be assumed uni- 
form throughout the boundary layer and equal to that for the 
main stream; i.e., assume that a reference temperature is 
sufficient to define the reference state. 

To complete the solution for e, the boundary condi tion 



229 


on 6 at the station x = Xg, denoted by must be specified. 
An approximate technique for the evaluation of that boundary 
condition was proposed in Ref. (25). However, subsequent 
analysis of the available experimental data suggests that 
that proposal is not satisfactory. In order to illustrate 
that point, an analysis of the data that are presented in 
Ref. (19) is discussed in Subsection A. 2. 5. However, to 
facilitate that discussion, it is convenient to first con- 
sider a special case for which the foregoing developments 
are simplified substantially. 

A. 2. 4 A Special Case 

Consider the special case that is defined by the 
following conditions. 

1. The gas flow is incompressible, constant density, and 
over an adiabatic flat plate. 

2. The skin-friction law is that due to Prandtl (48, Eq. 
22.7). Thus 

= 0.0128 Re Q ~ h (A-22) 

3. The velocity profile for the boundary layer is given 
by Eq. (A-3). Thus 



(A-3) 



230 


For this special case, Eq. ( A- 12 ) reduces to the fnrm 


aw 


T g - W 


mj; 

m ir 


( A- 23 ) 


Equation ( A- 13 ) reduces to (refer to Eqs. ( A- 16 ) and ( A- 17 ) ) 


m" = (' ' n) p u g 6 


Equation ( A- 18 ) 


_ de 

2 cfx 


simplifies to the form 

d Re e 

dR % 


where 


( A- 24 ) 


( A- 2 5 ) 


the Reynolds number for the main 
stream based on the characteristic 
length x 

(A- 26 ) 

Note that due to assumption (1), the distance x is measured 
from the "effective" starting point for the growth of the 
turbulent boundary layer over the gas-vapor cooled region. 

It follows from the discussion presented below that the 
specification of that starting point is equivalent to the 
specification of the parameter Gg. 



wm 
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Upon combining Eqs. ( A- 2 3 ) and ( A - 24 ) , one obtains 




( A— 27 ) 


In addition, substitution of Eq. ( A - 2 2 ) into Eq. ( A - 2 5 ) 
yields 


0.0128 dRe v = ReJ dRe fl 
x y d 


(A-28 ) 


Integrating Eq. (A-28) and :ubstituting the limit Re Q = Re ( 
at x - Xg, one can readily show that 


Re 6 9 0« x 9 

£ = _£. = (—) 

Re 0 0 ' x ' 


4/5 


( A- 29 ) 


Substituting Eq. ( A- 29 ) into Eq. ( A- 27 ) results in 




( A- 30 ) 


Equation (A- 30 ) presents, for the subject special case, a 
simple expression for the calculation of the adiabatic-wall 
temperature T... as a function of the distance from the end 
of the liquid film (x - Xg). The problem, therefore, is to 
determine x 2 . The data of Ref. (19) give some insight into 
the problem of the evaluation of that parameter. Those data 
are considered in the following subsection. 




m 
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A. 2. 5 Experimental Data for the Adiabatic-Wall 
Temperature From Reference (19) 

Kinney, Abramson, and Sloop (19) present some data for 
the temperature for an adiabatic wall immediately downstream 
of a liquid-film cooled wall region*. To the author's 
knowledge, these are the only such data that appear in the 
published literature. The experimental apparatus that was 
employed in the subject investigation was discussed previ- 
ously in Subsection 4.3.2. The reported data were obtained 
with the 4-inch diameter smooth surface tube at a nominal 
gas stream temperature of 1400F and a nominal gas stream 

5 

Reynolds number of 5.3 x 10 . The length of the liquid-film 
cooled region was approximately 19.0 inches. 

To obtain an evaluation of the present analysis, it is 
assumed that Eq. ( A- 3 0 ) that was developed for the afore- 
mentioned special case is applicable to the subject data. 

The determination of the value for (and thus the effec- 
tive starting point for the growth of the boundary layer 
over the gas-vapor cooled region) for these data gives some 
insight into the phenomena involved. 

Figure A. 2 presents a plot of the left-hand side of 
Eq. (A-30) as a function of (x - Xg), the distance measured 
from the end of the liquid film. The circles shown in the 
figure were calculated from the data presented in Fig. 4 of 
Ref. (19). The solid curves shown in Fig. A. 2 were calculated 


♦See Figure 4 of Ref. (19). 
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from Eq. (A-30) with values for *2 °f 0* 8 > anc * 1,0 

inches. For the wall region immediately downstream of the 
liquid film (say, x - X 2 < 10.0 inches), Fig. A. 2 shows that 
a value of 0.8 inches for X 2 results in an excellent corre- 
lation of the experimental data. The small value of 0.8 
inches for Xg* however, is difficult to interpret physically. 
This small value suggests that the boundary layer phenomena 
are complex, with perhaps a considerable change in the 
boundary layer structure occurring at (or very near) the 
terminus of the liquid film. The basic analysis, however, 
shows considerable promise. Even the complex case where a 
compressible gas is strongly accelerated over the gas-vapor 
cooled region can be reduced to the problem of determining 
the "effective starting point (x = 0 ) for the growth of 
the boundary layer over the subject wall region. The task 
that remains, therefore, is to obtain experimental data on 
the adiabatic wall temperature for a wide range of flow 
conditions so that the effective starting point can be 
correlated in terms of the pertinent variables. Unfortu- 
nately, since such data are not now available, the analysis 
proposed herein for the gas-vapor cooled region cannot be 
completely formulated. 




For the subject case, the basic analytical problem is 
that of relating the local values for the wall temperature 
and the wall heat flux to the main-stream conditions. It 


ca rii be inferred from several mass transfer cooling (and 
heating) studies that are reported in the literature (e.g.. 
Refs. 51-55) that heat- transfer data for the gas-vapor 
cooled region should correlate satisfactorily in the follow- 
ing manner. Thus 

St gv = St O (A " 31> 


where 




% 


(A -32) 


St = the Stanton number for heat transfer for the 
9 v gas-vapor cooled region 

St = the Stanton number for heat transfer in the 
0 absence of gas-vapor cooling, assuming that 
the heat transfer begins at the leading edge 
of the gas-vapor cooled region 

T * the temperature the wall would obtain if it 
aw was adiabatic 

T = tlie temperature the wall [obtains due to the 
w fact that it is cooled externally 

q = the heat flux from the hot gas stream to the 
w external coolant ^ 

Equation ( A-31) reduces the present problem to that of 
determining the Stanton number St Q as a function of the 
distance measured downstream from the end of the liquid fil 


The theoretical methods that were outlined in Subsection 
3.1.2 and in Subsection 3.3.1 can be employed to calculate 
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that Stanton number*. In that calculation, however, it is 
important to note that the following boundary condition 
applies (51-55)**: 

St =«> at x = x 0 
o c 

The task remains to obtain experimental data for the 
rate of heat transfer to an externally cooled surface that 
is located immediately downstream of the terminus of the 
liquid film so that the validity of Eq. (A-31) can be deter- 
mined. Although Emmons (20) and Guinn (45) did obtain some 
experimental data on that rate of heat transfer, the validity 
of Eq. (A-31) cannot be investigated with those data because 
neither Emmons nor Guinn determined the adiabatic wall tem- 
perature T flw that enters into the equation. There is, how- 
ever, an important observation that can be made with ref- 
erence to the experimental results due to Emmons. That ob- 
servation is discussed in the following subsection. 

* With F r = 1.0 and F * 1.0 

L 511 

**Guinn (45) developed a semi-empirical analysis for the 
gas-vapor cooled region. Guinn, in essence, utilized 
Eq. (A-31) in that analysis. However, he did not 
recognize that the above boundary condition applied In 
the evaluation of the Stanton number St,,. As a con- 

sequence, his analysis is not valid at the leading edge 
of the gas-vapor couied region. 
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/\ , 4 An Observation on the Data of Emmons_ 
for Gas -Vapor Cooling 

Emmons (20) investigated experimentally the rate of 
heat transfer to an externally cooled test section that was 
located immediately downstream of the liquid-film cooled 
test section. Those heat transfer rates, which were deter- 
mined cal ormetri cal ly at several locations downstream of the 
film cooled region, werfa compared with the heat transfer 
rates that were realized for the aforementioned 1. 48-inch 
convectively cooled test section that was located immedi- 
ately upstream of the liquid film*. As a result of that 
comparison, Emmons drew the following two conclusions (p. 67, 

Ref. (20)): ’ - ' ; 

•'A large portion of the film coolant vapor forms a 

mass of cool gas which remains adjacent to the wall for a 
considerable distance downstream from the liquid film. As 
a result the heat flux to the wall protected by the vapor 

is significantly reduced;" and 

"Heat transfer rates measured approximately 9,0 in. 
downstream from the liquid film were reduced by 30 to 70 
percent from the corresponding values measured upstream 
from the liquid film. Ammonia, with the lowest boiling 
temperature (160F) effected the largest heat transfer re- ^ 
duct ions for the cool ants investigated. 

These two conclusions imply that the rate o f mixing of 

*See Fig. 42. f.. > .hf '•¥ : 
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the relatively cool boundary layer gases with the mainstream 
gas is slow. However, the data of Ref. (19) that were dis- 
cussed in Subsection A. 2. 5 for the adiabatic wall temperature 
suggest that the rate of mixing is quite rapid. That is 
implied by the small value for (or e g ) that was found to 
result in a satisfactory correlation of the experimental 
data. That small value for x^ suggests that the rate of 
entrainment of the hot mainstream gas into the boundary 
layer (or the rate of growth of the boundary layer) is very 
rapid for the subject wall region. The data presented in 
Fig. A. 2 show that at a point approximately 10 inches down- 
stream of the liquid film, the cool boundary layer gas and 
the hot mainstream gas are essentially completely mixed. 

There appears to exist, therefore, some Inconsistency in the 
aforementioned conclusions that were drawn by Emmons and the 
experimental data that were reported in Ref. (19). However, 
an explanation for that inconsistency crn be suggested: It 

was argued in Subsection 4.4.5 that the measured rates of 
heat transfer to the 1.48-inch convectively cooled test 
section were apparently about 40 percent higher than those 
that would be predicted for fully-developed tube flow at the 
sa.re bulk flow conditions. Moreover, it was argued that 
part (or perhaps all) ♦ha*. difference could have been due 
to a discontinuous thermal boundary condition at the leading 
edge of the 1.48-inch test section. For the convectively 
cooled test section immediately downstream of the end of the 
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liquid film, however, there would be a continuous transition 
in the wall temperature from the end of the film cooled test 
section to the gas-vapor celled region. Therefore, because 
of the fundamental difference between the thermal boundary 
conditions for those two test sections, it would appear that 
the aforementioned conclusions that were drawn from the 
direct comparison of the heat transfer rates for the two 
test sections may not have been completely valid. If, in 
fact, an i sothermal -wal 1 boundary condition for the 1.48- 
inch long test section would have resulted in heat transfer 
rates that were 40 percent lower than those that were actu- 
ally realized, then the comparison of the experimental data 
for the two convectively cooled test sections would have 
showed that the heat transfer rates measured approximately 
9.0 downstream of the liquid film would have been reduced 
by 0 to 20 percent from the corresponding values measured 
up. tream from the liquid film, as compared to the 30 to 70 
percent reduction that was quoted above from Emmons (20). 
That particular result would appear to be in agreement with 
the data that were discussed in Subsection A. 2. 5. However, 
it cannot be argued conclusively that the significant con- 
clusions that were drawn by Emmons with respect to the sub- 
ject wall region are in error. That question cannot be 
completely answered until considerable more experimental 
data are obtained under well controlled experimental con- 
ditions for the gas-vapor cooled wall region.. Perhaps the 
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foregoing discussion will serve as an impetus for that 
research . 
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APPENDIX B 

DESCRIPTION OF APPARATUS 
B . 1 General Description 

The gas source employed in the subject investigation 
was an air blowdown wind tunnel in which a portion of the 
air flow was heated by combustion with a liquid fuel. The 
principal features of the apparatus were that it had a wide 
range of operating conditions, i.e., flow rate, gas temper- 
ature and pressure, reproducibility of operating conditions 
from run to run, low drift about the operating point over 
run times up to one hour, and was safe and reliable enough 
to work with in the test cell during operation of the 
apparatus . 

Figure B.l presents a photograph of the basic con- 
figuration of the experimental apparatus for that phase of 
the investigation where the gas stream velocity was essen- 
tially constant over the length of the wetted test plate. 
Figure B.2 presents a similar photograph of the basic con- 
figuration that was employed in that phase of the investi- 
gation where the gas stream was accelerated over the wetted 
test plate. The various components of the apparatus are 
described as follows. The numbers refer to either Fig. B.l 
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or Fig . B. 2. 

(1) The air and fuel flow rate control valves. A separate 
valve was utilized to control (a) the air flow rate to the 
primary combustor, (b) the air flow rate to the secondary 
mixer, (c) the fuel flow rate to the primary combustor for 
ignition, and (d) the fuel flow rate to the primary com- 
bustor during full-scale operation. 

(2) The primary combustor wherein a high- temperature com- 
bustion zone was established. 

(3. The secondary mixer wherein diluent air was added to 
the hot gases from the primary combustor so as to obtain 
the desired gas stream temperature in the experimental 
tunnel . 

(4) The 10-inch long mixing section which contained one 
10-mesh screen and one 30-mesh screen that were spaced 10 
inches apart, and which reduced the non-homogeneities of 
the gas stream flow substantially. 

(5) The 36-inch long adiabatic approach section. 

(6) The experimental test section which contained the 
film cooled test plate. 

(7) The variable area nozzle which controlled the gas 
stream pressure for a given gas stream Reynolds number. 

(8) The 35mm Nikkorex camera that was employed to photo- 
graph the surface characteristics of the liquid film. The 
camera lens was a Micro-Nikkar f5Emm lens with a close-copy 
adapter. 
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(9) The stroboscopic light source that was employed in 
conjunction with the camera (8). The strobe had a flash 
duration of 3 microseconds. 

(10) The cyclone separator that was employed to separate 
the liquid and the gas phases that were withdrawn at the 
capture slot. The details of the separator design is given 
in Appendix B.4. 

(11) The vent valves that controlled the rate at which the 
gas phase was vented to the surroundings (and, therefore, 
the rate at which the gas phase was withdrawn at the capture 
slot) . 

(12) The electronic timer tha>- was employed to measure the 
sample time for the mass transfer measurements. 

(13; The 500 psia air solenoid valve that served as the 
on-off control for the capture slot. The activation of the 
electric solenoid valve resulted in the flow of a high 
pressure jet of air in the withdrawal passage in the reverse 
direction through the capture slot into the test section, 
thereby stopping the capture of the liquid film on the test 
plate. This solenoid valve and the electric timer (12) were 
controlled simultaneously so that the length of time that 
the gas-liquid mixture was withdrawn from the test section 
through the capture slot was accurately controlled. 

(14) The heat exchanger that was employed to preheat the 
liquid phase before it was introduced onto the test plate. 
The heat exchanger consisted of two 3/4-in. O.D. stainless 
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steel tubes of length 36 inches with a calrod heater located 
along the centerline of each. Each calrod had a maximum 
capacity of 1000 watts, and the electrical energy that was 
supplied to each was controlled by a variable reostat. 

(15) The mercury manometers that were employed to measure 
the streamwise static pressure profile for those experi- 
mental tests in which the gas flow was accelerated over the 
wetted test plate. 

(16) The film coolant flowmeter i was employed to mea- 
sure the rate at which the coolant was introduced onto the 
test plate. The flowmeter was a 1/8-inch Potermeter. The 
frequency output of the flowmeter was displayed to the 
operator by the milliampmeter (19) and recorded on a strip- 
chart recorder. 

(17) The air flow meter. The total rate of air flow through 
the experimental apparatus was measured with a Brooks tur- 
bine type flow meter (model HP-321-6). The frequency output 
of the flow meter was displayed on a milliampmeter in the 
control room and was simultaneously recorded on a stripchart 
recorder. 

(18) The precision needle valve that was employed to con- 
trol the rate of film coolant flow through the heat ex- 
changer and then onto the test plate. 

(19) Milliampmeter utilized to monitor the coolant flow 
rate. 
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Figure B.3 shows the control room. The primary com- 
ponents of the control room are as follows. The numbers 
refer to Fig. B.3. 

(1) The automatic sequencer that was utilized to start 
the hot gas generator. The operation of the sequencer is 
discussed in greater detail in Appendix B.2. 

(2) The 3 small volume Nullmatic pressure regulators that 

were employed to control the flow rates through the three 
Annin Domotor control valves: (a) the primary air control 

valve, (b) the secondary diluent air control valve, and (c) 
the primary fuel control valve for ^ull-scale operation. 

(3) The stripchart recorder that was employed to continu- 
ously monitor the gas stream temperature. 

(4) The Bourdon tube pressure gage that was employed to 
continuously monitor the gas stream pressure. 

(5) The variable reostat that was utilized to control the 
temperature for the liquid coolant at the point of injection. 

(6) The mi 11 iampmeters that continuously m r nitored the 
pertinent system flow rates and system pressures. These 
various parameters were also recorded on stripchart re- 
corders located in an adjoining room. 

B.2 The Gas Generat or 

Figure B.4 n*®?®nts a simplified section* 1 drawing of 
the gas generator portion of the experimental apparatus. 




Fig. B.3. Photograph of Control Room 
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Secondary 






Table B.l lists the components of the gas generator 


and describes the function of each. 


Table B.l 


Components of the Gas Generator 


Component 
Primary Injector 


FI amehol der 


Primary Combustor 


Functio n 

To inject the fuei and air to 
be burned Into the primary 
combustor. 

To mix the fuel and air after 
injection and provide tur- 
bulence to hold the flame. 

To provide a contained volume 
for combustion. 


Secondary Injector To inject diluent air to re- 

duce temperature of working 
fluid to the desired value. 

Perforated cone and To thorough - ^ mix th„ diluent 

turbulence grids air with the combustion 

products from the primary 
combustor in order to estab- 
lish a uniform temperature 
throughout the gas flow. 
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Table ^.2 presents the output characteristics of the 
gas generator for the overall allowable range of operation 

Table B.2 

Output Characteristics of the Gas Generator 



Parameter 

Overall Range 

Gas 

Stream Pressure, psia 



Cold flow 

up to 400 


. Hot flow 

50 - 350 

Gas 

Stream Temperature, R 



Cold flow 

500 


Hot flow 

750 - 1100 

Gas 

Flow Rate, lb/sec 



Cold flow 

up to 15 


Hot flow 

0.75 - 10.0 

Fuel 

Flow Rate, lb/sec 

0.02 - 0.25 
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The air flow entering the gas generator was divided 
such that approximately 30 percent of the total flow entered 
the primary combustor to be burned with the fuel while the 
remaining 70 percent of the air flow was injected into the 
secondary mixing chamber. Approximately five percent of the 
fuel was injected from the ignition fuel injector located in 
the center of the primary injector while the remaining fuel 
was injected fro;, the main portion of the primary injector. 
The fuel and the primary air (representing approximately 200 
percent excess air) were burned in the primary combustor at 
a temperature of 2500F with a thermal efficiency of approxi- 
mately 50 percent. Smooth stable combustion rather than 
high thermal efficiency was desired. The low thermal effi- 
cency was the result of insufficient volume in the primary 
combustor to permit complete combustion of the fuel. 

The primary injector comprised three principal parts, 
as shown in Fig. B.4, namely, the face plate, the hat section 
air manifold and the cap plate. Ninety-five percent of the 
methyl alcohol fuel was injected through two concentric 
rings having a total of 30 separate fuel injector buds 
located on the face plate (see detail in Fig. B.2). The 
face plate was fabricated from a circular piece of stainless 
steel plate h in. thi>k On each side of the fuel injectors 
were two 0. 100-in. diameter injector orifices for air di- 
rected at 45° to the direction of the fuel flow, the com- 
bination forming a triplet, The primary air was introduced 
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into the injector .manifold by a single 1.0-in. tube fitting. 
The fuel was manifolded by %-in. O.D. stainless steel tube 
to six clusters of heavy wall tubing, each with tubes welded 
into the back of the face plate; these clusters fed fuel to 
the 30 fuel injectors. The ignition fuel which was required 
as part of the start sequence was manifolded directly to the 
center of the face plate and was injected through four small, 
orifices against a copper splash ring. Air for the com- 
bustion of the ignition fuel was supplied through a ring of 
small orifices concentric to the ignition fuel injector 
orifices. 

The flameholder consisted of a central perforated disc 
of stainless steel supported by eighteen ^-in. tubes of the 
same material which ran radially to an outside ring. Cool- 
ing water was introduced through a 1.0-in. tube fitting 
located 180° from the output fitting. The flameholder was of 
all welded construction. In the assembled gas generator the 
flameholder was sandwiched between the primary injector and 
the primary combustor. 

The primary combustor and secondary mixer chambers were 
of similar design and construction, differing only in length. 
The former was 10 in. long and the latter was 20 in. long. 

The inner wall was a section of 13-in. I.D. mild steel pipe 
having a wall thickness of h in. Welded at each end of the 
pipe section were one-inch thick mild steel flanges with 
twenty-four l-1n. nominal bolt holes equispaced on a bolt 
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circle of 20 in. The coolant passages were formed from ex- 
truded lead stripping with a %-in. by 3/4-in. cross section 
which was wrapped in a spiral around and soldered to the 
pipe. The outer jacket of the coolant passages was formed 
from half-sections of 1/8-in. wall, 16-in. I.D. stainless 
steel tube with an expansion joint located in the middle. 

The jacket was sealed all around by weldments. Inlet and 
outlet fittings for the coolant, which was water, were 
welded over holes in the outer jacket and located at either 
end of the lead spiral. All cooled parts were incorporated 
in a single pass coolant path. 

The secondary air injector was fabricated from seven 
parallel 3/4-in. O.D. stainless steel tubes set in a 13-in. 
I.D. outer ring which served as a distribution manifold. 

Each tube was drilled with a 0.125-in. O.D. hole every 3/8 
in. along .its length in such a way that two rows of orifices 
were formed directed downstream 45° off from the direction 
of gas flow from the primary combustor (see detail in 
Figure B.2). This arrangement enhanced turbulent mixing 
immediately downstream from the injector tubes. The injec- 
tor tubes were fed from the outer ring which in turn was 
fed secondary air through a 2-in. pipe fitting located on 
the outermost part of the ring. 

Downstream of the secondary air injector was located 
a stainless steel cone perforated with *s- 1 n . diameter holes. 
The cone pointed upstream and was clamped between the end 
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plate and the aft flange of the secondary mixer chamber. 
Turbulence grids were mounted just ahead of the 2-in. x 4-in. 
exit port located in the center of the end plate. The per- 
forated cone and turbulence grids induced mixing of the 
secondary and primary gases. The end plate was fabricated 
from 1-in. thick mild steel plate reinforced by gussets 
welded on the outside. 

The flanged sections of the gas generator were bolted 
together with twenty-four 1-8 studbolts. The gaskets be- 
tween the separate sub-assemblies were either 1/32-in. thick 
aluminum or 1/16-in. thick composition asbestos. Concentric 
serrations to aid sealing at the gaskets were machined on 
the injector face plate, the end plate and the flanges; no 
serrations were present on the flameholder or secondary air 
injector. 

Operation . Figure B.5 presents a simplified schematic 
diagram of the flow system for the experimental apparatus. 

The air used in the experiments was drawn from the high 
pressure supply tanks which comprise the blowdown air supply 
for the Jet Propulsion Center. The tanks were charged with 
a maximum pressure of 2400 psig. The compressed air was 
delivered to an isolating valve located on the roof of the 
test cell by a pipeline from the supply tanks. The Iso- 
lating valve on the roof was a 2-1n. 1500 ASA Annin pneu- 
matically operated open/closed valve. It was present as a 
safety feature to isolate the high pressure air from the 


high pressure 
air supply 
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test cell housing the apparatus in the event of a severe 
mechanical malfunction of the gas generator. 

The air supply was passed through a Fluid Dynamics high 
pressure air filter. The filter element removed particles 
from the gas stream that were larger than 10 microns. The 
air was then reduced from the supply line pressure by a 2-in. 
Grove Powereactor dome regulator located immediately down- 
stream from the air filter. The Grove regulator was loaded 
by means of a small volume Grove pressure regulator. The 
volume flow rate of air was measured with a 2-inch Brook's 
turbine flowmeter located after a straight section of pipe 
downstream of the regulator. The weight flow rate of air 
was determined from the frequency output of the flowmeter 
and the static pressure and temperature of the air measured 
immediately ahead of the flowmeter. The static pressure 
was measured with a Wiancko static pressure transducer and 
the temperature was measured with an 1 ron-constantan ther- 
mocouple. The three aforementioned parameters were re- 
corded during an experiment on Brown self-balancing strip- 
chart recorders. 

The air flow to the gas generator was controlled by 
two Annin Domotor valves, a 2-in. valve for the secondary 
air and a 3/4-1n. valve for the primary air. Control air 
for these pneumatic valves as well as for other pneumatic 
functions In the system was supplied at pressures of 400, 

100, 45, and 30 pslg from either Grove Powereactor 


258 


regulators or Nullmatic small volume regulators. 

The methyl alcohol fuel was stored in 55-gallon ship- 
ping drums and transferred to two 6-cu. ft. run tanks located 
outside the building housing the test cell. The fuel was 
pressurized with air by means of a small volume Grove pres- 
sure regulator. The fuel tanks were isolated from the 
building by a pneumatically operated Annin valve located 
beneath the fuel tanks. The flow rate of fuel was measured 
with a 1/8-in. Potter turbine flowmeter, the frequency out- 
put of which was recorded on a Brown recorder. The fuel 
flow rate was controlled by means of a modified 3/4-in. 

Annin Domotor valve which was rebuilt to control with pre- 
cision the relatively small flow rate of fuel. 

Before the gas generator was started in hot flow, the 
air pressure and fuel tank pressure were set to the desired 
value. Coolant water was supplied to the gas generator by 
a Worthington turbine pump driven by an electric motor. The 
three Annin Domotor valves were pre-set using Nullmatic 
control regulators (see Fig. B.3) to give the proper start 
flow rates of primary air, secondary air and fuel. 

The ignition source for the gas generator was a small 
combustor (see Fig. B.4) which burned gaseous methane and 
air ignited by du-’l r ;*rk plugs. 

The ignition sequence was controlled by an automatic 
sequencer (see Fig. B.3). The ignition sequence was as 
follows: 
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1. The methane-air ignitor was started. 

2. When the ignitor reached the desired temperature, 
the automatic sequence was initiated. 

3. The ignition fuel and primary air valves opened 
simultaneously and a pilot flame was established 
in the center of the primary injector off of the 
ignition fuel injector. 

4. Two seconds later, the primary fuel and secondary 
air valves opened simultaneously. 

5. The ignitor combustor was shut down and a second 
or two later the mainstage Ignition was accom- 
plished by the pilot ignition flame. 

The whole automatic start sequence required approxi- 
mately six or seven seconds. Once mainstage ignition was 
accomplished (usually at reduced flow rates), the operator 
brought the gas generator up to the desired run point using 
the controls for the Annin Domotor valves. 

During operation of the gas generator, monitoring of 
the various pressures associated with the system was accom- 
plished either by Bourdon tube pressure gages or Wiancko 
pressure transducers. The outputs of the latter were dis- 
played to the operator by ml 1 1 iampmeters* and were also 
recorded on Brown stripchart recorders. Flow. rates of fuel 
and air were also displayed to the operator by milliamp- 
meters* and recorded on stripchart recorders. The temper- 
ature of the gas leaving the gas generator was measured by 
two chromel-al umel thermocouples located immediately down- 
stream of the end plate at the entrance to the approach 
section. The output of one of the thermocouples was 
•Shown In Fig. B.3. 
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displayed to the operator on a direct reading mi 1 1 i ampmeter 
while the other output was recorded on a stripchart recorder. 
The pressure of the gas in the lead section was the value 
employed to set the pressure of the run point. The gas 
stream pressure was measured, displayed and recorded in the 
manner as previously described. 

After the desired operating point had been achieved, 
approximately ten minutes were required for the test section 
to achieve thermal equilibrium. Only then were the experi- 
ments begun. Luring run point operation the gas generator 
exhibited a random drift about the run point of ± 20F of gas 
stream temperature and ± 2 psi of gas stream pressure. 

B . 3 The Test Section and Te s t Plat e 

Figures B.6 and B.7 present sectional drawings of the 
experimental test section. The basic Internal dimensions 
of the test section were 2 in. by 5.5. in. cross section by 
18 in. long. The test section consisted of 4 side wal 1 * and 
the wetted test plate. The side walls were made from 9/16- 
In. thick stainless steel plate and the upper and lower 
walls were made from 3/4-1n. thick stainless steel plate. 

The test plate was *s-1n. thick stainless steel plate. All 
of the Internal wall surfaces were surface ground. A single 
Pyrex window of dimensions 11/16 x 1-5/16 x 9-7/8 In. long 
was mounted In the top wall of the test section, and 6 cir- 
cular l-1n. O.D. Pyrex windows 9/16-1n. thick were equi- 
spacetfat 2 inches in each side wall.. All of the Pyrex 
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windows were sealed by Viton-B o-rings, as indicated in the 
drawi ngs . 

The dashed line shown in Fig. B.6 represents the pro- 
file for the aluminum wedge for that phase of the investi- 
gation where the gas stream was accelerated over the wetted 
test plate. 

Figures B.6 and B.7 show that the test plate to which 
the liquid film was confined dropped into a recess that was 
machined into the bottom wall of the test section and 
"floated" on an o-ring seal. The liquid injection and liq- 
uid withdrawal manifolds are shown bolted directly to the 
test plate in Fig. B.6. The passage for the high-pressure 
air jet that was employed for on-off control of the knife- 
edge capture slot is shown in the withdrawal manifold. The 
locations of the two chromel -al umel thermocouples that were 
employed to measure the inlet and outlet liquid temperatures 
are indicated in the figure. Figure B.6 shows the porous 
stainless steel injector element that was constructed from 
a 1/8-in. thick by 2-in. long piece of "regimesh". The 
cross-sectional view of the knife-edge capture slot is also 
shown in Fig. B.6. Figures B.6 and B.7 show the 2-in, wide 
by 0.060-in. deep trough that was machined into the test 
plate to which the liquid film was confine The 2-in. by 
10-in. wetted test surface was hand polished. Figure B.7 
also shows the two 0.060-in. wide insulating grooves that 
were cut into the underside of the wetted test plate so as 
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to substantially reduce the lateral conduction of heat 
through the test plate to the liquid film. 

The four side walls of the test section were bolted 
together with 3/8-24 stainless steel bolts. Viton-B cord 
was employed to seal the test section along each of the four 
edges. Two flange pieces were welded to each of the four 
test section walls so that a continuous flange surface ex- 
isted at each end of the test section when the test section 
was assembled. The upstream approach section and the 
transition section between the test section and the exhaust 
nozzle were bolted to those flanged ends of the test section 
with Sg- 20 bolts. Extreme care was taken to align the in- 
ternal flow surfaces for the test section with those for the 
approach section. 

Figure B.8 presents two photographs of the assembled 
test section. Figure B . 8 ( a ) shows the configuration for 
that part of the investigation where the gas stream velocity 
was constant. The figure shows the location of the 35 mm 
camera that was employed to photograph the surface of the 
liquid film. Figure B.8(b) shows the test section configu- 
ration for the studies with an accelerated gas stream. For 
those experimental tests, the Pyrex window in the upper wall 
Of the test section v;:: removed so as to have access to the 
manometer leads employed to measure the static pressure 
along the aluminum wedge. The tubing for one of the 5 
static pressure taps is shown in the figure. 
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B . 4 The Separator 

Figure B.9 shows a sectional drawing of the cyclone 
separator that was employed to separate the liquid and gas 
phases. The gas and liquid mixture that was withdrawn from 
the test section through the capture slot entered the sepa- 
rator at top of the main chamber in a tangential direction 
through a Jj-in. O.D. stainless steel tube. The centrifugal 
forces on the liquid drops caused them to form a film on 
the outer wall of the separator and to eventually drain to 
the lower chamber. The gas exhausted to the surroundings 
through the center chamber. The rate at which gas was 
vented to the surroundings (and, therefore, the rate at 
which it was withdrawn at the capture slot) was controlled 
by two precision Hoke needle valves that were mounted on the 
center chamber. 

Figure B.10 shows a photograph of the separator and the 
support stand. Also shown in the figure is the graduated 
beaker employed to collect the liquid sample from the sepa- 
rator and the electric timer that measured the sample time. 
The on-off control for the capture slot (an electric sole- 
noid with a 500 psia air supply) is shown located under the 
test section. The aforementioned timer and the on-off 
control were simultaneously regulated by one of the switches 
that is shown on the front panel of the support stand for 
the separator. The procedure for operating the cyclone sep- 
arator was as follows. First the high-pressure air solenoid 
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was activated so that no gas or liquid entered the capture 
slot. Then the bottom drain valve on the separator was 
opened so that any liquid that remained in the separator was 
removed. The drain valve was closed. Then the timer was 
started and the high-pressure air jet was shut off simulta- 
neously, thereby beginning a mass transfer measurement. 

After a certain length of time, normally one minute, the 
switch was thrown to both stop the electric timer and to re- 
activate the electric solenoid for the high-pressure air 
jet. The drain valve on the bottom of the separator was 
then opened and the liquid in the lower chamber of the sep- 
arator was collected in the graduated beaker. The size of 
the liquid sample that was collected ranged from approxi- 
mately 100 to 800 ml/min, depending on the rate at which the 
liquid was injected onto the test plate. The temperature of 
the liquid in the beaker was frequently measured with a 
mercury thermometer so as to determine the specific gravity 
of the liquid at the time of collection. Knowing the spe- 
cific gravity for the liquid and the volume of the liquid 
that was collected, the average rate of liquid flow out of 
the test section during the sample time could be determined. 



270 


APPENDIX C 

EXPERIMENTAL PROCEDURE 
C . 1 Pre-Run Preparation 

The pre-run preparation of the experimental apparatus 
encompassed three major items, namely, (1) the preparation 
of the gas generator and instrumentation systems, (2) the 
calibration of the instruments, and (3) the establishment of 
the start-up configuration for the entire apparatus prior 
to the start of an experimental run. 

The air storage tanks having a total volume of 2500 cu. 
ft. were charged with air to a pressure between 1500 and 
2400 psi which was sufficient for operation of the gas 
generator for a period of from 2000 to 4000 seconds. The 
fuel tanks were filled with methanol. Also the film coolant 
tank was filled. Several hours prior to an experiment the 
instrumentation was turned on to allow for warm-up and sta- 
bilization. The camera. If employed, was loaded with film 
and mounted on the camera support. The strobe unit was 
checked to ensure proper 'functioning and was patched into 
the synchronization system of the camera. The combined sys- 
tem was then checked again to ensure correct functioning. 
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The pressure transducers monitoring the operation of 
the gas generator were calibrated in 5 psi increments around 
their normally recorded values employing a precision Bourden 
tube pressure gage as a standard. The precision gage was 
periodically calibrated employing an Ashcroft dead weight 
tester. The calibration of the turbine flowmeters monitor- 
ing the air, fuel, and liquid flow rates was accomplished by 
means of a 100 cps reference signal available at the Potter 
flowmeter electronic tachometers. The tachometers supplied 
a d.c. analog output proportional to the rotative speed of 
the flowmeter elements. All of the aforementioned pre-run 
calibrations were recorded on the same stripchart recorders 
which recorded the parameters during an experiment. The 
Pottermeter that measured the rate of liquid flow onto the 
test plate was normally calibrated before each run by the 
time-weight technique. 

The establishment of the start-up configuration for the 
entire apparatus was performed just prior to Ignition. The 
stripchart recorders, which were all set to run at 4 Inches 
per minute, were synchronized. The propellant air and fuel 
and liquid tank pressures were set to the desired values. 

The propellant air and fuel were then bled up to the main 
control valves. The flow rates desired at Ignition were 
established by opening one valve at a time and setting the 
desired flow rate of either fuel, primary air or secondary 
air. Five seconds prior to Ignition all of the stripchart 
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recorders were start 2 d and ran for the duration of the 
experiment. 

The start-up of the gas generator was accomplished with 
the aforementioned automatic start sequencer. 

C. 2 In-Run Procedure 

After mainstage ignition was secured, the operator 
brought the gas generator up to the nominal operating point. 
A period of time from five to ten minutes was required for 
the apparatus to reach thermal equilibrium after the operat- 
ing point had been reached. During the warm-up period the 
film-coolant flow rate was started and the heat exchanger 
for preheating the liquid was turned on. The film of liq- 
uid on the test plate was visually inspected to ensure 
proper coverage. After thermal equilibrium was established 
and the desired operating point was set, the experimental 
measurements were begun. The mass transfer measurements 
followed the procedure outlined in Section 2.5 and Appendix 
B.4. The following tasks were performed during each mass 
transfer measurement: (a) three photographs of the liquid 

film were taken (if the camera was employed during the 
test), (b) the static pressure in the approach section to 

the test section (that was displayed on the precision 

« 

Bourdon tube gage In the control room) was recorded, (c) 
the readings on the 5 mercury manometers for the static 
pressure taps In the aluminum wedge were recorded (If 
applicable), and (d) the inlet and outlet temperatures fof* 
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the liquid, that were displayed on their respective strip- 
chart recorders, were noted and the reostat for controlling 
the inlet temperature was adjusted accordingly. 

C . 3 Post-Run Procedure 

Upon the completion of an experimental test the gas 
generator was shut down. An afterbleed of propellant air 
was allowed to flow through the gas generator and test 
section to reduce the temperature of the equipment. The 
electric heat exchanger for preheating the film coolant was 
turned off. The liquid continued to flow onto the test 
plate during the cool-down period. 


\ 
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APPENDIX D 

PHYSICAL PROPERTIES OF FILM COOLANTS 

D . 1 Subjec t In vest i gati on 
Figures D.1-D.4 present the surface tension a, the 
liquid viscosity y-j, the heat of vaporization H ^ , and the 
vapor pressure p v as a function of the liquid temperature 
T-j. A separate plot is given for each of the film coolants 
of the subject investigation, namely, methanol, butanol, 
water, and RP-1. Because the film coolant RP-1 is a mixture 
of complex hydrocarbons, the properties for this •coolant can 
be specified only approximately. The property values that 
are given herein for RP-1 assume, where necessary, that none 
of the vol ati 1 es • have been removed. The circle that is 
presented for each curve in each of the figures represents 
the maximum withdrawal liquid temperature that was measured 
for that coolant during the course of the investigation. 

Each of the curves that is presented in Figs. D.1-D.4 was 
obtained from one or more of the references (56-60). 

The molecular of the subject film cuuiants 

(employed in the calculation of the maximum liquid tempera- 
ture T, --see Subsection 3.4.1) are as follows: metha- 
1 ,m 

nol : 32; butanol : 74; water : 18; RP-1 : 172. 
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Fig. D.2. Liquid Viscosity as a Function of Liquid 
Temperature for Liquids of Present 
Investigati on 
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Fig. D.3. Heat of Vaporization as a Function of Liquid 
Temperature for Liquids of Present 
Investigation 








D . 2 Investigation Due to Emmons 
Figures D.5-D.7 present the surface tension a, the heat 
of vaporization H y , and the vapor pressure p y as a function 
of the liquid temperature T-j for each of the film coolants 
that was employed by Emmons (20), namely, ethanol, water, 
ammonia, and Freon-113. The values for the calculated maxi- 
mum liquid temperature, T-j m> and for the liquid reference 
temperature, T, , are indicated for the coolants in each of 

I » r 

the figures. Figure D.5 shows considerable extrapolation 
of the data for the surface tension a. The extrapolation of 
the data for a was guided by the critical temperatures for 
the coolants; i.e., at the critical temperature the surface 
tension for the liquid is zero. The values for the critical 
temperature for the various coolants is given in Fig. D.5. 
Each of the solid curves that is presented in Figs. D.5-D.7 
was obtained from one or more of the references (56, 58-61). 
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Fig. D.6. Heat of Vaporization as a Function of Liquid 
Temperature for the Liquids Employed by 
Emmons 
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Fig. D.7. Vapor Pressure as a Function of Liquid 
Temperature for the Liquids Employed by 
Emmons 
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APPENDIX E 

A SAMPLE CALCULATION FOR THE MAXIMUM LIQUID TEMPERATURE 

To illustrate the manner in which the calculated re- 
sults presented in Fig. 34 and in Fig. 35 were obtained, a 
sample calculation for the maximum liquid temperature Tj m 
for an air-water system (Fig. 34) is presented below. 

For the flow conditions represented by Fig. 34, Eq. 

( 3- 76b ) reduces to the form 





=1 

M v 


± 2 ± 


V,$ 


« 1.61 (E-l) 

p v,$ 

Solving Eq. (E-l) for the partial pressure of the vapor at 
the Interface, one obtains 


P 


v,s 


1 

I7ST 



If 

IT 




+ 1 


(E-2) 
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For purposes of illustration, assume that 

T = 2000 R (c = 0.2772 Btu/lb-R) 
9 p air 

p = 500 psia 

For that case, Eq. ( E- 2 ) reduces to 


500 . , _ . 

p v , s ' 1724 H psia E " 4 

v jin ,t_ * 

(JOW-'T-J 1 

To initiate the trial-and-error calculations, assume that 

Ti m = 400F = 860R (the accuracy of the first approximation 

for T i « Is not critical). For T, M « 960R, the data in 
I l »tn 

Ref. 41 gives H » 826.0 Btu/lb and p = 247.3 psia. 

V f in V f s 

Substituting those data into Eq. (E-3) one obtains 
247.3 ps 1 9 l - r re°V l " si » 
f 190.5 psia 

For the second calculation, assume that p„ e « 190.5 psia. 

V t S 

Corresponding to p u ■ 190.5 psia, the data In Ref. 41 

gives H u * 846.7 Btu/lb and T, _ * 837. 6R. Substituting 
v f n» i f w 

those data Into Eq. (E-3) one obtains 

190. S psL 2 i P s, » 

f 189.9 psia 



Continuing this procedure, we assume for the next calcula- 
tion that p = 189.9 psia. Corresponding to p = 189.9 

psia, the data in Ref. 41 gives H y m = 846.9 Btu/lb and 

T, = 377. 4F = 837. 4R. Substituting those data into Eq. 

1 ,m 

(E-3), one obtains 

189.9 = 37633 + 1 pSia 

= 189.9 psia 

Thus, for T] m = 2000R and p = 500 psia, and for the flow 
conditions represented by Fig. 34, the calculated value of 
the maximum liquid temperature is T-j m = 377. 4F = 837. 4R. 


2 

1 
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APPENDIX F 

DATA REDUCTION 
F . 1 General Data Reduction 

The pressure transducers that were employed to measure 
the static and total pressures in the gas generator system 
and the test section were calibrated against a precision 
pressure gage prior to each experiment. The calibration was 
in the form of discrete pressure steps both ascending and 
descending; the maximum value of the largest calibration 
step was in excess of the value recorded during an experi- 
ment. The calibration and the corresponding run data were 
both recorded on the same Brown stripchart recorder. The 
numerical value of the particular pressure parameter re- 
corded during an experiment was determined by linear inter- 
polation of the recorded transducer output against the cali- 
bration. The estimated error resulting from the transducer 
hysteresis and data transcription was estimated at ± 2 per- 
cent at the nominal run value of pressure. 

All of the temperatures monitored during an ■experiment 
were recorded either on direct reading Brown recorders or 
on pre-cal ibrated Brown recorders. The values of tempera- 
ture recorded during an experiment were correspond ngly 
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read off directly from the recorder chart or determined with 
aid of a calibration curve. The estimated error in the 
temperatures read from the Brown recorders was from + 2 to 
+ 4 percent at the nominal run value of the temperature. 

The volume flow rates of air, fuel and film coolant 
were measured with turbine flowmeters. The flowmeters which 
measured the flow rates of the fuel and film coolant were 
calibrated with the respective liquid in terms of frequency 
output as a function of weight flow rate. In this manner a 
calibration constant was established for a particular flow- 
meter through which the weight flow rate and frequency out- 
put were related. The flowmeter employed to meter the flow 
rate of air to the gas generator was calibrated by the com- 
pany supplying the instrument as volume rate of flow versus 
frequency at a static pressure of 600 psig. The flowmeter 
was subsequently calibrated again at the Jet Propulsion 
Center, Purdue University, against ASME standard orifices 
at flow pressures of 400, 600, and 800 psig. The latter 
calibration and the calibration supplied with the meter 
agreed within ± 5 percent at the nominal run value. The 
weight flow rate of air measured during a run was evaluated 
as 




(F-l) 





W « ■ mm* 
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where 

Q = indicated volume flow rate at the frequency out- 
put of the flowmeter 

p = density of the air determined from the measured 
static pressure and temperature of the air 
immediately upstream of the flowmeter, and em- 
ploying the equation of state for a perfect gas 

The output circuit of the flowmeters was calibrated 
prior to each experiment against known frequencies available 
in the flowmeter instrumentation. The estimated error in 
the measurement of the film coolant flow rates was + 2 per- 
cent or less of the nominal run value. The estimated error 
in the measurement of the air flow rate was taken as the 
deviation of the two aforementioned calibrations, that is, 

± 5 percent. 


F. 2 Gas Stream Parameters 

The procedure that was employed to compute the values 
for the gas stream parameters G, u^, and M q at the leading 
edge of the experimental test plate is outlined as follows: 
(a) A first estimate for the gas mass velocity, G, was 
obtained by assuming one-dimensional flow in the experi- 
mental tunnel, in which case it follows from the conserva- 
tion of mass that (neglecting the fuel flow rate into the 
gas generator) 


W 


(F-2) 
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where 

W = measured air flow rate 

A cr = cross-sectional area of the experimental tunnel 

(b) The displacement thickness for the turbulent boundary 
layer on each side wall of the test section at the leading 
edge of the wetted test plate was computed from the classi- 
cal expression (36) 



0.048 x. Re" 0,2 
x x x 


( F-3 ) 


where 



p = viscosity for the gas stream evaluated at the 
® measured gas stream temperature 

= 40 inches 

(c) A corrected cross-sectional area for the flow at x = x^, 
denoted by A', was computed from 

V » 


A' 

A cr 


= A cr - 6,* (wetted perimeter of tunnel) 


(F-4) 


(d) The final value for the gas mass velocity was computed 
by multiplying the result of step (a) by the ratio of the 
uncorrected to the corrected cross-sectional area for the 
tunnel at x 5 x,. 

(e) The gas stream velocity, u g , was computed from the 
result of step (d) for the mass velocity 6 together with 
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the expression 



where the gas stream density was computed from the equa- 
tion of state for a perfect gas and the measured values for 
the gas stream pressure and temperature. 

(f) The momentum parameter, M Q , was computed from 

M rt = P„ u ? (F-6) 

o K g g ' ' 

The accuracy of the computed values for the foregoing 

parameters is estimated as + 5 percent for u and G, and + 5 

y 

to 10 percent for M q . 

For those tests in which the gas stream was not accel- 
erated over the wetted test plate, the values for the pa- 
rameters Ug, G, and M Q were assumed constant throughout the 
test section. For those tests where the gas stream was 
accelerated over the wetted test plate, the local values 

for the gas stream parameters u„, G, and M A for x > x. were 

g o l 

computed from the static pressure measurements that were 
made down the length of the test section. Those calcula- 
tions and measurements re discusseu previou c1 y in Sub- 
section 3.3.1 and will not be repeated here. The accuracy 
of the values that were computed for those local flow pa- 
rameters is estimated as ± 5 percent for u^ and G, and 
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± 5 to 10 percent for M . 

F . 3 Mass Transfer Data 

The procedure that was employed to measure the net rate 
of mass transfer from the liquid film to the gas stream has 
been described in detail in Section 2.5. Moreover, the 
operation of the separator that was utilized to measure the 
rate of liquid withdrawal from the wetted test plate has 
been, adequately described in Appendix B.4. 

During the course of any one experimental test, a 
number of consecutive measurements of the rate of liquid 
withdrawal was normally made for one or more fixed rates of 
liquid injection in order to check the reproducibility of 
the data. The measured rates of liquid withdrawal could 
always be reproduced with an accuracy of 5 percent or better. 

A concerted effort was not made to try to reproduce the 
experimental data of any one experimental test by conducting 
an independent experimental test at the same flow conditions. 
It was very difficult to operate the hot gas generator at 
the same run point for two different experimental tests. 

The primary problem was that of establishing the same rate 
of air flow (and, therefore, gas mass velocity or gas stream 
Reynolds number) on two different occasions. However, a 
matching of the run conditions was realized for two differ- 
ent sets of experimental tests, and in each case the measure- 
ments of mass transfer were found to reproduce to within an 
accuracy of 10 percent or better. 
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APPENDIX G 
EXPERIMENTAL DATA 

Table G.l lists the measured or computed values for 
the pertinent flow parameters for the experimental tests 
with a constant velocity gas stream. The letter that fol- 
lows the test number indicates the film coolant for that 
test; i.e., M - methanol; B - butanol; W - water; and 
R - RP-1. The parameter values listed in the table are the 
average of the values that were measured for each of the 
datum points for mass transfer during an experimental test. 
From 6 to 15 datum points were obtained for each test. Not 
listed in Table G.l are those experimental tests that were 
concerned only with photographic data, those tests which 
were purely exploratory as regards experimental technique, 
etc., or those tests for which some failure of the experi- 
mental apparatus voided the test. 

Table G.2 presents a similar listing of the average 
values of the pertinent flow parameters for the experi- 
mental tests with an accelerating gas stream. 

Figures G.l through G.13 present graphically the mass 
transfer data for the experimental tests listed in Table 
G.l. Figure G.14 presents the mass transfer data for the 


experimental tests that are listed in Table G.2. 



Average Values of the Pertinent Flow Parameters During the Experimental 
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OOHQOvOM-O^CT>NCMHM , OlCOhs 

oioiooocoouncMu.ooiooooNm 

ininmminini£x*eo^fco^f*M-u>min 


i i a i i i i i i i i i i i i « 

NCOO^HfOM-lftCMrOM-iniDNCOOCM 

OOPOrOM’M*M*M*HHHHHHHCM{M 


125-M 418 150 41 19 24 133 160 4.14 5.4 0.00691 
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HU‘)M^vDNCOO^cvjincoo>N^'Hvorocr»coHr)Ln , ^> 

couJo^csjCMCococvjr^ocnincsjLnLnrorocor^cof^fHLn 

<NJO>*H<*ir>cT>cocoooc\iir>csjr^ococoo><-oir>m*o«a- 

HHHHOOHHOHHCVlHfOOOOHHHHOH 

ooooooooooooooooooooooo 


ooooooooooooooooooooooo 


Kr<yi\o<sJoocvj^-ooco^cocNjoiHo>ir>a»oocnio^ooo> 


HO»M^O^Oir)HNOCOlOHWfOONOHOO>V£)^ 
rH rH rH rH rH rH CM HHHHCVJHHHC\JHfO rH 


rocvjror^oir>ro^j“COvo^c\jrH«x>vccjocvi^j“VDcnc\j^r 
r-*cococnLn«M'CMOM*comu. ooooM-*a-o>vocOM-corHo 


ooi/)cr>HcorvH^)Lnrs.Onioosoco^)Ocsju 30 Jvr)^a> 

rH rH rH rH H H rH H H rH 


(n<-oo^^o^oco^ 5 inoooo^^ocooco 
inu)ova>cocoa>a>rHOOrHrHcsjcor^coa>criOrHOPO 
iHHHHHHHH CVI CM CNJ C\J CM CM CM CM CM CM CM CO CO CO CO 


po^voPOHLncoHinmo^rsLncnoco^ovDCMooin^'co 

coco^^<ncocsifo^«tf*u->*f**r*a-t£>t£>'« 3 -vorN.iOCOcoco 

HHHHHHHHHHHHHHHrlHHHHrlHH 


mHin^oo>Hfnin^mc 3 >^cncooir>if)(Mroir)H 
OWCOCSJWCOCO^Cn^NC^NSOOONriNN^ 
HCOHCVJ CM ^ rH rH CO rH rH CM CO rH CO 


rHrHCSJrHCM<JOa>^VDCOO<J-^a><X>OOr^CMrHir>rHin 

^ 1 -fN.co^rrHcsicomcMcomcocM^cocM^in^r^cMCM 


^oioiocvjoinfo^ovDcsjvDHro^^H^ocNjn 
oOincor^ir)rHvocoLoooor^CMioinrHo>coooi0r-Hin 
rH rH rH I-H rH CM iHHrlCMrlrHHCM rH rH 


OCslNM’UHn^OKHHHNiOVO^NrsNMCOOO^’ 
inmKNNKNNlOininNNNNNNNlOinNin 
HH HHH HH rH 


^ONOOOOCOHOHMCninHfO<OOU)HU)MOt 

OHOOOOOHOOO^^OOOOOlftfOCSMCH 

w^'^'^fOvj 


XSttoamttfiOttcDcocottcoaBacocccococfftta 
i i i t i « i i i i i i « i i i i i i « i i 
«OO^MfO^U)\Dt^OHNn^OH(NJ( 0 ^in«ON 
NNOOOOOOOHHHHHOOO © CO O O 
rH rH CM CM CM CMCMCMCMCMCMCM^MCMCOCOCO CO CO PO CO CO 



308-R 554 77 172 35 188 


Table 6.1 (Continued) 
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Table G.2 


Average Values of the Pertinent Flow Parameters During 


the 

Experimental 

Tests 

With An 

Accel erati ng 

Gas Stream 

Test 

T 9 ° 

p° 

T i.i 

T l,2 



F 

psia 

F 

F 

lb 

f t-sec 

401-M 

442 

74 

139 

142 

0.0219 

402-M 

434 

52 

140 

126 

0.0168 

403-M 

440 

103 

144 

154 

0.0267 

405-M 

447 

126 

135 

173 

0.0314 

502-W 

444 

103 

118 

195 

0.0102 

503-W 

435 

54 

78 

163 

0.0062 

505-W 

405 

99 

78 

188 

0.0093 

507-W 

480 

125 

71 

213 

0.0132 
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Data for Mass Transfer. T =600F, p=75 psia. Methanol Film Coolant 
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Data for Mass Transfer. 7 =400F, p=150 psia. Butanol Film Coolant 
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